

















THE JOURNAL 


OF 


CHEMICAL PHYSICS 





Votume 14, NumBER 4 


Aprit, 1946 





An Improved Hot Wire Cell for Accurate Measurements of Thermal Conductivities 
of Gases over a Wide Temperature Range 


Results with Air between 87° and 375°K 


Wi.1aM J. TayLor! AND HERRICK L. JOHNSTON 
Department of Chemistry, The Ohio State University, Columbus, Ohio 


(Received December 26, 1945) 


An improved hot wire cell of the potential lead type is described, and theoretical treatments 
given for radiation, end conduction, potential lead conduction, and the ‘‘temperature jump” 
effect. The latter is found to be satisfactorily eliminated by linear extrapolation of reciprocal 
plots of “apparent conductivities” as a function of pressure. Convection effects are absent at 
pressures below about 20 cm of mercury. The emissivity from bright platinum was measured. 
The results confirm the measurements by Milverton. The thermal conductivity of air was 
measured at seventeen temperatures between 87° and 375°K. Our results, which are believed 
accurate to within +0.5 percent, are generally lower than other results which have been reported 
in recent years, for air. Tabulation of our results is given, together with a critical analysis of 
possible sources of error. The advantage of our method over either the “‘thick wire’’ or ‘‘com- 


pensating lead”’ type of cells is pointed out. 





RESEARCH program on the accurate 
measurement of the absolute thermal con- 
ductivities of gases, over a wide temperature 
range, is being undertaken in this laboratory to 
supplement investigations on measurements of 
heat capacity by the velocity of sound method. 
Investigation of the literature reveals that the 
only reliable thermal conductivity data on gases, 
at temperatures other than 0O°C, heretofore 
reported, are those of Archer? and of Sherrat and 
Griffiths? on carbon dioxide; of Gregory and 
Dock‘ on hydrogen; and the older data of 
Eucken® for several gases. The latter measure- 
‘Fellow of the Charles A. Coffin Foundation 1939-40, 
E. I. du Pont Fellow 1940-41. 
*C. T. Archer, Phil. Mag. 19, 901 (1935). 
* G. G. Sherrat and E. Griffiths, Phil. Mag. 27, 68 (1939). 


*H. Gregory and E. H. Dock, Phil. Mag. 25, 129 (1938). 
* (a) A. Eucken, Phys. Zeits. 12, 1101 (1911); (b) 14, 324 
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ments were relative to air, and were made at 
widely spaced intervals. 

Data on the thermal conductivities of gases is 
also of interest in connection with further de- 
velopment of the theory of transport properties 
of gases. The classical theory is well developed 
for spherically symmetric molecules possessing 
translational energy only,*® but very little prog- 
ress has been made on the more difficult problem 
of molecules with internal degrees of freedom. 
The Chapman-Enskog relationship 


K=ency (1) 


between thermal conductivity (K), viscosity (7), 
and heat capacity (C,) is of particular interest. 


6S. Chapman and T. G. Cowling, Mathematical Theory 
of Non-Uniform Gases (Cambridge University Press, New 
York, 1939). 








220 Ws J. 


Other work from this laboratory’ has supplied 
data on viscosity and on heat capacity, covering 
a wide temperature range. It was felt desirable 
to supplement these results with extensive data 
on thermal conductivity. 

Nearly all precise thermal conductivity meas- 
urements on gases have been made with some 
modification of the hot wire cell. The general 
method of operation of these cells is well known. 
In all of thém a wire, fixed along the axis of a 
cylindrical tube, is heated by an electric current. 
The rise in temperature of the wire, as deter- 
mined from its resistance increases, together with 
the power dissipated and the dimensions of the 
wire and tube, permit calculation of the thermal 
conductivity of the gas in the tube. The prin- 
cipal correction is for heat conducted longi- 
tudinally by the wire, since the thermal con- 
ductivity of the metal is much greater than that 
of the gas. Various modifications of the hot wire 
cell represent attempts to deal with this dif- 
ficulty. The principal types that have been used 
may be referred to as: (1) compensating cells; 
(2) the “‘thick-wire”’ cell; and (3) the potential 
lead type of cell. Compensating cells* involve the 
use of two cells that differ only in length, so that 
differential measurements refer, effectively, to 
the central portion of the longer cell. In the 
“thick-wire” cell, developed by Kannuluik and 
Martin,® no attempt is made to minimize end 
conduction by the use of fine wires. Instead, the 
ends of the cell are so constructed as to lead to 
relatively simple boundary conditions, so that 
the problem of the complete heat flow (along the 
wire and through the gas) may be solved. 

In the potential lead type of cell, which has 
been used by Weber,’ very fine potential leads 
are attached to the cell wire so as to give the 

7 Viscosities: (a) H. L. Johnston and K. E. McCloskey, 
J. Phys. Chem. 44, 1038 (1940); (6) H. L. Johnston and 
E. R. Grilly, ibid. 46, 948 (1942). 

Heat capacities: (c) H. L. Johnston and A. T. Chapman, 
J. Am. Chem. Soc. 55, 153 (1933); (d) H. L. Johnston and 
M. K. Walker, ibid. 55, 172 (1933); (e) H. L. Johnston and 
D. H. Dawson, ibid. 55, 2744 (1933); (f) H. L. Johnston 
and C. O. Davis, ibid. 56, 271 (1934); (g) C. O. Davis and 
H. L. Johnston, ibid. 56, 1045 (1934); (h) H. L. Johnston 
and E. A. Long, J. Chem. Phys. 2, 389 (1934); (z) H. L. 
Johnston and M. K. Walker, J. Am. Chem. Soc. 57, 682 
, investigation with compensating ceils is that 
by B. G. Dickens, Proc. Roy. Soc. A143, 517 (1934). 

®*W. G. Kannuluik and L. H. Martin, Proc. Roy. Soc. 


A141, 144 (1933); A144, 496 (1934). 
10S, Weber, Ann. d. Physik 54, 325 (1917). 


TAYLOR AND H. L. 


JOHNSTON 










ee 





iio 





ert LLL 
pees SYLLLL/ 






QIZIILTIIEIEIEOy 










. CLidbhddiidaiiiun 








ENNISES Vg 


\2 
a 





















aN 





SSOOASSS 





New Desion 
THermar Conovctrivity Cer 
A 


Giass 
BRaAss~. CoprPer 


mico Sree 


BS\NZ 


Sot DER 


ONE IncHy 





One INCH , 


Fic. 1.,Thermal conductivity cell. 


potential drop through a central segment. End 
conduction is thus reduced to a small correction. 


io i ee ee ee ee a a a. 


Loman 


THERMAL CONDUCTIVITY 


The characteristics of these three types of cell 
may be summarized by the statement that the 
end ‘conduction is canceled out in the com- 
pensating type of cell; is large but subject to a 
relatively exact correction in the ‘‘thick-wire’”’ 
type of cell; and is essentially absent in the 
potential lead type of cell. 

The potential lead type of cell was chosen for 
the present series of investigations because of the 
simplicity of the calculations in the treatment of 
the data and because the near absence of end- 
conduction removes the possibility of serious 
systematic error. Important improvements were 
made in the design of the cell. 

The present paper is limited to a description 
of the apparatus and method, and to results 
obtained with air. A subsequent paper by one of 
us will present data obtained with several other 
gases. 

Air used in this investigation was taken from 
outside the laboratory and freed from carbon 
dioxide and water vapor by slowly bubbling it 
through a train consisting of: (a) concentrated 
sodium hydroxide ; (b) concentrated sulfuric acid ; 
(c) stick potassium hydroxide; and (d) phos- 
phorus pentoxide. Finally, it was passed through 
a liquid nitrogen trap (at less than atmospheric 
pressure so that no air condensed). 


DESCRIPTION OF APPARATUS 


The apparatus consists of : (1) the cell; (2) the 
thermostated bath ; and (3) the electrical circuits. 
Their description follows. 


(1) The Cell 


An accurately bored steel tube, with a diam- 
eter sufficiently small to insure the absence of 


convection, was used for the cell. Positive 
centering of the 0.01-inch platinum cell wire was 
obtained by the use of cylindrical pistons sliding 
in the tube. Constant tension on the cell wire 
was maintained by means of a weight, so chosen 
that there was negligible vibration of the wire. 
Platinum potential leads, 0.0005 inch in diam- 
eter, were used. A more detailed description 
follows. 

The construction of the cell is illustrated in 
Fig. 1, where it should be noted that the vertical 
scale of the main drawing is compressed in the 
ratio 1:2 with respect to the horizontal scale. 
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The full length inset of the wire system is on 
the same scale as the cell, while the inset showing 
the cross-sectional view of the tube interior and 
wires is on an enlarged scale. 

The central portion of the cell consists of a 
mild steel tube, (a), 11 inches in length, ? inch 
in outside diameter, and approximately 0.800 
cm in internal diameter, which was accurately 
bored by the Winchester Repeating Arms Com- 
pany, New Haven, Connecticut. The exterior of 
the tube was chromium plated, to prevent cor- 
rosion. At the upper end this tube is soft soldered 
to a copper-Pyrex glass seal, (b), about 2.5 cm in 
diameter, and at the lower end to a long brass 
cap, (c), forming an extension of the tube. 

The cell wire was selected for freedom from 
irregularities from Bishop C. P. 0.01-inch 
platinum wire. A 22-cm length, (d), was fused 
to 1-mm platinum leads (e, e), which, in turn, 
were silver soldered through holes drilled cen- 
trally in small steel cylinders, (f, f), 0.43 cm in 
diameter by 1.5 cm in length. The cell wire was 
brought to the axis of these cylinders by slightly 
bending the heavy platinum wires until observa- 
tion with a micrometer comparator showed 
proper centering. Glass tubes (g, g), 1.5 cm long, 
fit closely over the steel cylinders (f, f), and 
fit closely within the steel tube. These tubes 
were specially moulded by the Fischer-Schurman 
Company, and were selected for concentricity. 
The glass served to insulate the cell wire, elec- 
trically, from the tube. In order to permit the 
free passage of gas, a 60° section was cut from 
each of the glass insulators. This use of glass 
insulated ‘‘pistons’” insures accurate centering 
of the cell wire within the tube and at the same 
time permits vertical motion made necessary to 
compensate for thermal expansion or contraction 
of either the upper lead wire or of the cell wire. 
The possible departure of the cell wire from 
center is estimated as not more than 0.014 cm, 
which is about the radius of the wire. 

A constant tension is maintained on the cell 
wire by means of an 18-gram weight, (h), sus- 
pended from the lower “‘piston’’ but insulated 
from it by a glass bead, (z). It is important to 
maintain constant tension on the cell wire since 
it is used as a sensitive resistance thermometer. 
It was determined that, with this amount of 
loading, vibration of the cell wire was negligible 
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under normal operating conditions (maximum 
amplitude 0.001 cm). 

The current lead to the upper end of the cell 
wire enters by way of a tungsten through Pyrex 
seal. The external portion of this lead is insulated 
from the bath by a glass tube fused to the seal. 
The lower end of the cell wire is connected to the 
steel tube by means of a short spiral, (j), of 
0.01-inch copper wire of small rigidity. The 
lower external current connection is made di- 
rectly to the metal cell. 

The system of potential leads (Rk, k) enters the 
cell by tungsten-through-Pyrex seals at the 
upper end of the copper-Pyrex seal. Like the 
upper current leads, the external portions of the 
potential leads are insulated from the bath by 
glass tubes. The external leads and the portions 
of these leads that extend between the tungsten 
to Pyrex seals and the upper piston are of 
0.01-inch platinum. The 0.01-inch platinum 
leads pass downward through 0.02-inch holes in 
the upper piston, from which they are insulated 
by means of glass capillaries. Between the two 
pistons the leads continue as 0.003-inch platinum 
wires parallel to and 0.3-cm distant from the 
cell wire. These 0.003-inch wires are held under 
a 2-gram tension by means of springs (/, /) con- 
sisting of 0.15-cm helices of 0.003-inch piano 
wire, hooked to the lower piston but insulated 
from the platinum leads by glass beads (m, m). 

From the 0.003-inch platinum lead wires con- 
nections are made to the central portion of the 
cell wire—6.3 cm in length—by 0.5-cm lengths 
of 0.0005-inch diameter platinum Wollaston 
wire (n, m), the junctions. being gold soldered. 
The junctions with the cell wire each required 
only 0.015 mg of gold (a 1X5-mm strip of gold 
leaf) but inspection with a 20 microscope 
showed that the Wollaston wires emerged cleanly 
from the gold and normal to the cell wire. 

In order to prevent injury to the Wollaston 
wires by an accidental movement of the parallel 
“side wires’ it was found necessary to add the 
two semi-circular braces (0,0), of 0.003-inch 
platinum wire shown in the inset of Fig. 1. 
These braces are 1 cm distant from the Wollaston 
wires and are at no point nearer than about 
0.25 cm to the cell wire. They are insulated from 
the ‘“‘side wires” by 0.5-mm glass beads. 

Tests with the completed cell showed that no 
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TABLE I. Dimensions (at 300°K) that enter directly into 
the calculations of thermal conductivity. 








Internal radius of tube (re) 
Radius of cell wire (r1) 
Effective length of cell wire (L) =6.327 cm 


=0.01277 cm (+0.18%) 
(+0.03%) 








anomalous effects were caused by the use of the 
fine Wollaston wires, even though the voltages 
were measured with a precision of 0.1 microvolt 
(or 1 part per million). It was determined, by 
reversing the current in the cell, that no per- 
ceptible thermoelectric force was produced in the 
potential leads by heating of the cell wire. 

Dimensions of the cell which enter directly 
into the accurate calculation of the conductivity 
of the gas are the internal radius of the cell tube, 
the radius of the cell wire, and the effective length 
of the cell wire (that is, the length between the 
potential leads). The radius of the tube was cal- 
culated from its length and volume; the length 
was measured on a Société Genevoise micrometer 
comparator and the volume was determined by 
weighing the mercury required to fill the tube 
(using plane end pieces). The radius of the 
platinum cell wire was determined by weighing 
the segment used in the cell. The length of the 
segment was measured with a Gaertner cathe- 
tometer, after annealing the wire for eight days 
at 600°C under a tension of 20 grams. The value 
selected for the density of platinum, which was 
required in the calculation, was 21.40+0.06 
grams cm~ at 20°C." 

The effective length of the cell wire was 
measured with a Gaertner micrometer com- 
parator after the potential leads were attached. 
The final values obtained for the principal dimen- 
sions of the cell at 300°K, and the estimated 
uncertainties, are given in Table I. 

Dimensions that are required only for cal- 
culating corrections to the conductivity were 


1! Kahlbaum and Sturm, Zeits. f. anorg. Chemie 46, 217 
(1905), report a density of 21.44 grams cm~* for pure 
annealed platinum wire; Stas, Oeuvres Complétes, Brux- 
elles-Leipzig 2, 742 (1894) found platinum to approach an 
upper limiting density of 21.46 grams cm= on working 
and annealing. The widely quoted value of 21.37 g cm™ 
is due to a misinterpretation of the work of Richards and 
Stull, Carnegie Inst. Pub. 76, 55 (1907), who reported a 
density of 21.31 grams cm™ but also tabulated the specific 
gravity referred to 20°C, namely 21.37. Since Richards 
and Stull worked with a platinum cylinder it might be 
expected that their result would be lower than the maxi- 
mum densities quoted above. 


=0.40017 cm (+0.02%) 
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determined with less accuracy. They are listed 
in Table II, together with the estimated uncer- 
tainties. 

The lengths of the end segments of the cell 
wire differ by only 0.04 cm and it is permissible 
to consider them equal, and equal to (/)—/)/2. 


(2) Thermostat 


In the completed experimental arrangement, 
the cell is mounted in a thermostat, containing 
a suitable bath liquid, which may be adjusted to 
maintain any temperature in the range 120° to 
400°K. When the thermostat is carefully ad- 
justed the oscillations of the temperature about 
its mean value may be reduced to +0.003°. The 
mean temperature may vary on the order of 
0.01° per hour, but these slow variations in- 
troduce no error into the measurements. The 
thermostat is similar to the cryostat described 
by Scott and Brickwedde,” but the smaller Dewar 
tube available made necessary a different place- 
ment of the cooling and heating surfaces. An 
improvement is the use of a _ continuously 
variable heater control circuit which decreases 
the tendency of the temperature to oscillate 
about its mean value. 

A Dewar tube 4 inches in diameter and 36 
inches deep holds the bath liquid. It stands 
beneath a circular hole cut in a heavy (400- 
pound) stone table top, to which the support for 
the thermal conductivity cell is bolted; the cell 
makes no direct contact with the other objects 
in the thermostat. Sections of red-fiber were 
introduced into all metal supports descending 
into the bath to minimize heat leak, and the 
cell is further insulated from its support by red- 
fiber washers. 

The heavy (20-pound) brass stirring well, 2 
inches in internal diameter, is placed at one side 
of the Dewar. A ?-inch stainless steel shaft, 
turning in brass bearings mounted on the outside 


TABLE II. Dimensions (at 300°K) that enter only into 
correction terms. 








Total length of cell wire (Jo) =22.09 cm (+0.2%) 
xternal radius of tube (73) = 0.795 cm (4+1%) 


Radius of potential leads 
6.25X10-* cm (+10%) 


(Wollaston) (rz) 
ngth of potential leads (Jz) = 0.5 cm (410%) 








" R. B. Scott and F. G. Brickwedde, J. Research Nat. 
Bur. Stand. 6, 401 (1931). 
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of the stirring well, drives a 23-inch propeller at 
‘tthe bottom of the Dewar tube. The shaft is 
coupled by means of a red-fiber rod to the drive 
shaft above, which is equipped with oversized 
ball bearings. With the cooling tube (to be 
described) in place, the bath liquid is circulated 
(up through the stirring well) two to four times 
per minute at 400 r.p.m. 

The cooling tube, an unsilvered Pyrex Dewar 
tube 13 inches in external diameter and 30 inches 
deep, with a separation of 1 mm between the 


- walls, is mounted so that the lower 24 inches 


descends into the stirring well. When cooling is 
required, liquid oxygen (or nitrogen) is main- 
tained at a constant level in the cooling tube by 
a float and relay arrangement which controls the 
flow of liquid oxygen from a reservoir. Control of 
the cooling rate is effected by varying the gas 
pressure between the double walls of the cooling 
tube. A pressure of 0.5-mm Hg of hydrogen gives 
sufficient cooling at 125°K while it is possible to 
thermostat a few degrees below room tempera- 
ture with liquid oxygen in the Dewar by reducing 
the pressure to 0.002-mm Hg. With 1-cm Hg of 
hydrogen in the tube, the bath is cooled at an 
initial rate of 50 degrees per hour; the rate drops 
to 20 degrees per hour at 160°K and to 5 degrees 
per hour at 120°K. 

The heater consists of 30 ohms of bare 
Nichrome wire wound bifilarly directly on the 
part of the cooling tube in the bath. The bath 
liquid thus circulates through the annular space 
between the cooling tube and the stirring well, 
in intimate and simultaneous contact with the 
cooling and heating surfaces. It is convenient to 
adjust the pressure in the cooling tube to give 
somewhat greater than optimum rate of cooling 
at the desired temperature, and then to reduce 
the cooling rate to the small optimum value by 
supplying a constant heating current to the 
heater. Similarly, in thermostating above room 
temperature, the constant heating current is 
adjusted to give a small rate of cooling at the 
desired temperature. 

The actual thermostating, in both cases, is 
then accomplished by automatic regulation of a 
few watts of variable heating. A six-junction 
copper-constantan thermocouple, with the junc- 
tions exposed directly to the bath, is balanced 
against a Leeds and Northrup Type K poten- 
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tiometer. A beam of light is reflected by a Leeds 
and Northrup Type R galvanometer onto a 
gas-filled photoelectric cell, which in turn con- 
trols the variable component of the heating 
through a phase shift Thyratron circuit, as 
described by Strong." 

Several bath liquids were employed. Above 
room temperature a low viscosity mineral oil was 
used ; water was not satisfactory as it was ab- 
sorbed by the red-fiber. Difluorodichloromethane 
was used from a degree or two above its melting 
point, 120°K, to about 160°K. This bath gave 
excellent thermostating because of its low vis- 
cosity. A mixture of 32.7 percent ethyl bromide, 
26.5 percent methylene chloride, 15.6 percent 
trichloroethylene, 13.6 percent chloroform, and 
11.6 percent transdichlorethylene (by volume), 
as described by Scott and Brickwedde,” was 
satisfactory in the range 150° to 250°K. In the 
range 205° to 300°K, a mixture of equal parts by 
weight of chloroform and carbon tetrachloride 
was employed. 

Temperatures near 80° and 90°K were ob- 
tained by filling the thermostat with liquid 
nitrogen or liquid oxygen, respectively. It was 
necessary to stir these baths vigorously to 
maintain a steady temperature. 


(3) Electrical Circuits 


The electrical circuit used to measure both the 
resistance of the ‘effective length” of the cell 
wire, and its energy dissipation, is quite simple. 
Its resistance—which ranged from about 0.029 
ohm at 80°K to about 0.177 ohm at 380°K— 
was determined by dividing the potential dif- 
ference between the potential leads by the 
current through the wire. The latter was deter- 
mined by measuring the potential difference 
across a standardized manganin resistance of 
0.128477+0.000025 ohm, at 25°C, which was 
inserted in the external circuit of the current 
leads. Potential differences were measured with 
a White double potentiometer of 100,000 micro- 
volt range, with a precision of 0.1 microvolt. The 
current through the cell, and standard resistance, 


13 J. Strong, Procedures in Experimental Physics (Pren- 
tice-Hall, Inc., New York, 1938), p. 447. Note that the 
connection shown by Strong between the cathode of the 
photoelectric cell and the plate of the Thyratron should 
be eliminated. 
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was supplied by two 2-volt heavy duty lead 
storage cells, of 450-ampere hour capacity, in 
parallel, and amounted to about 0.1 ampere 
during (resistance thermometer) calibration of 
the cell wire and to about 0.5 ampere during 
thermal conductivity measurements. The current 
was varied in steps, as desired, by a set of 
manganin resistors, with mercury contacts, 
which were in series with the cell wire and were 
adjustable in steps of 0.01 ohm. 

Since it is important that the standardized 
resistance in series with the cell wire be prac- 
tically independent of the cell current special 
attention was given to its design and its con- 
struction. It consists of 15 parallel strands of 
bare manganin wire (0.25 ohm per cm) soldered 
between copper bars (5X 10~-* ohm per cm). The 
entire unit is immersed in yellow ceresin wax. It 
was calibrated by comparison with National 
Bureau of Standards certified resistances of 
nominal 10 ohms, 100 ohms, and 10,000 ohms, 
respectively. A ratio method was employed, and 
potential differences were measured on the White 
potentiometer. The final value of this calibration 
is given above, its uncertainty being fixed by the 
small uncertainties assigned to the National 
Bureau of Standards certification of the standard 
resistances. The resistance of the standardized re- 
sistance proved constant, to within 0.01 percent 
for currents ranging from 0.06 to 0.6 ampere. 

Intercomparisons of the decade coils of the 
White potentiometer were made in the manner 
recommended by the manufacturer (Leeds and 
Northrup) and the resistance ratios were found 
to lie within permissible tolerances. In the course 
of the measurements, small contact and thermal 
electromotive forces in the potentiometer circuit 
and in external circuits were determined and the 
proper corrections applied. 


TEMPERATURE MEASUREMENTS 
(1) Temperature Scale 


The temperature scale for this research was 
established by means of a standard thermo- 
couple, thermocouple C, which was mounted in 
the thermostat. The junction was about 1 cm 
from the central portion of the steel cell tube. 
Thermocouple C was calibrated against Univer- 
sity of California thermocouple No. 101, which 
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had been calibrated by Greensfelder and Milner'* 
in terms of the hydrogen scale of Giauque, 
Buffington, and Schulze.* The original calibra- 
tion of couple C extended over the range 110° to 
300°K. The correction curve was extrapolated to 
80°K. The calibration was extended to higher 
temperatures by comparison, at 323°, 343°, and 
373°K, with a mercury-in-glass thermometer 
certified by the National Bureau of Standards. 
The three new calibration points lie on a smooth 
extrapolation of the original table. It is believed 
that the temperature scale is reliable to about 
+0.05 degree over the entire range. The electro- 
motive force of the thermocouple is read on a 
Leeds and Northrup, White double potentiom- 
eter, with a range of 10,000 microvolts in steps 
of 1 microvolt, and a galvanometer sensitivity of 
6cm on the scale per microvolt. The sensitivity 
of the thermocouple varies from 20 to 50 micro- 
volts per degree. The precision of measurement 
of the temperature was 0.003 degree, in a series 
of measurements at a given temperature. 


(2) Calibration of Cell Wire 


The calibration of the resistance of the cell 
wire was carried out concurrently with the 
thermal conductivity measurements on air. The 
wire had previously been annealed at 600°C by 
passing a current through the evacuated cell. 
The calibration involved a measurement of the 
resistance of the wire when heated by a current 
of about 0.1 amp. This value was a compromise 
between the requirement for minimum heating 
of the cell wire and the necessity for sensitivity 
in the measurement of the temperature of the 
wire. The resultant heating of the cell wire was 
about 0.50 degree (approximately independent 
of temperature), and the precision of measure- 
ment was 0.003 degree. 

The initial correction for the heating of the 
cell wire was effected by a linear extrapolation 
of the wire resistance against the square of the 
current, to zero current. This was _ possible 
through the use of the conductivity measure- 
ment, made concurrently with about 0.5 ampere. 
A preliminary correction curve for the wire was 
then computed, and preliminary thermal con- 

* (a) B. S. Greensfelder and R. T. Milner, J. Am. Chem. 


Soc. 50, 2205 (1928); (b) W. F. Giauque, R. M. Buffington, 
and W. Schulze, ibid. 49, 2343 (1927). 
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ductivities of air calculated. The latter then per- 
mitted an accurate calculation of the heating of 
the wire in the 0.1 ampere runs, and thus the 
calculation of the final correction curve. The 
final curve differed from the preliminary curve 
by 0.07 degree at the lowest temperature and by 
0.01 degree at the highest, in just the way it was 
anticipated the dependence of the resistance and 
the gas conductivity on the absolute temperature 
would affect the linearity of the aforementioned 
extrapolation. 

The corrected resistance readings were com- 
pared with the temperature, read on the standard 
thermocouple by means of a deviation plot, on 
which deviations of the observed resistances from 
an equation in T were plotted against tempera- 
ture. The average deviation of the 20 calibration 
points from a smooth curve on the deviation 
plot was +0.015 degree. However, the average 
deviation of 17 of the 20 points was only +0.008 
degree. | 


CALCULATION OF THE CONDUCTIVITY 
(1) General 


The heat produced in the effective length of 
the cell wire is dissipated principally by con- 
duction through the gas to the cell wall. Cor- 
rections are required for longitudinal conduction 
in the cell wire, for conduction of heat from the 
cell wire by the potential leads, and for radiation 
from the cell wire. Correction must also be made 
for the temperature drop through the cell wall, 
and for temperature discontinuities between the 
gas and the surfaces of the cell wire and cell wall 
—the well-known ‘‘temperature-jump”’ effect. 

The temperature-jump correction requires 
that measurements be carried out at two or more 
gas pressures. 

The heat flow in the gas in the central region 
of the cell may be regarded as radial—that is, 
the effect of longitudinal heat flow in the gas is 
negligible.'* The differential equation governing 
the heat flow is, therefore 


—0T/dr =(Q—2qi)/2eLKr, (2) 


where T is the temperature at the radial distance 


1% E. H. Kennard, Kinetic Theory of Gases (McGraw-Hill 
Book Company, Inc., New York, 1938), Chap. VIII. 

16 W. G. Kannuluik and L. H. Martin, Proc. Roy. Soc. 
(London), have given the necessary equations for calcu- 
lating the latter effect (see reference 9). 
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r; Q is the total heat produced per unit time in 
the effective length, L, of the wire; g; is the heat 
flowing per unit time from the effective length 
of the cell wire by any mechanism other than 
conduction through the gas; and K is the thermal 
conductivity of the gas at the temperature T. 


(2) Temperature-Jump Correction 


In order to account for an apparent small 
dependence of the conductivity, as measured in 
a hot wire cell, on the pressure of the gas, it is 
found necessary to allow for the temperature- 
jump effect at the surface of the wire and, to a 
lesser extent, at the cell wall. The existence of 
such a temperature discontinuity when heat 
flows across a gas-solid interface is well estab- 
lished, and the temperature jump, 67, is known 
to be proportional to the temperature gradient 
in the gas and inversely proportional to the 
pressure of the gas.'® 


8T =gdT /dr =(g'/p)dT/dr. (3) 


When account is taken of the temperature 
jump, Eq. (2) may be integrated to yield the 
relation 


1/Ka=(1/K)+(A/p), (4) 


where K is the true average conductivity of the 
gas over the range 7; to J2, and Ka, the apparent 
conductivity at the pressure , is calculated from 
the relation 

Ko=K.u—)d: Ki. (S) 


K,, is the uncorrected conductivity of the gas 
calculated from the experimental data and the 
equation 


K.=QI1n (r2/r1)/2aL(T1—T>2). (6) 


K; is an effective conductivity for the 7th sub- 
sidiary mechanism of heat conduction, as defined 
by 

Ki=4: In (r2/r1)/2eL(T,—T>2). (7) 


The constant A in the small correction term in 
Eq. (4) is given by the expression 


A= [1/In (ro/r1) IL (gi’/ Kurs) +(g2’/ Kore) J, (8) 


where K, and K, are the conductivities of the gas 
at 7, and 7», respectively. 

It follows from Eq. (4) that a plot of (1/K.<) 
against (1/p), for constant 7, and T», yields a 
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straight line, and that extrapolation of the line 
to its intercept at (1/p) =0 yields (1/K), where 
K is the true conductivity. It is clear from Eq. 
(3) that the temperature jump effect vanishes at 
(1/p) =0, but actual measurements at very high 
pressures are not practicable because of the 
great increase of convection at high pressures. 
The use of Eq. (4) is thus a method of correcting 
for the temperature-jump effect and at the same 
time avoiding convection in the gas. Careful 
experimental tests of Eq. (4), as applied to a hot 
wire cell, have been made by B. G. Dickens. 
The results of the present investigation also 
confirm the correctness of the linear relation 
predicted by Eq. (4), as will be shown later. 


(3) Correction for End Conduction 


The correction for end conduction, or for the 
heat flowing from the effective length of the cell 
wire into the end-segments of the cell wire, may 
be calculated with considerable accuracy. The 
differential equation for the longitudinal tem- 
perature distribution in the cell wire may be 
written 





a | 30+B,=0, (9) 
dx? 

:-| 3 ~(BRs Ro), (10) 
Ar? In (72/71) 

By=(Rol®/ar:2LyJ), (11) 


where 0 =(7,—T2), at constant T2, and x is the 
longitudinal distance along the wire. The radii 
of the cell wire and of the cell are r; and 1, 
respectively, and L is the length of the effective 
portion of the cell wire. K is the average thermal 
conductivity of the gas, and ) is the thermal con- 
ductivity of the (platinum) cell wire. Ro is the 
resistance of the effective length of the cell wire 
at 0=0, (T,1=T:), and Ry’ is the temperature 
coefficient, (dRo/dT2). I is the current flowing 
through the cell wire, and J is the mechanical 
equivalent of heat. 

The boundary conditions on Eq. (9) are that 
@=0 at the ends of the wire, x=0 and x=Ly, 
where Ly is the total length of the cell wire. The 
nature of the solution depends upon whether the 
constant B, is real or imaginary; B, is real for 
the present cell, under the conditions of opera- 
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tion. Instead of calculating directly the heat 
flowing from the effective length of the cell wire, 
@ may be integrated over the range x = (Iy)—L)/2 
tox=(Lo+L)/2 to give the average temperature, 
@ av., of the effective length, which is related to 
the observed resistance R. The result is: 


0 av. = [(R—Ro)/Ro’ ] 





si (B./B,)|1 ae (2/B,L) 


sinh (B,L/2) 
| (12) 
cosh (B,L/2) 


As Lo increases without limit, the correction for 
end conduction vanishes, and © av. approaches 
the constant value (B,/B,’). The effective con- 
ductivity for end conduction, K,, is therefore: 


K.=((2/B,L) sinh (B,L/2)/ 
cosh (B,Lo/2) |Ku, (13) 


where K, is given by Eq. (6). 

For the present cell, and for J=0.5 ampere, 
K, decreases from 0.0086 X 10- to 0.0006 x 10-> 
cal. deg.-' cm=! sec.-! as K, increases from 
2x10 to 5X10-* cal. deg.—! cm= sec.—!. The 
maximum value of the correction in the present 
investigation was 0.5 percent. 


(4) Lead Wire Correction 


The correction for the heat flowing from the 
cell wire to the platinum potential leads is not 
susceptible to as accurate a calculation as the 
other corrections. The following treatment is a 
modification of that of S. Weber.’® The tem- 
perature distribution in the potential lead, as a 
function of the distance from the cell wire, x, 
should satisfy Eq. (9) approximately, with B,=0 
and 

By,=[2K/X(r1’)? In (r2’/ri’) J}, (14) 


where 7,’ is the radius of the potential lead, and 
ry’ isan effective radius, of the order of magnitude 
of the radius of the tube, 72. The boundary con- 
ditions are O=(71—T2) at x=0 and O=0 at 
*=L’, where L’ is the length of the potential 
leads. The heat flowing into the potential leads 
may be calculated from the value of (dOQ/dx) at 
*=0 (only half of this heat contributes to the 
correction). The result may be expressed as the 
effective conductivity of the potential leads, K>, 
by means of Eq. (7), 


K,=(r1’)?B, In (r2/r1)/2L tanh (B,L’). (15) 


It is clear that K, is not sensitive to the value of 
(ro’/ry’). It seems safe to assume that 


10°< (ro’/ry’) <10* 
since these limits correspond to 
0.06 cm <ro’ <6 cm 


and r2’ should be of the order of magnitude of 
r2=0.4 cm. For (r2’/r:’) in this range, and for 
the values of the gas conductivity that occur 


0.99 <tanh (B,L’) <1, 


so that the correction is essentially independent 
of L’. On setting tanh (B,L’) equal to unity, and 
introducing the above mentioned limits for 
(ro’/ry’) and the numerical values of the other 
constants, Eq. (15) reduces to 


K,X10°=(0.012+0.002)(K X 10°). (16) 


The largest value of the correction in the 
measurements on air was 0.80 percent, and the 
average value was 0.50 percent. 


(5) Radiation Correction 


A correction is required for the transfer of 
heat by radiation from the cell wire. Milverton!’ 
has measured the radiation loss from a polished 
platinum wire in a hot-wire cell, over the range 
0° to 100°C. His results are in essential agree- 
ment with the theoretical expression for radia- 
tion between two coaxial cylindrical surfaces." 
A measurement of the radiation loss from the 
present cell was carried out at 317°K. Conduc- 
tion of heat by the residual gas in the evacuated 
cell was negligible. Under the conditions of the 
measurement, Eqs. (10)—(12) applied, with an 
obvious change in the meaning of K, since radi- 
ation replaced thermal conduction. The result of 
the measurement, expressed as the emissivity of 
the platinum wire, was 5.3 percent higher than 
the value found by Milverton at the same 
temperature. 

For the calculation of the radiation correction 
in the thermal conductivity measurements, the 
value of the emissivity of platinum measured at 
317°K in the present cell was adopted, together 


17S. W. Milverton, Phil. Mag. 17, 414 (1934). 

18M. Jakobs and G. A. Hawkins, Elements of Heat 
Transfer (John Wiley and Sons, Inc., New York, 1942), 
pp. 125-126. 
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TABLE III. Thermal conductivity of air at 273.263°K. 








Time (A.M.) 2:08 
Pressure, p (cm Hg) 19.55 
Reciprocal of pressure, (1/p) 0.051 
E.m.f. of thermocouple C (microvolts) 254.2 
Wall temperature, T» °K 266.457 
E.m.f. across 0.128477 ohm re- 

sistor (microvolts) 69437.3 
Cell current, J (ampere) 0.54046 
E.m.f. across cell (microvolts) 70047.5 
Cell resistance, R (ohm) 0.129606 
Wire temperature, 7; (°K) 280.068 
Average temperature, 

Tay =43(T1+T2) (°K) 273.263 
Temperature interval, 

@=(7:-—T:2) (°K) 13.611 
Uncorrected conductivity, (cal./cm 

K,,=0.020720(RI?/0) sec. deg.) 5.763 
Sum of effective conductivities of 

end conduction, potential leads (cal./cm 

and radiation sec. deg) 0.105 
Apparent conductivity, (cal./cm 

Kea (at Tay and p) sec. deg.) 5.658 
Apparent conductivity, (cal./cm 

Ka (at 273.263°K and p) sec. deg.) 5.658 
(1/K,), at pcm Hg) 1.7674 
[1.7650+-0.0496(1/p) ] x 104 1.7675 


2:50 3:26 4:06 4:39 7:46 
8.94 4.37 1.70 0.615 8.30 
0.112 0.229 0.588 1.63 0.120 
254.3 254.2 254.3 254.4 254.1 
266.454 266.457 266.454 266.452 266.460 
69433.6  69429.2 69422.9 69410.3 69435.8 
0.54044 0.54040 0.54035 0.54025 0.54045 
70049.4 70056.5 70084.4 70175.7 70053.8 
0.129616 0.129638 0.129701 0.129894 0.129620 
280.089 280.134 280.265 280.663 280.097 
273.272 273.296 273.360 273.558 273.279 
13.635 13.677 13.811 14.211 13.637 
5.753 5.735 5.681 5.528 5.752 X10> 
0.105 0.105 0.105 0.105 0.105 x 10> 
5.648 5.630 5.576 5.423 5.647 X 1075 
5.648 5.629 5.574 5.418 5.647 X10 
1.7705 1.7765 1.7940 1.8457 1.7709 X10! 
1.7706 1.7764 1.7942 1.8458 1.7710 X10! 








True conductivity, at 273.263°K, K =(1/1.7650 X10‘) =5.666 X10-5 cal. cm-! sec.-! deg.-1. 


with the temperature coefficient found by Mil- 
verton. The effective conductivity for radiation, 
K,, was then calculated from Eq. (7). In the 
thermal conductivity measurements on air, the 
correction amounted to 0.20 percent at 100°K, 
0.74 percent at 200°K, 1.65 percent at 300°K, 
and 2.59 percent at 375°K. 


(6) Wall Correction 


A small correction for the temperature drop 
through the wall of the steel cell tube is readily 
calculated from the internal and external radii 
of the tube, and the thermal conductivity of the 
steel, 0.11 cal. cm™ sec.—! deg.—'. The greatest 
value of the correction in the measurements on 
air was 0.01 percent: 


SAMPLE CALCULATION 


The complete calculation of the thermal con- 
ductivity of air at 273.263°K is given in Table 
III for purposes of illustration. In this particular 
experiment the conductivity of air was measured 
at six pressures, over a period of about six hours. 
The uncorrected conductivity, K.u, at each 
pressure, is first calculated from Eq. (6). The 
effective conductivities for end conduction, K,, 
the potential leads, K,, and radiation, K,, were 
0.0004, 0.0286, and 0.0764 10- cal. cm sec. 


deg.—!, respectively. Subtraction of the sum of 
these corrections from K,, as in Eq. (5), yields 
the apparent conductivity, Ka, at the pressure 
p and the average temperature for the experi- 
ment. The values of K, are then corrected to a 
constant temperature, which is taken here as the 
average temperature of the measurement at the 
highest pressure, 273.263°K. These minute cor- 
rections may be made with an approximate 
temperature coefficient. According to Eq. (4) a 
linear relation should exist between (1/Ka,) and 
(1/p). As shown in the table, the six points show 
a maximum deviation from a linear relation in 
(1/p) of about 0.02 percent. Extrapolation of 
this linear relation to (1/p)=0 yields the true 
thermal conductivity of air at 273.263°K, 
5.666 X 10— cal. cm sec. deg.—. 


EXPERIMENTAL RESULTS 


Results obtained with air, in the temperature 
range 87° to 376°K, are given in Table IV. 
All data were obtained in August 1941 and the 
dates, listed in the first column, indicate the 
order in which measurements were made. Tem- 
peratures, in the second column, are the mean 
of the wire and the cell wall, which usually dif- 
fered by about 13° in temperature. Pressures, 
in the third column, are those at which apparent 
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conductivities were measured. The true con- 
ductivities, in the fourth column, were obtained 
by plotting (1/K.) against (1/p) and extra- 
polating to (1/p) equal zero as described in the 
Sample Calculation. 

The average percentage deviations from a 
straight line, in the (1/p) plots, are shown in the 
fifth column and the percentage deviations from 
a smooth curve of K versus temperature, in the 
sixth. These indicate the precision of the data. 

Values of the slopes in the (1/Ka,) versus (1/p) 
plots are given in the seventh column. These will 
be referred to in a subsequent paper dealing with 
accommodation coefficients. For immediate pur- 
poses, they serve as a sensitive check on the 
reliability of the K’s determined by extrapola- 
tion. Thus, the abnormally high slope in the 
149.12° runs—which is clear out of line with 
the values at adjacent temperatures—makes the 
accuracy of that K somewhat suspect. Indeed, 
this is confirmed by the relatively large deviation 
of this K from the smooth curve of K versus 
temperature, as shown in Fig. 2. 

The runs at 268.72° and at 279.45°K were 
taken with a single pressure only—namely 8.30 
cm of Hg—and were made for the purpose of 
investigating the possible effect of the magnitude 
of the temperature interval. Their K’s were ob- 
tained by using the observed d(1/K)/d(1/p) at 
273.26° for extrapolation to (1/p) equals zero. 
Neither these two runs nor the one at 149.12° 
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Fic. 2. Thermal conductivities of air. 


were given any weight in obtaining the smoothed 
curve of K versus T. 
Smoothed values of the thermal conductivities, 
at even temperatures, are given in Table V. 
The table includes temperature coefficients, 
together with the smoothed conductivities, for 
use in making interpolations. 


ACCURACY OF THE RESULTS 


The average deviation of our 14 experimental 
values (the 149.12, 268.72, and 279.45° runs 
excepted for the reasons previously given) from 
the smooth curve on which Table V is based is 


TABLE IV. Thermal conductivities of air. 








Pressures 


‘ee 7 (cm Hg, 0°C) 


d(1/K) 


% dev = 
d(1/P) 


(1/Ka) 


% dev. 
K X105 K 





87.50 

97.04 

132.33 

147.61 

149.12 

160.41 

164.13 

191.30 
217.42 
240.78 .24, 9.79 
262.99 5, 5.28, 10.50 


268.72 
273.263 , 8.94, 4.37, 1.70, 
615, 8.30 
279.45 
3 310.75 
14 337.56 
14 375.76 


45, 10.32 


3 
87, 11.39 
7 
16, 2.66, 5.68, 12.30 


—0.05 378 
+0.19 145 
—0.17 254 
+0.09 299 
—0.45* 1444 
—0.03 505 
+0.17 383 
—0.31 313 
—0.11 405 
+0.04 434 
+0.20 

+0.34* 


1.922 
2.141 0.03 
2.924 0.04 
3.268 
3.283 
3.542 
3.628 
4.171 
4.684 
5.114 
5.503 
5.607 


5.666 

5.774 

6.268 

6.705 

7.319 0.01 
Average deviation 


+0.04 
+0.17* 
—0.05 
+0.00 
—0.01 
+0.104 








* These points were given no weight in obtaining the smoothed values of K, and have been omitted from the average deviation. 
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TABLE V. Smoothed values of the thermal conductivities 
of air as a function of the temperature. 











res 4 K X105 (dK /dT) 
80 1.755 0.0224 
90 1.979 224 
100 2.203 224 
110 2.427 224 
120 2.651 224 
130 2.875 223 
140 3.097 221 
150 3.317 2185 
160 3.534 215 
170 3.747 2105 
180 3.955 2055 
190 4.158 2005 
200 4.356 1955 
210 4.549 1905 
220 4.737 1855 
230 4.920 1805 
240 5.098 176 
250 5.272 172 
260 5.442 1685 
270 5.609 1655 
280 5.773 163 
290 5.935 162 
300 6.097 162 
310 6.259 162 
320 6.421 162 
330 6.583 1615 
340 6.744 161 
350 6.905 161 
360 7.066 161 
370 7.227 161 
380 7.388 161 
273.1 5.660 1647 








+0.104 percent. This corresponds to approxi- 
mately +0.014° in the temperature interval and 
is about the precision with which we were able 
to maintain, and to measure, the temperature of 
the thermostated bath. Our precision is therefore 
limited, principally, by the evaluation of bath 
temperature. 

To evaluate the accuracy of the smoothed 
values given in Table V consideration must be 
given to the influence of possible systematic 
errors. These are as follows: 


(1) Shape Factor 


A small systematic error, which may amount 
to 0.12 percent as a maximum, may enter through 
slight uncertainties in the dimension of the cell 
or slight inaccuracy in the centering of the wire. 
This limit is arrived at by considerations of the 
figures presented in Table I, in conjunction 
with application of the formula given by Kan- 
nuluik and Martin® for the effect of the possible 
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0.014-cm inaccuracy in the centering of the wire 
(cf. infra). 


(2) Energy Input 


A maximum systematic error of +0.06 percent 
may enter through slight errors in (a) standard 
cell voltage, (b) standard resistance, and (c) po- 
tentiometer calibration—all of which contribute 
to the measured energy input. 


(3) Temperature Head Between Cell Wire and 
Cell Wall 


Systematic errors associated with evaluation 
of the radial temperature difference between the 
cell wire and the cell wall may contribute +0.20 
percent at most to errors in thermal conduc- 
tivities. This figure is regarded as a liberal 
estimate and takes into account: (a) inaccuracies 
in the potentiometer coils (0.01 percent) which 
may contribute both to thermocouple and wire 
resistance measurements; (b) drift in the stand- 
ard cell used with the potentiometer (<0.005 
percent); (c) drift in the standard 0.128477-ohm 
resistance during the month in which measure- 
ments were made; (d) error in calibration of the 
cell wire (<0.01°); (e) drift in calibration 
of the cell wire (<0.01°); (f) error in the tem- 
perature coefficient of e.m.f. of the thermocouple 
(<0.05 percent) ; (g) the very small temperature 
difference between the inner and outer surfaces 
of the cell wall (<0.01 percent); and (hk) the 
possible existence of a small temperature dif- 
ference between the outer surface of the cell wall 
and the thermocouple. The relatively large 
external surface of the metal cell and the rapid 
stirring of the bath liquid make it unlikely that 
any significant temperature difference of the 
type (h) existed. 

TaBLE VI. Limits in systematic errors introduced 
through determination of average gas temperature and 


through radiation and conduction corrections, as functions 
of the temperature (°K). 








Total 





magnitude 
of correction Limiting systematic errors 
(percent of (percent of conductivity) 
Source of error conductivity) 80° 180° 280° 380° 
Gas temperature +0.08 +0.03 +0.02 +=+0.01 
Radiation correction 0.1 to 2.5 0.02 0.05 0.07 0.18 
End conduction cor- 
rection 0.5 0.08 0.01 0.00 0.00 
Potential leads cond. 
corr. 0.15 to 0.8 0.23 0.15 0.13 0.11 
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TABLE VII. Systematic error limits (total). 








Total systematic error (percent) 





ye Random Additive 
80 +0.35 +0.79 
180 +0.29 +0.62 
280 +0.28 +0.60 
380 +0.32 +0.68 








Drift in the standard 0.128477-ohm resistance 

amounted to only 8 parts in 13,000 over a period 
of eight months and could not have much ex- 
ceeded 1 part in 13,000 during the month that 
measurements were under way. 
, The runs at 268.72°, 273.26°, and 279.45°K 
(previously mentioned) were made for the 
express purpose of finding if there might be any 
systematic trend in the results as a function of 
the magnitude of the temperature difference 
between cell wire and cell wall. To accomplish 
this, all three runs were made at a pressure of 
8.30 cm of Hg and with a bath temperature of 
266.45°K. All other conditions except the wire 
temperature were kept comparable. Tempera- 
ture heads for the three runs were, respectively, 
4.535°, 13.638°, and 25.985°. Corrected by the 
same d(1/K)/d(1/p) to zero (1/p) the K’s of 
these three runs deviated from the smooth K 
versus T curve by the respective amounts +0.34 
percent, +0.04 percent, and +0.17 percent. No 
systematic trend is evident and the 0.34 percent 
deviation corresponds to an error of only 0.015° 
in the temperature interval, which is within the 
recognized range of experimental error. We 
conclude that no systematic error has been 
introduced by our adoption of an approximate 13° 
temperature head for the majority of the runs. 


(4) Average Gas Temperature 


While temperature differences are measured 
much more accurately, we cannot claim much 
better than +0.05°, in an absolute sense, for our 
temperature scale. With account taken of some 
further error in calibration of the cell wire, we 
regard +0.06° as a liberal estimate of the uncer- 
tainty in the average temperature of the gas in 
the cell. This is equivalent to a systematic error 
of +0.08 percent in the thermal conductivity of 
air at 80°K and to +0.013 percent at 380°K. 
Equivalent errors at intermediate temperatures 
are shown in Table VI. 


(5) Systematic Errors in Corrections for Radia- 
tion, End Conduction, and Conduction of 
Heat by the Potential Leads 


The corrections referred to here have been 
discussed in a previous section, and magnitudes 
of the corrections indicated. Like the error in 
average gas temperature, systematic errors in 
making these corrections are dependent on the 
temperature at which runs are made. Table VI 
summarizes the magnitude of these corrections 
at several temperatures and of the errors in the 
thermal conductivity which may enter through 
their presence. In arriving at these estimates 
consideration has been given to the error limits 
in both the correction equations and in the data 
they employ. 


(6) Temperature-Jump Correction 


The short extrapolation in the (1/K4) versus 
(1/p) plot, which corrects for the temperature 
discontinuities between gas and wire and between 


TABLE VIII. Comparison of our own results on the thermal 
conductivity of air with those of other observers. 











This 
research Published values 

T,°K KX105 KX105 Method* Worker Year 

81.6 1.791 1.80 HW(C) Eucken*® 1913 
194.6 4.250 4.26 HW(C) Eucken*® 1913 
273.1 5.660 5.62 HW(PL) Schleiermacher> 1888 
273.1 5.660 5.72 HW(PL) Schwarze* 1901 
273.1 5.660 5.69 HW(PL) Schwarze® 1902 
273.1 5.660 (5.66) HW(C) Eucken* 1913 
273.1 5.660 5.68 HW(PL) Weber? 1917 
273.1 5.660 5.67 Putzkie 1918 
273.1 5.660 5.92 HW(PL) Schneider! 1926 
273.1 5.660 5.74 HW Webere 1927 
273.1 5.660 5.85 HW(C) Gregory and Archer» 1933 


273.1 5.660 5.81 HW(C) Milverton 1934 
273.1 5.660 5.84 HW(C) Dickensi , 1934 
273.1 5.660 5.76 HW(TW) Kannuluik and Martin® 1934 
273.1 5.660 5.72 PP Hercus and Sutherland! 1934 
273.1 5.660 5.78 HW(C) Northdurft™ 1937 
313.1 6.308 6.84 HW(PL) Schneider! c 1926 
323.1 6.470 6.68 HW(AH) Sherratt and Griffiths® 1939 
353.1 6.953 7.18 HW(C) Milvertoni 1934 
373.1 7.275 7.18 HW(C) Eucken* , 1913 
373.1 7.275 7.59 HW(AH) Sherratt and Griffiths® 1939 








* Key to symbols: 
HW(C), compensating type of hot wire cell. 
- HW(PL), potential lead type of hot wire cell. 
HW(TW), thick wire type of not wire cell. 
HW(AB), hot wire cell with auxiliary heaters, 
PP, parallel plate method. 
® See reference 5. 
b Schleiermacher, Wied. Ann., 34, 623 (1888). 
¢ Schwarze, Phys. Zeits., 3, 264 (1901, 1902); Diss. Halle, 1902; Ann. 
d. Physik 11, 303 (1903). 
4 See reference 10. , 
¢ Putzki, Diss. Halle, 1918; cf. Trautz and Ziindel, Zeits. f. tech. 
Physik 12, 273 (1931), Table I. 
f FE, Schneider, Ann. d. Physik 79, 177 (1926). 
«S. Weber, Ann. d. Physik 54, 437 (1917); 82, 479 (1927). 
b H. S. Gregory and C. T.*Archer, Phil. Mag. 15, 301 (1933). 
i See reference 17. ; 
i See reference 8. 
k See reference 16. 
1E. C. Hercus and Sutherland, Proc. Roy. Soc. A145, 599 (1934). 
m W. Northdurft, Ann. d. Physik 28, 137 (1937). 
® See reference 3. 
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gas and cell wall, is made with great accuracy. 
In the 273.26° series of runs the linearity of the 
reciprocal plot was well established in runs at 
six pressures that varied between 0.615 and 
19.55 cm of Hg (cf. Table III), and the precision 
was high (+0.02 percent) on individual runs. 
The average difference between the apparent 
conductivity measured at the highest pressure 
and the true conductivity obtained by extra- 
polation was 0.18 percent, and the maximum 
difference was 0.37 percent (except for the 149.12° 
run for which the slope was anomalously high 
and the difference 0.49 percent). We estimate 
that the error introduced by the extrapolations 
—where runs were made at two or more different 
pressures—cannot exceed about 0.03 percent. 


(7) Convection 


The linearity in the reciprocal plots referred 
to in the previous paragraph indicates freedom 
from effects of convection. This was further con- 
firmed in a separate series of runs made with 
oxygen in the cell, at 322.06°K. Seven runs 
were made in this series. Five runs in the pressure 
range 0.75 to 11.29 cm of Hg gave an average 
deviation from a straight line of only 0.01 percent, 
while runs at 28.5 and 75 cm, respectively, 
deviated from the linear plot by amounts of 
0.08 percent and 0.14 percent, respectively. The 
deviations at these higher pressures were the 
result of convection. Convection was suspected 
in these two runs, even before the results were 
calculated and plotted because a small but con- 
tinuous fluctuation of the cell resistance was 
observed that disappeared -upon reducing the 
pressure. In the measurements on air at 273.26°K 
no convection was observed at 19.55-cm pressure. 

We are confident that there is no significant 
influence of convection in the conductivities ob- 
tained for air. 


(8) Over-all Systematic Error 


The combined effect of the several sources of 
systematic error referred to in the paragraphs 
above has been calculated, as a function of tem- 
perature, on the basis of two limiting assump- 
tions: (a) that the several errors are random in 
sign and (6) that they are all additive. Results 
are shown in Table VII. 


TAYLOR AND H. L. 


JOHNSTON 


The true situation must be somewhere between 
these extremes (for the maximum error) and can 
probably be regarded as 0.5+0.1 percent. We 
thus believe that our smoothed values of the 
conductivities in Table V are correct to within 
about one-half of one percent in an absolute 
sense. Probable errors are, of course, even less 
than this. 


COMPARISON WITH THE RESULTS OF 
OTHER INVESTIGATORS 


Comparison of our own results (smoothed 
values of Table V) with those of other workers 
is shown in Table VIII. 

As is apparent from the table, very few 
measurements have been made, by other workers, 
at temperatures other than 0°C. Below the ice 
point the only published data are the two values 
reported by Eucken, which are relative to air at 
0°C. These agree well with our own measure- 
ments. 

Our own value at 273.1, which we believe to 
be reliable to within 0.5 percent, lies 2.5 percent 
below the average of values obtained with com- 
pensating cells; about 1.5 percent below Kan- 
nuluik and Martin’s result with the thick wire 
cell; 1 percent below Hercus and Sutherland’s 
determination by the parallel plate method ; and 
from 0.2 percent to 1 percent below most other 
values by the potential lead type cell. 

The errors in our own measurements have 
been carefully considered (cf. infra) and it is 
difficult to see how our results can be low beyond 
the +0.5 percent claimed for accuracy unless we 
have overcorrected for end conduction, radiation 
or conduction by potential leads. However, the 
total correction for end conduction, at 273.1°, is 
only 0.01 percent. The radiation correction is 
1.34 percent but can be applied accurately. It 


TaBLE IX. Temperature coefficient of the thermal con- 
ductivity of air at 273.1°K. 





(dK /dT) X10° 





Investigator 





This research 2.91 
Schneider 3.9 
Gregory and Archer 2.9 
Dickens 2.9 
Milverton 2.8 
Schwarze 2.5: 
Weber 3.6 
Schleiermacher 2.8 
Eucken 27 
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was determined experimentally by us and is in 
agreement with the careful work of Milverton.!’ 
Furthermore, it corresponds to a reflectivity of 
96 percent, which is in reasonable agreement with 
optical measurements for the wave-lengths of 
maximum intensity at 273.1°K (~10 microns). 
While the potential lead correction is less certain 
its full magnitude at 273.1° is only 0.51 percent. 

From the standpoint of absolute accuracy the 
advantage of our method is that the energy and 
temperature measurements are made directly on 
the effective length of the cell wire, with a cor- 
rection of only 0.5 percent necessary for the cell 
leads. The use of compensating cells assumes 
cancellation of end effects in the two cells 
amounting to 5 or 10 percent of the heat pro- 
duced in the effective length of the wire. In the 
thick wire cell as much as 50 percent of the heat 
flows from the ends of the cell, and is evaluated 


by means of an auxiliary determination of the 
thermal conductivity of the wire. It would seem 
that the possibility of serious systematic error is 
least with the potential lead method. 

The temperature coefficient of the thermal 
conductivity of air at 273.1°K has been much 
in dispute.!? The value obtained in this research 
is 0.00291. A tabulation of values obtained, by 
various workers, is shown in Table IX. 


ACCOMMODATION COEFFICIENTS 


Accommodation coefficients may be computed 
from the slope of the (1/K,) versus (1/p) plots. 
A theoretical treatment of this will be reserved 
for a subsequent paper from this laboratory. The 
paper will include tabulations of accommodation 
coefficients on bright platinum for air, from the 
data of this research, and for nine other gases, 
from data obtained by Johnston and Grilly. 
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Thermal conductivities of O2, Ne, CO, NO, He, He, N20, COez, and CH, have been measured 
between 80°K and 380°K, with the potential lead type of hot wire cell developed by Taylor 
and Johnston. Results—which have a precision generally better than 0.1 of one percent and 
are believed to be accurate to +0.5 of one percent—are tabulated for the eight gases. Com- 
parisons with the results of other investigators are also shown in tabular form. 


E have measured thermal conductivities 

of oxygen, nitric oxide, hydrogen, carbon 
monoxide, carbon dioxide, nitrous oxide, meth- 
ane, and helium over the temperature range 
80°K (or temperatures at which the vapor 
pressures amounted to about 10 cm of Hg) to 
380°K. Measurements were made at intervals of 
about 15°. 


APPARATUS AND METHOD 


Apparatus and method have been described 
in an earlier paper by Taylor and Johnston.! The 
following slight modifications were made for this 
work: 


'W. J. Taylor and H. L. Johnston, J. Chem. Phys. 14, 
219 (1946). 


(1) Thermocouple C, used by Taylor and 
Johnston, was replaced by two new copper- 
constantan thermocouples, designated couples 
V and W. These were calibrated at 90°K (b.p. 
Oz), 195°K (s.p. COs), 234°K (f.p. Hg), 273°K 
(m.p. H.O), 305°K (t.p. NaeSO,), and 373°K 
(b.p. H2O) by the method suggested in the 
Symposium on Temperature.? The two thermo- 
couples were used simultaneously to measure the 
bath temperature. They always agreed with 
each other to within 0.1 microvolt, which is 
within 0.005°, throughout the course of the 
measurements. 

* Temperature, Its Measurement and Control in Science 


and Industry, American Institute of Physics (Reinhold 
Publishing Corporation, New York, 1941). 
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Fic. 1. Shift of cell wire calibration with time. 


(2) The multiple junction thermocouple used 
in thermostating the cryostat bath was replaced 
by a resistance thermometer, which gave steadier 
thermostating. 

(3) The float system for controlling the level 
of liquid air in the cryostat was replaced by an 
air thermometer, which actuates a relay by 
making and breaking contacts with the mercury 
in a U tube with which it is joined. 

(4) Two sets of high capacity batteries were 
substituted for the one set used by Taylor and 
Johnston to supply energy to the cell wire. One 
set was used for the higher heating current (0.5 
ampere) and was kept stabilized at this discharge 
rate during the full day when measurements were 
in progress, while the other set was used and 


TABLE I. Thermal conductivity of oxygen. 











Date \ ie 4 K X105 % dev. from curve 
6/19/43 86.53 1.851 0.00 
6/19/43 94.55 2.035 0.00 

11/4/42 98.02 2.113 —0.02 
6/20/43 127.22 2.774 —0.08 
10/28/42 130.99 2.861 0.00 
6/20/43 144.56 3.166 0.00 
10/28/42 151.24 3.316 +0.06 
6/20/43 167.28 3.676 +0.22 
6/20/43 185.17 4.051 —0.22 
6/17/43 205.38 4.488 0.00 
6/17/43 229.36 4.981 +0.01 
6/17/43 250.87 5.408 —0.01 
6/17/43 272.07 5.819 0.00 
6/17/43 290.90 6.175 —0.01 
6/22/43 323.86 6.828 +0.02 
6/22/43 341.99 7.201 —0.10 
6/22/43 357.63 7.547 +0.05 
6/22/43 376.30 7.956 +0.05 


Average deviation from smooth curve +0.047% 








stabilized for the lower current (0.1 ampere) used 
in calibrating the wire. 

(5) Kerosene was substituted, as thermostat 
bath liquid, for the mineral oil employed by 
Taylor and Johnston. Because of its lower vis- 
cosity, kerosene provides more uniform ther- 
mostating. 

To eliminate “temperature jump” apparent 
conductivities were always measured at three 
different pressures. These were usually taken 
close to one, five, and ten cm of Hg, respectively. 

Calibrations of the cell wire were made simul- 
taneous with the thermal conductivity measure- 
ments, in the manner described by Taylor and 
Johnston. During the period of our work the 
resistance of the cell wire decreased slowly but 
steadily, indicative of a continuing relief of 
strain. Figure 1 shows this effect over the eight 
months period that encompassed our measure- 
ments. The radius of the circles is equivalent to 
approximately 0.02°. New calibration curves were 
constructed as often as necessary, from the cal- 
ibration data taken simultaneous with the 
conductivity measurements. 


PREPARATION OF GASES 


Pure gases were prepared by the methods out- 
lined below. Where reference is made to frac- 
tional distillation the procedure followed was the 
efficient ‘‘bubbler’” method described by John- 
ston and Giauque.* Only the middle third of the 
liquid was retained in each of the several cycles. 


3H. L. Johnston and W. F. Giauque, J. Am. Chem. Soc. 
51, 3194 (1929). 
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TABLE II. Thermal conductivity of nitric oxide. 


TABLE IV. Thermal conductivity of carbon dioxide. 











Date ie K X105 % dev. from curve Date "ie 8 K X105 % dev. from curve 
7/1/43 129.86 2.786 +0.04 11/15/43 186.38 2.073 +0.02 
7/1/43 147.82 3.165 —0.06 11/15/43 200.86 2.285 —0.02 
7/1/43 174.57 3.719 +0.01 11/15/43 216.49 2.527 —0.03 
7/1/43 191.51 4.077 +0.13 11/16/43 230.68 2.752 +0.01 
6/26/43 208.40 4.420 +0.02 11/16/43 245.57 3.002 —0.02 
6/26/43 227.76 4.793 —0.04 11/16/43 259.60 3.245 0.00 
6/26/43 245.51 5.139 —0.04 11/16/43 274.22 3.505 0.00 
3/8/44 266.20 5.544 +0.04 11/16/43 287.75 3.754 0.02 
3/8/44 283.03 5.860 —0.08 11/16/43 302.38 4.031 +0.05 
3/8/44 298.80 6.168 +0.02 11/16/43 311.24 4.199 0.00 
3/8/44 313.38 6.444 +0.01 11/13/43 333.87 4.649 0.00 
3/8/44 330.56 6.752 —0.21 11/13/43 347.32 4.932 —0.04 
3/8/44 346.25 7.056 —0.07 11/13/43 363.36 5.273 0.00 
3/8/44 361.97 7.361 +0.08 11/13/43 379.49 5.614 +0.09 
3/8/44 376.77 7.630 0.00 


Average deviation from smooth curve +0.057% 








TABLE III. Thermal conductivity of hydrogen. 











Date he 4 KxX10®° = &% dev. from curve 
11/5/43 84.70 13.90 —0.57 
11/5/43 93.99 15.41 +0.33 
12/23/43 125.60 19.91 —0.35 
12/23/43 142.23 22.46 +0.22 
12/23/43 158.70 24.89 +0.36 
11/2/43 175.95 27.19 —0.22 
11/2/43 192.28 29.38 —0.67 

3/31/44 200.40 30.68 —0.03 
11/2/43 207.66 31.62 —0.15 
3/31/44 218.18 33.09 —0.02 
3/31/44 236.56 35.56 +0.02 
3/31/44 254.38 37.68 +0.20 
3/31/44 274.02 39.75 +0.01 
11/7/43 277.69 40.05 —0.17 
11/7/43 293.28 41.64 —0.03 
11/9/43 324.28 44.53 +0.11 
11/9/43 341.27 45.97 —0.15 
11/9/43 357.31 47.51 —0.04 
11/9/43 374.25 49.10 +0.04 


Average deviation from smooth curve +0.194% 








Oxygen: Pure oxygen was obtained by passing 
commercial oxygen (99.5 percent) over potas- 
sium hydroxide, hot copper oxide, potassium 
hydroxide, and phosphorous pentoxide in turn. 
This was followed by fractional distillation. 

Nitric Oxide: Nitric oxide was prepared and 
purified by the method of Johnston and Giauque.* 

Hydrogen: Hydrogen was prepared by the 
electrolysis of potassium hydroxide solution. The 
product was passed over platinized asbestos 
followed by phosphorus pentoxide and was 
passed through a liquid air trap as the final step. 
A spectroscopic analysis by P. G. Wilkinson‘ 





*P. G. Wilkinson, Ph.D. Dissertation, The Ohio State 
University, 1941. 


Average deviation from smooth curve +0.020% 








TABLE V. Thermal conductivity of nitrous oxide. 











Date i K X105 % dev. from curve 
12/18/43 190.25 2.175 0.00 
12/18/43 206.88 2.438 —0.04 
12/18/43 225.80 2.760 +0.04 
12/18/43 232.06 2.870 0.00 
12/18/43 249.05 3.172 — 0.04 
12/18/43 262.40 3.421 +0.10 
12/19/43 275.60 3.669 +0.03 
12/19/43 290.59 3.958 0.00 
12/19/43 303.43 4.217 0.00 

2/21/44 317.55 4.511 +0.03 
2/21/44 332.42 4.829 +0.10 
2/21/44 348.08 5.150 —0.20 
2/21/44 364.11 5.505 —0.12 
2/21/44 378.45 5.836 +0.05 


Average deviation from smooth curve 








indicated that the mole percent of oxygen im- 
purity was less than 0.01 percent. 

Carbon Dioxide: Carbon dioxide was prepared 
by thermally decomposing sodium bicarbonate 
as described by Wilkinson,‘ who showed the 
oxygen content of gas so prepared to be less than 
0.001 percent. 

Carbon Monoxide: Carbon monoxide was pre- 
pared by dropping formic acid into concentrated 
sulfuric acid. The product was passed through a 
50 percent potassium hydroxide solution, dried 
over phosphorus pentoxide, and fractionally 
distilled. Wilkinson‘ obtained gas with less than 
0.001-mole percent impurity, by this method. 

Nitrous Oxide: Cylinder nitrous oxide (‘‘Puri- 
tan” hospital quality) was passed over phos- 
phorous pentoxide, and fractionally distilled. 
According to Quinn and Wernimont® and Blue 


5 E. L. Quinn and G. Wernimont, J. Am. Chem. Soc. 51, 
2004 (1929); 52, 2725 (1930). 
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TABLE VI. Thermal conductivity of carbon monoxide. 











Date ry 4 K X105 % dev. from curve 
12/20/43 87.41 1.802 —0.02 
12/20/43 97.87 2.039 +0.02 
12/21/43 130.06 2.717 —0.50 
12/21/43 145.31 3.054 —0.33 
12/21/43 155.07 3.261 +0.07 
11/19/43 176.06 3.682 —0.16 
11/19/43 192.28 4.013 0.00 
11/19/43 207.36 4.308 —0.03 
11/19/43 222.62 4.604 +0.03 
11/19/43 237.21 4.876 +0.01 
11/19/43 250.12 5.162 —0.03 
11/19/43 267.08 5.412 —0.03 
11/20/43 282.15 5.683 0.00 
11/20/43 296.90 5.950 +0.12 
11/20/43 312.00 6.210 0.00 
11/22/43 328.89 6.502 —0.09 
11/22/43 344.66 6.788 0.00 

2/23/44 357.58 7.021 0.00 
2/23/44 376.96 7.372 0.00 


Average deviation from smooth curve +0.076% 








TABLE VII. Thermal conductivity of methane. 











Date a * K X10 % dev. from curve 
1/13/44 96.76 2.451 0.00 
1/10/44 12055 3.195 —0.33 
1/10/44 140.49 3.608 0.00 
1/10/44 155.29 4.005 +0.05 
1/10/44 169.93 4.406 +0.25 
1/12/44 185.41 4.805 —0.16 
1/12/44 201.24 5.245 —0.11 
1/12/44 216.69 5.684 0.00 
1/12/44 231.71 6.117 +0.08 
1/12/44 246.86 6.554 +0.04 
1/12/44 262.85 7.025 0.00 
1/12/44 278.72 7.510 0.00 
1/12/44 294.49 8.026 +0.25 
1/14/44 324.00 8.986 —0.25 
1/14/44 338.60 9.539 +0.02 
1/14/44 353.96 10.146 +0.10 
1/14/44 368.33 10.729 0.00 
1/14/44 383.57 11.376 +0.01 


Average deviation from smooth curve +0.092% 








and Giauque® the impurities in this grade of 
commercial nitrous oxide are nitrogen, oxygen, 
and water. These should all be easily removed 
by fractional distillation. 

Methane: Methane was obtained from natural 
gas’ by passing the latter through a trap at 


6 R. W. Blue and W. F. Giauque, J. Am. Chem. Soc. 57, 
992 (1935). 

7 The following analysis was furnished us by Mr. A. M. 
Hutchison of the Ohio Fuel Gas Company, as typical of 
the natural gas in the Columbus mains: carbon dioxide, 
0.5 percent; oxygen, 0.0 percent; hydrogen, 1.3 percent; 
carbon monoxide, 0.6 percent; nitrogen, 5.0 percent; 
illuminants, 0.0 percent; methane, 80,6 percent; higher 
hydrocarbons, 12 percent. 
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193°K, followed by two consecutive tubes of 
charcoal’ at 193°K. After this the product was 
fractionally distilled. Storch and Golden® used 
this method on natural gas of similar composition. 
In their experience they found the following 
concentration of impurities after the passage 


TABLE VIII. Thermal conductivity of helium. 











Date re, K X105 % dev. from curve 
2/1/44 82.79 15.72 0.00 
2/1/44 92.71 16.96 +0.15 
2/3/44 122.67 19.91 —0.02 
2/3/44 137.90 21.41 — 0.42 
2/3/44 154.30 23.13 —0.16 
2/3/44 169.89 24.72 +0.05 
2/3/44 184.32 26.09 0.00 
2/4/44 200.24 27.50 —0.32 
2/4/44 216.89 29.12 +0.04 
2/4/44 231.74 30.42 — 0.04 
2/4/44 248.42 31.87 +0.02 
2/4/44 264.36 33.19 —0.03 
2/4/44 280.53 34.50 0.00 
2/4/44 296.18 35.74 +0.08 
2/18/44 328.79 38.05 —0.04 
2/18/44 344.43 39.12 +0.06 
1/31/44 359.59 39.99 —0.07 
1/31/44 375.68 40.94 0.00 


Average deviation from smooth curve +0.083% 








TABLE IX. Smoothed thermal conductivities (K X 105°). 








T,°K Or co He H2 CHa NO CO: NO 





80 1.701 1.652 15.34 13.31 

90 1.930 1.870 16.43 14.78 2.272 
100 2.159 2.086 17.51 16.25 2.536 
110 2.387 2.302 18.58 17.71 2.800 
120 2.614 2.516 19.64 19.17 3.065 2.580 

130 2.840 2.730 20.68 20.63 3.331 2.792 
140 3.064 2.942 21.71 22.09 3.595 3.003 

150 3.287 3.153 22.73 23.54 3.860 3.214 
160 3.508 3.361 23.73 24.98 4.128 3.423 

170 °3.728 3.565 24.72 26.42 4.396 3.631 

180 3.946 3.768 25.68 27.84 4.667 3.838 1.984 2.021 
190 4.162 3.968 26.63 29.25 4.940 4.042 2.126 2.173 
200 4.375 4.166 27.56 30.64 5.216 4.245 2.272 2.330 
210 4.584 4.361 28.49 32.00 5.496 4.445 2.424 2.493 
220 4.790 4.553 29.39 33.33 5.778 4.643 2.580 2.661 
230 4.993 4.741 30.28 34.64 6.063 4.840 2.741 2.823 
240 5.194 4.926 31.15 35.89 6.351 5.035 2.907 3.011 
250 5.392 5.107 32.00 37.09 6.643 5.229 3.077 3.192 
260 5.586 5.286 32.83 38.23 6.940 5.423 3.251 3.377 
270 5.780 5.465 33.65 39.32 7.242 5.615 3.429 3.564 
280 5.970 5.644 34.45 40.35 7.549 5.807 3.611 3.755 
290 6.159 5.821 35.23 41.33 7.862 5.999 3.796 3.949 
300 6.350 5.998 36.00 42.27 8.186 6.189 3.984 4.149 
310 6.547 6.175 36.74 43.19 8.518 6.379 4.175 4.352 
320 6.748 6.351 37.46 44.11 8.862 6.568 4.371 4.562 
330 6.954 6.528 38.15 45.02 9.219 6.756 4.571 4.775 
340 7.164 6.705 38.81 45.94 9.590 6.944 4.777 4.991 
350 7.378 6.884 39.44 46.85 9.978 7.131 4.988 5.206 
360 7.594 7.065 40.05 47.77 10.372 7.318 5.202 5.425 
370 7.812 7.246 40.62 48.69 10.797 7.504 5.416 5.644 
380 8.033 7.427 41.17 49.60 11.220 7.690 5.630 5.861 


273.1 5.839 5.521 33.90 39.65 7.336 5.674 3.485 3.623 
293.1 6.218 5.876 35.48 41.64 7.962 6.057 3.854 4.011 
298.1 6.314 5.964 35.86 42.10 8.123 6.153- 3.948 4.110 








8 The charcoal was baked under vacuum at 450°C for 
forty hours and was maintained at 193°K while the gas 


was passed through. 
9H. H. Storch and P. L. Golden, J. Am. Chem. Soc. 54, 


4662 (1932). 
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TaBLE X. Comparisons of our own results on thermal conductivities at 273.1°K with those reported by other observers. 











% dev. % dev. 
from from 
Gas KxX10° J andG Method* Author Year Gas KX105 J and G Method* Author Year 
Oo 5.839 Pi This research 1944 He 39.65 at Fs This research 1944 
5.89 +0.9 3 Northdurft® 1937 41.82 +5.4 e Archeri 1938 
5.90 +1.1 > Dickins» 1934 42.45 +7.4 c Northdurft« 1937 
5.83 —0.2 TW Kannuluik and Martine 1934 42.69 +7.5 CG Gregory* 1935 
5.89 +0.9 ‘ Gregory and Marshall 1928 41.70 +5.1 ¢ Dickins> 1934 
5.768 —1.2 a Webere 1917 41.3 +4.0 TW Kannuluik and Martine 1934 
(5.839) 0.0 = «C Euckenf 1911 40.62 +2.5 c Hercus and Laby! 1922 
Av. +0.0 42.38 +6.8 Pls Weber™ 1927 
40.43 +1.9 C¥* Gregory and Archer® 1926 
NO 5.674 PL This research 1944 41.75 +53 P Ds Schneider® 1926 
5.680 +0.1 Cc Euckent 1913 41.65 +50 PL Webere 1917 
. : 40.8 +2.7 C Euckene 1913 
40.6 42.5 C Eucken! 1911 
co §.521 P i This research 1944 Av. +4.4 
5.58 +1.0 § Dickins> 1934 
5.37 —3.0 TW Kannuluik and Martine 1934 COz 3.485 FP ft This research 1944 
5.633 +2.0 Cc Gregory and Archer» 1928 3.47 —0.4 i Archer° 1935 
5.557 +0.6 & Euckens 1913 3.51 +0.7 c Dickins» 1934 
Av. +0.1 3.43 —1.5 TW Kannuluik and Martine 1934 
- 3.604 +3.4 C*¥* Gregory and Marshall 1927 
CH. 7.336 PL This research 1944 — a= -- eo +4 
a “aR a , 3.393 —2.7 PL Webere 1917 
7.21 1.7 Cc Dickins and Manni 1931 . “tens 
oa f 3.442 —1.2 & Euckené 1913 
7.200 —-1.8 Pf, Webere 1917 3.421 -~18 C Tesclaon? 1911 
7.319 —0.2 C Euckent 1913 tees 06 —— 
Av. 1 ee 
N20 3.623 PL This research 1944 
He 33.90 PL This research 1944 3.64 +0.4 c. Dickins> 1934 
35.10 +36 C Dickins> 1934 3.61 —0.3 j Kannuluik and Martine 1934 
34.30 +1.2 TW Kannuluik and Martine 1934 3.740 +3.2 C¥* Gregory and Archer» 1928 
34.38 +1.4 ri, Webere 1917 3.530 —2.6 PL Weber® 1917 
34.21 +0.9 c Euckenf 1911 3.601 —0.6 Gc Eucken¢ 1913 
Av. +1.5 Av. 0.0 














® W. Northdurft, Ann. d. Physik 28, 137 (1937). 

+ B. G. Dickins, Proc. Roy. Soc. (London) A143, 517 (1934). 

eW. G. Kannuluik and L. H. Martin, Proc. Roy. Soc. (London) 
A144, 496 (1934). 

4H. Gregory (Spencer-Gregory) and S. Marshall, Proc. Roy. Soc. 
(London) A114, 354 (1927); A118, 594 (1928). 

¢S. Weber, Ann. d. Physik 54, 437 (1917). 

f A. Eucken, Physik. Zeits. 12,°1101 (1911). 

« A. Eucken, Physik. Zeits. 14, 324 (1913). 

» H. Gregory (Spencer-Gregory) and C. T. Archer, Proc. Roy. Soc. 


(London) Allo. 91 (1926). 
iB. G. Dickins and W. B. Mann, Proc. Roy. Soc. (London) A134, 


77 (1931). 
ic. Tf. Archer, Proc. Roy. Soc. (London) A165, 474 (1938). 


through charcoal: oxygen, 0.2 percent; nitrogen, 
2.0 percent; higher hydrocarbons, less than 0.05 
percent. We believe that the conditions that per- 
tained during our purification of methane were 
such as to yield at least as high purity at corre- 
sponding stages of the process. Our final product 
showed a triple point temperature of 90.63 
+0.04°K which compares with the value 90.67 
+0.03°K reported by Clusius and Weigant,!° for 
pure methane. It yielded a vapor pressure of 
82.4 mm of Hg at 90.06°K. There is no data in 
the literature which can serve as a comparison 
for the vapor pressure. However the triple point 
comparison sets an upper limit of 0.15 percent 
for the mole fraction of impurity in our methane. 
The probable purity is considerably higher than 
this. 

Helium: The helium was bought from the 





© K. Clusius and K. Weigand, Zeits. f. physik. Chemie 
B46, 1 (1940). 


an S. Gregory (Spencer-Gregory), Proc. Roy. Soc. (London) A149 
35 (1935). 
1 E. O. Hercus and T. H. Laby, Phil. Mag. 3, 1061 (1927). 
mS, Weber, Ann. d. Physik 82, 479 (1927). 
» E. Schneider, Ann. d. Physik 80, 215 (1926). 
°C. T. Archer, Phil. Mag. 19, 901 (1935). 
* Key to symbols: P L Potential lead type of hot wire cell, 
4 Compensating type of hot wire cell, 
T W Thick wire cell. 
** Temperature jump correction omitted or uncertain. 
*** The hydrogen prepared by Weber was purified with charcoal at 
80°K, which might have increased the parahydrogen content and thus 
the thermal conductivity. 


Linde Air Products Company as “spectro- 
scopically pure.’’ A spectroscopic analysis made 
on a portion of this gas, in connection with a 
previous investigation,!! gave evidence of barely 
perceptible traces of argon and of neon. A repeat 
analysis, just before our thermal conductivity 
measurements were begun, yielded a similar 
result. 


EXPERIMENTAL RESULTS 


Experimental values of the thermal conduc- 
tivities obtained for the several gases are pre- 
sented in Tables I to VIII. The first column in 
each of these tables gives the dates on which the 
data were taken. The second column gives the 
average gas temperatures (cf. reference 1). The 
third column gives the thermal conductivities. 
These are obtained by extrapolating the re- 


H. L. Johnston and E. R. Grilly, J. Phys. Chem. 46, 
948 (1942). 
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TABLE XI. Comparisons of our own results on thermal 
conductivities with those reported by other investigators, 
at temperatures other than 273.1°K. 











K X105 % dev. 
This Other from 
Gas T,°K res. res. J andG Author Year 
COz 285.1 3.703 3.69 —0.3 Archer® 1935 
309.7 4.169 3.99 —4.3 
347.7 4.940 4.65 —5.8 
373.6 5.493 5.12 -—7.0 
194.6 2.193 2.24 2.5 Eucken> 1911 
373.0 5.480 5.08 —8.0 
339.3 4.763 4.95 +4.0 Sherrat and 1939 
341.7 4.812 4.76 +1.1 Griffiths 
342.1 4.821 4.80 —0.4 
380.0 5.630 (5.62) —0.2 
He 94.6 15.46 15.3 —1.0 Spencer-Gregory 1938 
96.6 15.75 15.5 —1.5 and Dock4 
96.3* 15.71 15.6 —0.7 
99.1 16.11 15.9 —1.3 
99.3 16.15 16.0 —0.7 
202.5 30.98 31.1 +0.3 
209.4 31.92 31.7 —0.7 
222.7 33.68 33.9 +0.7 
224.5 33.92 34.5 +1.8 
226.0 34.12 35.0 +2.6 
280.4 40.39 42.2 +4.5 
81.6 13.54 13.56 +0.2 Eucken® 1913 
81.5 13.53 13.63 +0.7 Eucken> 1911 
83.4 13.81 14.16 +2.5 
90.4 14.84 15.16 +2.2 
194.6 29.89 31.40 +5.0 
373.0 48.96 51.16 +4.4 
Oz 81.6 1.738 1.763 +1.5 Eucken> 1911 
90.4 1.939 1.991 +2.6 
194.6 4.260 4.397 +2.7 
373.0 7.878 7.608 —3.4 
co 82.1 1.698 1.690 —0.5 Euckene 1913 
91.6 1.905 1.890 —0.8 
198.4 4.134 4.118 —0.4 
He 81.4 15.49 15.20 —2.0 Euckene 1913 
80.3 15.39 14.92 —3.0 Eucken> 1911 
194.6 27.06 26.96 —0.4 
373.0 40.78 40.82 +0.1 
NO 201.7 4.279 4.262 —0.4 Euckene® 1913 
NO 201.3 2.351 2.776** +-18.1 Euckene 1913 
CH, 91.5 2.312 2.303 —0.4 Euckene 1913 
197.5 5.147 5.061 —1.7 











® See Table X, refernce o. 
>See Table X, reference f. 
eG. G. Sherrat and E. Griffiths, Phil. Mag. 27, 68 (1939). 
4H. Spencer-Gregory and E. H. Dock, Phil. Mag. 25, 129 (1938). 
¢See Table X, reference g. 
* Probably intended for 97.3°. 
** This abnormally large value must be due to the presence of air. 


ciprocals of the apparent conductivities (cf. 
reference 1) against the reciprocals of pressure— 
for pressures of approximately one, five, and ten 
cm as referred to above—to (1/p) equals zero. 
The fourth column lists the percentage devia- 
tions from smooth curves of the K’s versus tem- 
perature. Thermal conductivities are expressed 


H. L. JOHNSTON AND E. R. GRILLY 


in the usual units, namely calories per second 
per centimeter per degree. 

Smoothed values of the thermal conductivities, 
at even temperatures, are given in Table IX. 

A good estimate of the precision of the 
measured conductivities is obtained by reference 
to the average percentage deviations from the 
respective smooth curves, in the plots of K versus 
temperature. Except for hydrogen, for which the 
average deviation amounts to almost 0.2 percent, 
these range from 0.02 percent (carbon dioxide) 
to 0.09 percent (methane). The absolute accuracy 
of the smoothed values should be about the 
same as obtained for air, by Taylor and Johnston! 
—namely +0.5 percent. The careful error 
analysis made in that paper should be applicable 
here, with little modification. 


COMPARISONS WITH THE RESULTS REPORTED 
BY OTHER INVESTIGATORS 


Few measurements of thermal conductivity 
have been reported at temperatures other than 
the ice point. Comparisons with results reported 
by other investigators at 273.1°K are listed in 
Table X, for the percentage deviations from our 
own values, and with a reference to the method 
employed. 

Comparisons at temperatures other than 273.1° 
are given in Table XI. Our own values in this 
table are from our smooth curves. 

All of the values in Tables X and XI are ob- 
tained by methods which yield absolute values 
for the conductivity, except those of Eucken, 
who expressed his data relative to air at 273.1°K. 
For our purposes, his values have been shifted to 
correspond to Ko73.1= 5.839 X10- for oxygen. 


ACCOMMODATION COEFFICIENTS 


Accommodation coefficients may be computed 
from the reciprocal plots of apparent conductivity 
versus pressure. In a subsequent paper we will 
give accommodation coefficients on bright plati- 
num for the several gases included in this 
research. 
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Formulas are derived which make the tables and equations of the first paper of this series 
applicable to this more general class of molecules. Additional approximations are involved and 


these are examined carefully. 





INTRODUCTION 


N the first paper of this series,! hereafter 

designated I, a treatment was given for all 
molecules which could be regarded as rigid 
frames with attached symmetrical tops. The 
present paper will consider molecules composed 
of a rigid frame with attached, unsymmetrical 
tops. A third paper will discuss compound 
rotations—rotating groups attached in turn to 
other rotating groups. 

With the more complex situation it does not 
seem feasible to attain the same accuracy and 
rigor as in I. However, it seems probable that the 
methods proposed will give thermodynamic 
results with errors not appreciably greater than 
the best current experiments. Also the energy 
separations of the lowest and first excited levels 
should be given without excessive error. 

Unsymmetrical, attached rotating tops fall 
into three classes. First, balanced groups, where 
the center of gravity of the top is on its axis of 
rotation. The balance may be accidental or due 
to a twofold symmetry axis (if a three- or higher- 
fold axis is present the top is symmetrical). A 
phenyl group is an example. Second, groups with 
a small “‘off balance’’ factor such as OH, NHz, 
CH.D, CCL**,Cl37. For each of these classes it 
seems unlikely that the over-all moments of 
inertia will change greatly with internal rotation. 
Since the properties of interest are not sensitive 
too small changes in moments of inertia, it is 
permissible to ignore them and treat the molecule 
in terms of the various quantities appropriate to 
the equilibrium configuration. 





1K. S. Pitzer and W. D. Gwinn, J. Chem. Phys. 10, 
428 (1942). ; 


The third class comprises groups with large 
off-balance factors such as CHCl, CHBr2, CHO, 
etc. Here the changes in moment of inertia are 
too large to dismiss. However it must be realized 
that vibrational motions are also affected. Con- 
sequently, any attempt to correct for these 
effects on rotation alone will have little practical 
value. Indeed it can be shown that the changes 
in the mass or inertia factor for the various 
normal modes of vibration just cancel the 
changes in moment of inertia provided all 
motions are treated classically. Since vibration 
frequencies are usually known only for the 
equilibrium configuration, it seems best again 
to take all quantities appropriate to the equi- 
librium configuration and to ignore changes with 
angle in the first approximation. However this 
leaves the question of how nearly classical are 
the vibrations that effectively compensate the 
changes in moment of inertia. Study of indi- 
vidual cases indicates that the answer will almost 
always be favorable, however each case should 
be studied from this point of view. 

In molecules of these types the potential 
energy function is likely to be much more com- 
plex than with symmetrical groups. The direct 
numerical integration of the classical partition 
function is a practical procedure and will be of a 
good approximation for cases with large moments 
of inertia. The quantum mechanical energy 
levels are formally given by Eq. (18)-I and the 
associated discussion for a reasonably general 
case. However, the labor involved in this method 
may be large. Once the energy levels are ob- 
tained the partition function sums are readily 
computed. 
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DEFINITIONS 


For the entire molecule, of mass M, including 
_all rotating groups in their equilibrium positions, 
find the principal axes of rotation through the 
center of gravity and the principal moments of 
inertia,? J,, Is, and J;. 

For each top, m, take the rotation axis as z 
and make x pass through the center of gravity 
of the top. If the top is balanced, define x 
arbitrarily and in any case make y perpendicular 
to both x and z. Although it makes no real dif- 
ference, it is simpler not to include atoms on the 
axis of rotation as part of the top. Thus the 
rotating group is just the off axis H atom for an 
OH radical. If m; is the mass of the 7th atom, then 


An=DLim(x?+y2) the top moment of inertia 
about z. ° 
Bn= di mx 2; the xz product of inertia. 
Cn = >i miyz; the yz product of inertia. 
Um = >>; mx; the off-balance factor. 


In each case the sum is over only those atoms in 
mth top. The array of quantities: 


9 
n= n*® Cin* 
Qn 4 Om?4 Qm>Y 
On! 2 Om? z Om? z 


are the direction cosines between the axes of the 
mth top (x,y,z) and the axes of the whole 
molecule (1, 2,3). (All sets of axes must be 
uniformly either right-handed or left-handed. 
To check this verify that the determinant of a’s 
is +1.) The vector from the center of gravity of 
the whole molecule to the origin of coordinates 
of the top is r, with components f»!, 7m’, 7m? on 
the principal axes. 
mam 4 UnmUm  Bm'Bms* 


ja = ’ 


‘ M I; 





where 
s..* _ Om'?A m~ tin**B.. =a rm 
+ U,.(an** vr,,*t 1 Om't! *7,,.**). 


The superscripts i—1 and i+1 refer to cyclic 


? The calculations can be made without finding principal 
axes by using the products of inertia. However, the equa- 
tions are much more complex so it is doubtful if that 
method would ever be much easier. One principal axis is 
obvious in most all molecules, either a rotational sym- 
metry axis or perpendicular to a plane of symmetry. 


S. PITZER 


shifts of axes so that if 7<=1, —1=3 and if 7=3, 
1+1=1. 

The reduced moment of inertia for internal 
rotation if there is only a single rotation (or the 
first approximation if there are several) is: 


(m*Y =)? (Bm*)? 
In=Am—Amm=An— + ‘ 
2 M I; 





(1) 


A second approximation is given by Eq. (1b)-I 
but using the new definitions of Amm’ given 
above. 

For the case of balanced tops, all of these 
additional terms are zero, so that the old for- 
mulas apply. 


SUMMARY OF RESULTS 


All essential results are now given by the 
equations and tables of paper I using the new 
definitions of J and A. 

If the potential energy can be represented by: 


7 = , 3 roth — CUS NnPm), (2) 


then the tables in paper I apply. However, if the 
symmetry number, ”», corresponds to an ap- 
proximate symmetry of the potential curve 
rather than an exact symmetry of the rotating 
group, then (RInm»,) must be added to the 
functions S and (H,°— Fr°)/T. 

Examples: 
Sym. ethane d, (CH2)—CH2D). One CHD 
is the frame; the other is the rotating group with 
a small off-balance factor and an approximate 
symmetry number of 3 for the potential energy. 
Use mm=3 for the tables but add RIn3 to 
entropy and free energy function. 
Ethylbenzene. With the CH as the rigid frame, 
the methyl group is symmetrical; the phenyl is 
balanced; the internal symmetry numbers are 
3, and 2, respectively. Presumably there is no 
large cross term in the potential energy. 
m-Chloro ethylbenzene. Now the phenyl group 
has an off-balance factor and its symmetry 
number of 2 is only an approximation for the 
potential energy. The tables of paper I may be 
used with ~=2 (for the phenyl rotation) but 
add RiIn2 to the entropy and free energy 
function. 

If the potential energy cannot be expressed by 
Eq. (2), then the tables in I cannot be applied 
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but the methods described in the last section of 
I must be used. 


PROOFS 


Two mathematical arguments remain to be 
given. First, that with the new definitions, the 
kinetic energy is given -by Eq. (7)-I for motion 
at the equilibrium configuration. 

Let there be a frame of axes of fixed direction 
with origin at the center of gravity when all 
particles are in original positions. Then (0) is a 
set of axes with the same origin but rotating 
with angular velocity wo. Then for each rotating 
group there is a set of axes (m) whose origin is 
fixed in (0) and which rotates with angular 
velocity , with respect to (0). Actually the 
direction of @,, though arbitrary is also fixed 
in (0). Each particle, of mass m;, is located in 
one of these frames by a vector oo; Or omni. The 
vector fr, is from the origin of (0) to the origin 
of (m). 

The velocity of a particle 7 in the main frame is 


Vi=@0 X01, (3) 
while that of a particle mj in the group m is 
Ving =@0 X (Lm+ mj) +m X Om;- (4) 
The kinetic energy is then 


27 =D MNi-Vit Dm Dj MmVnj Vins (5) 





where the sums cover all appropriate particles 
such that each particle is included once. 

Now in general the center of gravity may 
move with a velocity 


v=(LimviXdn >; Mm iWVmj)/M. (6) 


However, this can only arise from the group 
rotations because the main frame rotation is 
about its center of mass. Hence 


v= ae Wm x > 5 Mm j0mj)/M. (7) 


The kinetic energy of this center of gravity 
motion is 
2T” = Mv-v, (8) 


and must be subtracted from 7” to give the 
true kinetic energy of rotation T. 


2T={d; m(@oXo,)? 
+ Yom 05 Mm iLooX (m+ Om;) ]*} 
+ Yom 205 Mmj{2L@o X (Lm+ Omi) J 
‘Lom X Omj J+ Lom X ems ]*} 


1 
—={ Lin Om X Lj MnjOmj}*. (9) 
M 
Let us now expand the vectors and assume the 
principal axes are used for the system (0). We 
replace wm by the scalar ¢,, used in I, and add 
primes to the components of ao. 


3 3 N 
2T= 3 T(w;’)?+2 7 Zz w UmPm(am'| ¥,*t _— ay *** Y7»,*—) 
i=l 


i=1 m=1 


N 


3 } 
+2 .&  ® wi Pm(Qm* 7A mm Am" *5.. a Om" ¥Cm) 


i=l m=1 


N 1 N N 3 
+ si AnGmn? —— > _ GmUmeOm' Um , > Qm* Yotmr” ¥ ° (10) 


m=1 
Introduce the substitution, 


w:’ =wi-—Dimm Bom'/Ti, (11) 
then 


3 N N N 
2P=) Iw ?+ > Ann? — > os PmPm 


i=1 m=1 m=1 m’=1 


Om” YOtmr™ YU mn Um" Bm'Bm'* 
Xx » az 
CS * 





which is Eq. (7)-I when the new value of Anm: 
is used. 


M m=1 m’=1 i=l 





The second mathematical argument concerns 
the compensation of changes in moment of 
inertia with changes of vibration frequency. Let 
us first take as an example a chain of four suc- 
cessive atoms at tetrahedral bond angles. The 
internal rotation is about the central bond and 
the rotating top is grossly unsymmetrical. 
Kassel® has given the needed equations for this 
case. The appropriate classical partition function 

3L. Kassel, J. Chem. Phys. 3, 326 (1935); 4, 276 (1936). 


In the latter article the terms for the hydrogen atoms are 
removed from the equations for n-butane. 
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for vibration and rotation is: 


(=) (—) **[ S}? 
J = e-Viviktd y, 
o 


Il»: 


i= 1 








(13) 


Where a is the number of modes of vibration, 
vy; the frequency of the ith mode, S is the 
effective product of moments of inertia for 
external and internal rotation, and ¢ is the angle 
of internal rotation. 

From Kassel’s work we find 


a kh? 
ann(A) CY 
inl 2a 
where 
Y?= (9+ 247? —8y*) (33 —8y?—16y*) (15) 
and 
7 =cos (¢/2) 
(Kassel’s @ or a is half of ¢ as used here). 
[ S}t=8(mr_?/9)?Y. (16) 


Thus the factor Y which expresses the change 
in the product of vibration frequencies with 
internal rotation exactly compensates the same 
factor in the product of moments of inertia. 
Consequently, the vibration frequencies and 
moments of inertia could have been computed 
for any particular configuration and would have 
given the correct result. 

It would be quite impractical to work out a 
general case by this method; the generality of 
the result may be seen by using the method 
introduced by the writer for long chain mole- 
cules. Here one avoids normal coordinates 
completely. 

Let us assume there are ~ atoms in the rigid 
frame of which the 1st and 2nd atoms are on the 
axis of internal rotation. Then the atoms of the 
rotating top, +1 to m, are attached to the 1st 
atom or to others in the rotating top. The 
partition function is given by an equation 
analagous to Eq. (12) of reference (4). 





=f. P sf ewrtde, -+dZn, (17) 
QakT \ 32/2 2 

A= = II m;3, (18) 
h? i=l 


4K. S. Pitzer, J. Chem. Phys. 8, 711 (1940). 
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where m; is the mass of the ith atom, V is the 
potential energy, and x,—z, are coordinates 
describing the positions of the various atoms. 
The integrations are, strictly, over the volume 
of the container. 

Selecting arbitrarily atom 1, it may be placed 
anywhere in the container of volume v without 
change in potential energy. Hence the integra- 
tion gives just the factor v. Converting to polar 
coordinates for the second atom (relative to 
atom 1), integration over the angles (with no 
change in potential energy) gives 4m so that now 


Q=B f+ fev*tdy: «dan, 


R 
B=Adro [ exp [— V(n1)/RT ]ri2dri. (20) 


(19) 


The upper limit of integration in (20), R, is 
any number several times larger than the equi- 
librium value of 7:1. The whole value of the 
integral comes from the region near 7; so the 
exact value of R is not important, indeed for a 
harmonic force law it may be infinity. For other 
force laws, however, the integral might not 
converge at infinity. 

Selecting as 3 an atom bonded to atom 2 and 
at an angle less than 180° from the 1-2 bond, 
one again used polar coordinates. Now the 
potential energy depends on 7 and @ but not ¢. 


Q=Cf + ferrdrs don, (21 
C= Bax { f exp [ — V(rs, 03)/RT ] 
Xr? sin 63073403. (22) 


From now on each atom in the frame can be 
treated in terms of coordinates defined with 
reference to lower numbered atoms with the 
potential energy depending on 7, 6, and ¢. 


o=Df. ’  fevrnrde,: + +dZn, (23) 
Pp 2r ¥ R 
p-cii( f f f exp [— V(r:0:¢;)/kT ] 
ima \J0 Jo Jo 
Xr? sin 6,drd0d¢;). (24) 


The first atom in the rotating top defines the 
position of the top. Hence the potential energy 
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of.internal rotation is V(¢g p41) while r and @ have 
typical bond stretching and bending energies. 
The remaining atoms in the top give terms 
similar to the later atoms of the frame. 


Q=DEF J "exp [—V(epe)/kT Meni, (25) 


e=-{ f * exp C— V(rps1pus)/kT 


Xo? sin Op4107 p41 8 p41, (26) 


P= 11 (fff eet- V(r:0:¢:)/RT ] 


Xr? sin 6.drd0d¢;). (27) 


Thus we see from Eq. (25) that the factor for 
internal rotation is independent of the other 
coordinates and has essentially the same form as 
with symmetrical groups. This proves that the 
result obtained with the four atom chain is 
general and that for a classical treatment it is 
permissible to calculate vibration frequencies 
and moments of inertia for just a single position 
with respect to internal rotation. Of course, if 
desired this method may be employed to actually 
calculate the classical partition function.‘ 

This argument was based, for simplicity, on a 
reasonably simple valence force potential energy 
function in order to obtain the explicit form of 
Eq. (25). However, the critical point is the 
integration of the kinetic energy terms to give 
Eqs. (17) and (18) in which, so far as inertia or 
kinetic energy terms are concerned, all incre- 
ments of internal rotation angle are weighted 
equally. Thus, the basic conclusion, that the 


changes in inertial quantities cancel one another, 
is quite independent of the force field assumed 
here. However, if there are cross terms in the 
potential energy between internal rotation and 
vibrational motions, then the use of vibration 
frequencies for the equilibrium configuration 
will represent an approximation. This could be 
removed by combining Eq. (25) with others as 
necessary and integrating the multiple integral 
with the cross terms in the potential energy 
present. 

It is true that for most molecules of interest 
only some of the vibrations can be treated clas- 
sically. Very roughly, vibrations may be classified 
as (1) hydrogen vibrations, (2) bond stretching 
motions of heavy atoms, and (3) bending motions 
of heavy atoms. The frequencies in the first two 
groups are largely independent of the arrange- 
ment of more distant parts of the molecule. 
Consequently, the third group must be prin- 
cipally responsible for the changes which com- 
pensate the moments of inertia. However, the 
vibration frequencies of this third group are 
usually below 400 cm and consequently nearly 
classical motions at or above room temperature. 

Taking again the four atom chain as studied 
by Kassel,* we find that the product of stretching 
frequencies varies with internal rotation only by 
5 percent while the product of bending fre- 
quencies changes by 45 percent. While we have 
no general proof of this point, it seems probable 
that this example is typical. If so, the method 
proposed above will be a good approximation. 

The writer wishes to thank Professor W. D. 
Gwinn for checking these derivations. 
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The rate of permeation of a gas through a polymer as a function of temperature may be 
represented as P = Poe~#/#7, All available data on the permeability of gases through polymers 
show that for a given gas there is a linear relation between log Po and E (the energy of activation 
for permeation). An explanation is offered for this apparent relationship. The effect of plasti- 
cization on permeation of water vapor is studied experimentally and it is shown that the 
lowering of the heat of solution is the predominant effect. From the data the entropy of solution 
may be calculated and interpreted as showing that water molecules dissolved in the unplasti- 
cized polymer exhibit much less freedom than when they are dissolved in plasticized polymer. 
Other calculations show (in agreement with Barrer) that a large number of degrees of freedom 
are simultaneously operative in the process of the diffusion of a gas molecule through a polymer. 





INTRODUCTION 


HE studies of Barrer! together with his 

compilation and analysis of existing data? 
have contributed most significantly to our funda- 
mental understanding of the transmission of 
gases and vapors through polymeric substances. 
These investigations demonstrated that trans- 
mission results from activated diffusion and that 
the molecular diameter of the diffusing gas is of 
only secondary importance. The process of 
permeation is visualized as follows when a film 
of polymeric substance separates two quantities 
of gas or vapor at different pressures. The gas 
molecules dissolve in the polymer and move 
about by jumping into holes which form and 
disappear in their immediate neighborhood. 
These holes are the result of fluctuations in the 
thermal motion of segments of the polymer 
molecules. The concentration gradient in the 
film causes a drift of the diffusing gas toward 
the lower concentration region. 

The rate of transfer of vapor across a film of 
unit thickness per unit area per unit pressure dif- 
ference is known as the permeability constant and 
characterizes the effectiveness of the film as a bar- 
rier. It has been shown that the permeability con- 
stant, P, may be resolved into the product of the 
diffusion constant, D, and the solubility coeffi- 
cient, S. The dimensions of D are area per unit 
time ; S has dimensions of volume of gas (at stated 


1R. M. Barrer, Trans. Faraday Soc. 34, 849 (1938); 
35, 628 (1939); 35, 644 (1939); 36, *1235 (1940). 
2R. M. Barrer, Diffusion in and through Solids (Cam- 
bridge University Press, New York, 1941). 


conditions) per volume of polymer. With a 
proper choice of units we may write: P=DS. 
This relation has been verified for permanent 
gases’? and water vapor® in a number of poly- 
meric substances. 

It is important to note that it has likewise 
been verified that P, D, and S vary exponentially 
with reciprocal absolute temperature for perma- 
nent gases. Thus: 


P=P, exp (—E/RT), (1) 
D=D, exp (—Ep/RT), (2) 
S=S,exp (—AH/RT). (3) 


Since P and D are essentially rate constants, 
E and Ep may be interpreted as energy of 
activation for permeation and diffusion, respec- 
tively. The solubility coefficient is an equilibrium 
constant, hence AH may be identified as the 
heat of solution. Moreover, Py and Do correspond 
to frequency factors in the Arrhenius equation 
of chemical kinetics. Barrer’? has interpreted 
his measurements of these quantities in terms of 
the transition state theory to show that the 
activated state for diffusion involves many 
degrees of freedom indicative of the formation 
of a “hole” for the diffusion molecule to move 
into. 

In this paper it is desired to present some 
further interpretations of this problem based 
upon investigations of water vapor-polymer 
systems, some of which have been previously 


3P. M. Doty, W. H. Aiken, and H. Mark, Ind. Eng. 
Chem. (in press). 
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published.’ Specifically, we should like to point 
out an apparent relation between Pp» and E£ in 
Eq. (1) and report some measurements on the 
water vapor-polyviny! chloride system with and 
without plasticizer from which some significant 
conclusions may be drawn. Finally, it may be 
shown that our data confirm Barrer’s views 
summarized in the preceding paragraph. 


THE APPARENT RELATION BETWEEN P, AND E 


The temperature dependence of the permea- 
bility constant has been determined for a number 
of polymers with water vapor.* From these data 
P, and E in Eq. (1) can be evaluated. In Table I 
these constants are listed for all the materials 
measured. 

In Fig. 1, log Po is plotted versus E. Two 
other measurements of polyvinyl chloride-ace- 
tate, which are mentioned later, are included in 
the graph. The diameters of the circles are rough 
indications of the probable experimental error 
involved. The closeness of fit of the experimental 
data to a straight line, the large range of variables 
covered (ten decades in Py and 18,000 calories 
in E), and the complete randomness with which 
the materials tested were selected clearly indicate 
that in a very general sense log Py has been 
demonstrated to be a linear function of the 
energy of activation. Such a conclusion recalls 
the controversies that have occurred in chemical 
kinetics over apparent relations between A and 
E in the Arrhenius equation (k= Ae-#/"7) for 
reaction rates.‘ A number of these claims were 
shown to be illusory because the rate constant, 
k, was nearly constant in the series of reactions 
considered, thus forcing a change in E to be 
reflected in log A. This situation hardly exists 
here because the permeability constant (corre- 


TABLE I. Energies of activation for permeation and log Po 
for water vapor. 











Polymer E(cal./mole) log Po 
Polystyrene 0 — 6.08 
Polyvinyl chloride-acetate 2350 —5.10, —4.80 
Polyvinyl chloride 4000 —3.51 
Polyethylene (calendered) 8000 — 1.34 
Polyethylene (cast) 10200 +0.13 
Rubber hydrochloride 12800 +1.63 
Polyvinylidene chloride 17500 +4.20 








*C. N. Hinshelwood, The Kinetics of Chemical Change 
(Oxford University Press, New York, 1940). 
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Fic. 1. The correlation of log Po and E (cal./mole) for 
the permeation of water vapor through various polymeric 
materials. 


sponding to k) for the films measured varied 
one thousand-fold. 

If we accept the relation as real, it is necessary 
to attempt an interpretation, consider its range 
of application, and test it on other data if any 
exist. In the few cases in chemical kinetics where 
a bona fide relation between log A and E has been 
found, an explanation involving the frequency 
of vibration in the activated state has been 
postulated. Evidence described in the next 
section indicates that there is no direct relation 
between heat of solution and E. It seems there- 
fore unwise to seek an analogy with this case in 
chemical kinetics. However, it is relevant to note 
that Hammett e¢ a/.5 have demonstrated that in 
a number of cases log A varies only with that 
portion of the energy of activation which is 
kinetic energy. In view of the picture of the 
permeation process described in the Introduction, 
it would seem that, since no bonds are broken in 
the diffusion process, the energy of activation is 
largely kinetic energy, and we may reason as 
follows. Imagine that in some particular reference 
system the permeation was characterized by a 
certain Py and E. On passing then to a system 
with a lower energy of activation E’, it seems 
reasonable, in view of the picture of permeation 
described above, that the lowering of E has been 
due to the simultaneous cooperation of a number 
of factors. That is, in this latter case, more 
degrees of freedom are activated in hole forma- 
tion. A hole is considered to form when several 
segments simultaneously move outward from a 
given region. The translational motion of each 


5L. P. Hammett, Physical Organic Chemistry (McGraw- 
Hill Book Company, Inc., New York, 1940). 
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Fic. 2. The correlation of log Po and E (cal./mole): for 
the permeation of hydrogen through various polymeric 
materials. © = Bakelite; + =ebonite; |+|=Cellophane. 


segment would involve at least three degrees of 
freedom. Each additional degree of freedom 
operative in the latter case can be thought of as 
reducing E by a given amount e, and if m degrees 
of freedom are involved the change is given by 
E-—E’=AE=ne. However, the chance that each 
of these degrees of freedom is favorably activated 
at the moment of activation is considerably less 
than unity. Thus for ” of them to be favorable 
at once the probability of reaction will be reduced 
in the ratio P’/P =a", where a is less than unity 
and is the average probability of one degree of 
freedom being favorably activated at the time 
of the activation. If we take logarithms and 
eliminate n there is obtained 


log P—log P’ = (AE/e) log 1/a 
A log P=AE[(1/e) log 1/a]=RAE. (4) 


Thus it is found that log P will be a linear 
function of E. The validity of the relation just 
derived requires that log (1/a) is directly propor- 
tional to e. That is, the more the energy of 
activation is lowered by the proper activation 
of a certain degree of freedom, the less likely 
(logarithmically) is this activation to occur. This 
concept is statistical in nature: it cannot be 
expected to hold for every degree of freedom 
concerned with permeation in all polymers. 
However, it is expected that the average behavior 
of the various degrees of freedom would conform ; 
consequently, the essential point in this loga- 
rithmic relation is the assumption that the pri- 


mary act of permeation* depends upon a number 
of rarely fulfilled conditions corresponding to a 
large number of degrees of freedom. The con- 
formance of experimental data to the derived 
logarithmic relation lends credence to the as- 
sumption on which the derivation was based. 

It is necessary to note the probable limitations 
of this general behavior. All of the materials 
investigated were pure, unplasticized, linear 
polymers. On the basis of the present information 
nothing can be stated regarding plasticized 
polymers, cross-linked polymers, or polymers 
with which the diffusing gas has a specific 
interaction, such as the water-cellulose system. 

The only other data available for testing this 
apparent relation between Py and E are those of 
Barrer! and deBoer and Fast* on the permeation 
of hydrogen and of Barrer on helium and nitrogen 
through several organic membranes. In Table I] 
values of log Py calculated from these data are 
listed together with the energy of activation for 
all films studied. These values are plotted in 
Figs. 2-4. Unfortunately these data do not 


TABLE II. Energies of activation for permeation and log Po 
for some permanent gases. 











Material E(cal./mole) log Po 
—Hydrogen— 
Bakelite 5100 —6.14 
Celluloid 5600 —4.35 
Nitrocellulose 5700 —4.49 
Ebonite 6500 —4.64 
Cellulose triacetate 7500 —2.85 
Cellophane 7600 —5.23 
Polyethylene 8050 —2.00 
Butadiene-acrylonitrile copolymer 8200 — 1.99 
Polychloroprene 8300 —2.08 
Neoprene (vulcanized) 8300 — 1.65 
—Helium— 
Ebonite 5900 —4.52 
Rubber (vulcanized) 6300 — 3.59 
Polyethylene 7700 —72.44 
Neoprene (vulcanized) 8800 —1.55 
Cellophane 10100 —3.13 
—Nitrogen— 

Butadiene-styrene copolymer 7900 —2.67 
Butadiene-methyl methacrylate 

copolymer 9500 — 1.68 
Butadiene-acrylonitrile copolymer 9800 — 1.89 
Neoprene 10500 —1.27 
Bakelite 10500 —3.01 








* This argument could be presented equally well in terms 
of diffusion. See Discussion. 

6 J. H. deBoer and J. Fast, Rec. Trav. chim. Pays-Bas. 
57, 317 (1938). 
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Fic. 3. The correlation of log Po and E (cal./mole) for 
the permeation of helium through various polymeric 
materials.-++- = ebonite; | +|= Cellophane. 


extend over a large range of the variables and 
hence they cannot be considered decisive. How- 
ever, it is encouraging to note that a linear 
relation is exhibited by all systems measured 
except those of Cellophane, ebonite, and Bake- 
lite. Since none of these materials are pure linear 
polymers, the relation is confirmed with the 
restrictions originally pointed out. Vulcanized 
rubbers also fall in the linear relationship. It thus 
appears as though moderate cross-linking does 
not cause deviations from linear behavior, but 
excessive cross-linking such as is found in Bake- 
lite and ebonite does cause deviations. Conse- 
quently, we may conclude that all available data 
on the temperature dependence of permeability 
of gases in pure polymeric substances conform 
to the equation P = Pye—#/"7 and exhibit a linear 
relationship between log Py and E, except for 
very highly cross-linked materials. 


THE EFFECT OF PLASTICIZATION 
ON PERMEATION 


In order to determine in detail the effect of 
plasticization and composition on the permea- 
bility constant, the diffusion constant, and the 
solubility coefficient, some experiments have 
been performed on pure polyvinyl chloride, a 
polyvinyl chloride-acetate copolymer (Vinylite 
VYNW), and the latter polymer with the two 
different plasticizers added. The permeability 
constants and diffusion constants were measured 
in the manner described in a previous publica- 


tion.’ The data for polyvinyl chloride are taken 
from a previous publication.® 

Polyvinyl chloride was chosen as the pure 
polymer with which comparisons could be made. 
The chloride-acetate copolymer differs from the 
polyvinyl chloride in that five percent of the 
chlorine is replaced by acetate groups. For 
plasticization, two plasticizers were chosen which 
represent rather extreme degrees of interaction 
with water. Tricresyl phosphate is typical of 
plasticizers which are immiscible with water, 
whereas the other plasticizer chosen was a 
mixed glycol ester of 6, 8, and 10 carbon fatty 
acids (SC plasticizer), which is fairly miscible 
with water but is also quite miscible with the 
polymer. 

The results of our measurements with these 
four films at various temperatures are summar- 
ized in Table III. In Figs. 5-8 the logarithms of 
P, D, and S are plotted against reciprocal 
absolute temperature. From the slopes of these 
graphs the energies of activation for permeation 
and diffusion and the heat of solution have been 
calculated and are recorded in Table IV. 

There are several interesting conclusions to be 
obtained from the data as summarized in Table 
IV. First of all, the resolution of P into the 
product of D and S appears to be verified under 
the rather diverse conditions existing in the films 
measured. Moreover, the variation of these 
quantities with temperatures is in all cases of the 
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Fic. 4. The correlation of log Po and E (cal./mole) for 
the permeation of helium through various polymeric 
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7P. M. Doty, W. H. Aiken, and H. Mark, Ind. Eng. 
Chem. 16, 686 (1944). 
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Fic. 5. A plot of logarithm of P, D, and S versus re- 
ciprocal absolute temperature for the permeation of water 
vapor through polyvinyl chloride. 
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Fic. 6. A plot of logarithm of P, D, and S versus re- 
ciprocal absolute temperature for the permeation of water 
vapor through polyvinyl chloride-acetate. 


form previously found for the transmission of 
permanent gases through polymer membranes. 
In other words, the inclusion of a rather large 
amount of solvent into the polymeric system 
does not alter the fundamental nature of the 
permeation. 

Examination of Table IV indicates that E has 
the same value for the two pure polymers, but 
the value is considerably increased by the addi- 
tion of either of the plasticizers. Moreover, this 
increase due to plasticization results almost 
entirely from the change in the heat of solution 
that has taken place. That is, the energy of 
activation for diffusion is relatively constant in 
all four films, but the heat of solution is quite 
sensitive to composition both of the polymer 


itself and that of the added solvent. Similarly, 
P, is greatly increased by the addition of plasti- 
cizer, but Dy is not significantly changed. 

The heat of solution of water vapor in poly- 
vinyl chloride is very large. This probably results 
from the strong interaction of the C—Cl dipoles 
with the OH groups. When some of the chlorine 
is replaced by acetate in the copolymer, the heat 
of mixing is somewhat lowered; the addition of 
25 percent trioctyl phosphate or SC plasticizer 
reduces the heat of mixing very considerably. 
Either the plasticizer molecules solvate the C—C] 
dipoles to a great extent, thus preventing their 
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Fic. 7. A plot of logarithm of P, D, and S versus re- 
ciprocal absolute temperature for the permeation of water 
vapor through chloride-acetate plus 25 percent tricresyl 
phosphate. 
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Fic. 8. A plot of logarithm of P, D, and S versus re- 
ciprocal absolute temperature for the permeation of water 
vapor through chloride-acetate plus 25 percent mixed 
glycol esters of 6, 8, and 10 carbon fatty acids (SC plas- 
ticizer). 
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lasLe III. P, D, and S (calc.) at different temperatures 
for the permeation of water vapor through 
polyvinyl chloride. 











Temperature P X108 DX108 S 
Polyvinyl chloride 
0° 8.58 -232 3.7 
0 10.95 218 5.02 
10 11.5 BY 2.19 
10 12.2 .476 2.56 
25 14.1 1.63 .863 
25 11.9 Be .882 
30 14.9 1.30 .96 
35 15.5 1.76 .88 
45 18.7 2.61 .716 
55 20.3 5.41 375 
Polyvinyl chloride-acetate-copolymer 
0° 16.0 .309 5.18 
10 22.9 , .65 k RY 
25 29.2 1.45 2.01 
38 32.4 2.31 1.40 
Polyvinyl chloride acetate and 25 percent 
tricresyl phosphate 
10° 11.1 1.58 70 
10 9.9 1.4 706 
25 19.7 3.1 .635 
38 42.0 5.14 82 
38 35.2 6.98 504 
55 103.0 16.6 .621 
55 76.0 15.0 .506 
Polyvinyl chloride acetate and 25 percent 
mixed glycol ester 
0° 26.6 15.1 1.76 
5 RY Be 15.3 2.20 
10 44.2 26.2 1.69 
15 47.1 31.9 1.48 
20 70.0 39.7 1.76 








interaction with water, or the plasticizer may be 
so hydrophobic as to almost neutralize the 
hydrophilic nature of the polymer. The former 
assumption seems more probable since both 
hydrophobic and hydrophilic plasticizers increase 
AH by about the same amount. 

The above data may also be used to calculate 
the entropy change which occurs when a mole of 
water vapor dissolves in the polymeric material. 
Thus one can estimate the amount of freedom 


which the dissolved water molecules possess.’ 


The free energy change for the transfer of one 
mole of water from the vapor (at 1-cm Hg 
pressure) to the polymer, is AF= —RT In S. AH 
for the same process is listed in Table IV. AS, 
the change in entropy, is then found from the 
relation, AF=AH—TAS. The entropy changes 
at 20° and 50° are recorded in Table V. 

It can be seen from Table V that there is a 


considerable entropy decrease involved in the 
process of solution of water in polyvinyl chloride. 
Its magnitude indicates that the water molecules 
have lost almost as much-entropy as that in- 
volved when the vapor condenses to ice. They 
must, therefore, be rather tightly held in the 
film, possessing little freedom except when mov- 
ing from one point of attraction to another. In 
the copolymer almost the same conditions pre- 
vail. However, when plasticizer is introduced, 
the entropy change is lowered to that corre- 
sponding to the ordinary solution of a gas in a 
liquid. Thus water molecules in the plasticized 
film are about as free to move as they would be 
in an ordinary liquid solvent. It appears then 
that the effect of the plasticizer has been twofold. 
First, it has greatly lowered the heat of solution 
of water in the film and this in turn has increased 
log Py and E. Second, the freedom of movement 
of the water molecules is enormously increased 
when plasticizer is present. 


DISCUSSION 


If the values of log Py and E£ for the plasticized 
samples of polyvinyl chloride-acetate (Table IV) 
were plotted on the graph in Fig. 1, it is seen 
that these two points would fall above the 
straight line joining the points for pure polymeric 
substances. In considering Barrer’s data for 
permanent gases it was noted that the points for 
Bakelite and ebonite were considerably lower 
than the straight line determined by the other 
materials. The data for Cellophane must be 
excluded from consideration because it is a 
composite film. Since it is a coated regenerated 
cellulose material the actual thickness which is 


TABLE IV. Energies of activation and heats of solution. 








Film material E log Po Ep log Do 4H 


Polyvinyl chloride 2600 —4.96 10,000 —1.17 —7400 
Polyviny! chloride-acetate 2600 —4.62 9,100 -—0.60 —6500 
Polyvinyl! chloride-acetate 





+25% tricresyl phosphate 8700 -—0.30 9,300 -0.72 — 600 
Polyviny]! chloride-acetate 
+25% SC plasticizer 7800 +0.30 7,800 —0.61 0 








used in evaluating the permeability constant is 
meaningless. 

An obvious suggestion may be advanced to 
explain the deviations just summarized, if we 
interpret the behavior as the increase or decrease 
of the energy of activation at constant value of 
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log Po, that is, at a constant value of the entropy 
of activation, AS{. Viewed in this manner the 
addition of small molecules (liquid) merely 
lowers the energy of activation for permeation 
by acting essentially as a lubricant for the matrix 
of polymer chains through which the gas mole- 
cules diffuse. Conversely, excessive cross-linking 
produces such rigidity even on a sub-microscopic 
scale that the attainment of the activated state 
requires a greater energy. 

The explanation offered here for the observed 
relations between log Py and E are admittedly 
not compelling and have clearly been motivated 
by the data and not by a priori reasoning. At 
least two objections could arise. First, a linear 
relation between log Py and E may seem unreal 
since it has been demonstrated that Po is a 
composite quantity equal to the product of Do 
and So and that E is equal to Ep plus AH. If 
Do, So, Ep, and AH were all independent quanti- 
ties, it would be practically impossible to believe 
that a simple linear relation had been found 
between log DoSo and (Ep +AH). This situation 
could be resolved if it could be shown that So 
and AH change very little from system to system 
in comparison to Dy and Ep and hence it is only 
the relation between log Dy and Ep that is 
actually being considered. Since data on these 
quantities are not available for the measurements 
with water vapor, we can only examine the data 
on the permanent gases. Here the AH for the 
various systems vary up to 500 cal./mole and 
log So up to 0.4 from the average. Since the 
points in Figs. 2-4 seem to follow a linear relation 
With just about these tolerances, we can conclude 
that in view of the available experimental data 
the apparent linear relation between log Po and 
E may be a reflection of a real linear relation 
between log Dy and Ep. 

Secondly, in the argument used to obtain a 
relation between log Po and E it may appear 
that the energy of activation of a particular 
degree of freedom will surely differ as the polymer 
is changed, and hence the argument would be 
useful only in explaining a trend—not a relation- 
ship. As previously pointed out this objection 
becomes less and less valid the larger the number 
of degrees of freedom involved, for the deviation 
of the quantity from the average is less signifi- 
cant the larger the number of degrees of freedom. 


DOTY 


TABLE V. Entropies of solution of water in 
various polymers. 











Polymer 20° 50° 
Polyvinyl chloride —24.4 — 24.2 
Polyvinyl chloride-acetate — 20.8 —20.1 
Polyvinyl chloride-acetate+25% 
tricresyl phosphate —3.6 — 3.6 
Polyvinyl chloride-acetate+25% SC 
plasticizer —2.7 —2.7 








It is just this condition of a large number of 
degrees of freedom operating in the activated 
state that Barrer’? has postulated on the basis 
of both kinetic and statistical considerations. 
Using the following equation due to Wheeler,® 
the number of degrees of freedom used in ac- 
quiring energy of activation for diffusion may be 
calculated on the basis of a reasonable model 
assuming the vibration frequency of the diffusing 
particle to be approximated by 2.5 X 10°. 


d?=Dye*-») 





od RT ~ (5 
2 Ep+(n—1)RT is 


d is the mean free path and m is the number of 
degrees of freedom involved. Using data con- 
tained in Table IV this equation can be solved 
for d by assigning various values to . Thus for 
n=6, d=4.5A, and for n=11, d=1.6A, it is to 
be expected that the mean free path of a water 
molecule would certainly be somewhat less than 
the distance between the carbon-chlorine dipoles 
along the polyvinyl chloride chain. Thus a value 
of 1.6A is a good estimate and it thus appears 
that Wheeler’s treatment of diffusion may be 
interpreted as postulating a large number of 
degrees of freedom for activation. 

Wynne-Jones and Eyring? have derived an 
expression for the rate constant which may be 
transformed into the following expression for Do 
with only moderate approximations. 


Do=exp (AS{/R)(RT/h) (d?/2.72). (6) 


If we assume that the mean free path of the 
water molecule in polyvinyl chloride systems is 
1.6A, Eq. (6) can be solved to obtain AS{= 13.9 
e.u. Slightly larger values were obtained by 


( 5a} S. Wheeler, Trans. Nat. Inst. Sci. India 1, 333 
1938). 

9 W. F. Wynne-Jones and H. Eyring, J. Chem. Phys. 3, 
492 (1935). 
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Barrer using 5.0A for the mean free path in the 
case of the permanent gases. Even if these 
estimated values of d are doubled, AS{ is only 
reduced by about 3 e.u. This large value of AST, 
if it were all due to the diffusing molecule, would 
correspond to more than its entropy of solution. 
Since this is not possible, the medium must share 
in the entropy change. Thus the data on the 
diffusion of water vapor corroborate the conclu- 
sions reached by Barrer for permanent gases: 
that activation energy is shared in part or 


wholly in various degrees of freedom of the 
polymeric system. 
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The zeta-potential of vitreous silica in contact with several solutions of barium chloride 
has been determined by the measurement of streaming potentials. Solutions of seven concen- 
trations between 10-* M and 10 M, inclusive, were studied, and the zeta-potential was 
found to decrease in a regular manner with increasing concentration. It is evident that there 
is nO maximum in the curve of the zeta-potential vs. concentration, as most experimenters 


have reported. 


T is the purpose here to report measurements 
of the zeta-potential of vitreous silica in 
contact with solutions of barium chloride. The 
zeta-potential is not directly measurable, but is 
obtained conveniently from measurement of 
streaming potentials arising from the flow of the 
solutions through vitreous silica capillaries. The 
physical significance of the streaming potential 
measurements and their relation through Helm- 
holtz’s equation to the zeta-potential have been 
critically examined elsewhere.! The apparatus 
and technique for obtaining the necessary reliable 
streaming potential data for the computation of 
zeta-potentials have been described previously,’ 
together with measurements of the zeta-potential 
of silica in contact with dilute potassium chloride 
solutions. Added interest is lent to this investi- 
gation by the fact that the data obtained make 





*L. A. Wood, J. Am. Chem. Soc. 68, 432 (1946). 
*Grinnell Jones and L. A. Wood, J. Chem. Phys. 13, 
106 (1945). 


possible a test of Langmuir’s interpretation’ of 
the Jones-Ray effect‘ observed with barium 
chloride. 


EXPERIMENTAL 


The apparatus and experimental technique 
employed in this investigation have been de- 
scribed previously? and no important changes 
have been made. Duplicate sets of streaming 
potential apparatus with capillaries of different 
diameters were operated concurrently so that 
observations with one could be corroborated 
with the other. Previous work has demonstrated 
that initial measurements of the zeta-potential 
are not usually reliable, because the surfaces of 
the capillaries sometimes require many days to 
attain equilibrium with the solution. It has not 


3 J, Langmuir, Science 88, 430 (1938); J. Chem. Phys. 
6, 894 (1938). 

4 Grinnell Jones and W. A. Ray, J. Am. Chem. Soc. 59, 
187 (1937); 63, 288 (1941) ; 63, 3262 (1941); 64, 2744 (1942). 
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TABLE I. Examples of data obtained with 1.00107? M@ 
BaCl, (thirteenth day). 
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TABLE III. Examples of data obtained with 
1.00 10-* M BaCl, (fourteenth day). 









































E R ohme-t E R ohms=1 t 
milli- ohms cm Ex/P milli- P milli- ohms cm! Ex/P milli- 
Flow cm volts x10-6 106 104 volts Flow cm volts x<10-6 106 X10 volts 
Capillary I Capillary I 
LL’ 24.90 33.2 87.71 266.4 «3.55 —34.2 [eL’ 25.20 626.3 833.3 28.04 697 —67.3 
IL’ 24.70 34.0 87.87 = 265.9 3.66 —-35.4 [-+L’ 24.90 623.2 826.5 28.27 7.08 —68.4 
Capillary II Capillary II 
LL’ 24.50 33.1 327.9 265.1 3.58 —-34.6 T->7’ 24.50 680.7 3225 26.95 7.49 —72.3 
I—L’ 2410 32.8 326.8 266.0 3.62 —35.0 Fez’ 24.26 660.7 3145 27.644 7.53 —72.7 
TABLE . eS ne with tration by immersing them in hot nitric acid 
e 2- ° . 

(redistilled) for a day or more. Then steam and 
Capillary I Capillary II hot water from a silica condenser was drawn 
Time devia- dem through them for about an hour, in order to 
days Nos «Xion fmvy ton Xi smv tion remove the nitric acid, followed by a final rinse 
2 : rong’ “a5 a oy ae 2.0% while hot with about 20 cc of the solution about 
4 18 2668 —357 0.6 268.2 —336 12 to be measured. This cleaning treatment was 
5 18 267.9 “we 2 pone “a rs given the capillaries preceding a change in the 
7, = 66.3 34.3 0.6 267.6 33.9 0.6 concentration of the solution. It was found 
: . as me 7 ae 74 “a3 7: effective in establishing conditions for eventual 
10 12 2684 —34.7 0.6 268.2 —340 21 agreement of the zeta-potential of the two 

= : os me 74 a ee er = capillaries within a few millivolts. 
13. 12 2663 —346 17 265.3 —348 0.3 In obtaining streaming potential measure- 








8 “No.” denotes the number of individual measurements constituting 
the series summarized. Since each measurement was duplicated with 
Capillary II, the same number of measurements was made with each 
capillary. Four individual measurements are given in Table I. 

b A new sample of the solution was placed in each set of apparatus 
preceding this measurement. 


been found possible to control the previous 
treatment of the surfaces so that the potentials 
are immediately reproducible, and the practice 
in this laboratory has been to study the variation 


of the zeta-potential with time until it was’ 


reasonably certain that the equilibrium potential 
was reached. The two capillaries were adjudged 
to have reached equilibrium when they reached 
the same value and remained unchanged for 
several days, but in general it was necessary to 
be satisfied with somewhat less than this ideal 
behavior. With solutions of barium chloride the 
approach to a constant value was much more 
rapid than with any of the solutions which have 
been studied in this laboratory at the date of 
writing, including potassium chloride, thorium 
chloride, and lanthanum chloride. 

The capillary surfaces were prepared for 
measuremerts with solutions of a given concen- 


ments, the experimental procedure was to apply 
a measured pressure difference, P, between the 
two reservoirs at either end of the connecting 
capillary, and while the solution was flowing 
through the capillary, the following measure- 
ments were made in rapid succession: the electri- 
cal resistance, R, of the capillary was measured, 
then three measurements of the streaming po- 
tential, E, and of the pressure, P, were made, 
followed by a second resistance measurement. 
Then the direction of flow was reversed under 
approximately the same pressure, and the above 
procedure was repeated. These measurements 
were usually repeated two or three times with 
somewhat different pressures. No dependence of 
the quantity Ex/P on the pressure was ever 
observed. 

In presenting the data obtained, it is impracti- 
cal to give all of the large number of individual 
measurements, but in order to furnish examples 
of the experimental data, a representative ex- 
ample, selected at random, of the data obtained 
with each capillary will be given for each solu- 
tion. Then a table summarizing the results of 
each series will be given. The measurements are 
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TABLE IV. Summary of measurements with 
1.00 10-* M BaCle. 


TABLE V. Examples of data obtained with 
1.00 10-* M BaCl, (eighth day). 

















Capillary I Capillary II kK 
Av. Vv. E R ohms~! 

Time devia- devia- P milli- ohms cm-~! Ex/P milli- 
days No «x X106 ¢ mv tion «X108 g¢mv tion Flow cm volts x10-6 106 x 104 volts 
9> 18 27.80 —67.6 04% 27.38 —72.3 0.4% Capillary I 

10 12 27.80 —67.2 0.3 27.28 —72.2 0.3 I—L’ 23.80 7420 6983 3.346 1043 —100.8 
11 18 27.88 —67.1 0.4 27.62 —73.2 0.3 L<L’ 23.50 7330 7113 «3.285 = 10.25 —99.0 
12 24° 27.98 —67.7 0.4 27.71 —73.9 1.9 

13 18 27.87 —67.4 0.4 27.26 —72.5 1.8 Capillary II 

14 18 28.09 —67.7 0.4 27.43 —72.9 0.8 L—L’ 23.20 6735 24,390 3.564 10.35 —100.0 
154 12 27.85 —66.8 0.7 27.06 —72.0 2.2 L<—L’ 23.10 6395 24,390 3.564 9.87 —95.3 














® See footnote (a) in Table II. 

b Preceding these measurements were earlier series. The first sample 
of the solution was found to give higher values for «x so it was thought 
probable that it had inadvertently become slightly contaminated. 
Therefore a new sample of the solution was placed in the apparatus on 
the eighth day. The zeta-potentials determined with the first sample 
were somewhat lower than those above. 

¢ As a result of experimental difficulties, only 18 measurements were 
obtained for Capillary II in this series. 

4 Following the preceding series, both capillaries were warmed to see 
whether any change in the zeta-potential would result. When there 
was no change, the measurements were discontinued. 


discussed in the chronological order in which 
they were made, except for those with 10-* M 
and 10-* M, which are discussed at the end. The 
temperature was 25°C+0.03° in every instance. 
The values for the fundamental physical con- 
stants are taken from the tabulation of Birge,5 
except the data for the dielectric constant and 
viscosity coefficient, which were assumed to be 
the same as for pure water, as given by Dorsey.® 
The zeta-potential, ¢, is obtained from the 
equation of Helmholtz :!7 


f= (4rnEx)/(DP), 


in which 7 is the viscosity coefficient, «x is the 
specific conductance of the liquid in the capillary, 
and D is the dielectric constant, the other 
quantities being defined above. 


Measurements 


The behavior of the capillary surface with the 
10-* M BaCl, was unusual in that there was 
little variation and almost perfect agreement. 
Two individual measurements with each capil- 
lary are given in Table I. The large reservoirs 
at either end of the capillaries were identified as 
L and L’, and the direction of flow is indicated 





°R. T. Birge, Rev. Mod. Phys. 13, 233 (1941); Am. J. 
Phys. 13, 63 (1945). 

°N. E. Dorsey, Properties of Ordinary Water Substance 
(Reinhold Publishing Corporation, New York, 1940), pp. 
187 and 364. 

7H. Helmholtz, Wied. Ann. 7, 337 (1879). 


TABLE VI. Summary of measurements with 
1.00 10-> M BaCle. 











Capillary I Capillary II 
Av. Av. 
Time devia- devia- 
days No* «x X106 ¢mv tion «x X106 ¢mv tion 
2 18 3.523 —94.3 1.1% 3.697 —90.4 2.3% 
3 12 3.521 —97.3 2.0 3.668 —90.4 1.4 
4 24 3.489 —974 1.4 3.706 —90.5 2.5 
5 18 3.462 —98.3 3.1 3.688 —92.8 3.7 
6 6 3.290 —96.8 2.8 3.477 —88.8 5.1 
7> 24 3.299 —100.8 0.9 3.574. —97.3 1.7 
8 18 3.292 —99.3 0.7 3.556 —96.7 2.1 
9 18 3.306 —97.7 0.8 3.478 —93.6 2.2 
10 18 3.400 -—99.7 1.5 3.681 —97.3 2.3 








® See Note (a) of Table II. 

b Preceding the measurements on the seventh day, the solution in 
both sets of apparatus was mixed by forcing all the solution into 
reservoir L and bubbling a small amount of nitrogen through it. 


accordingly. It is believed that the precision of 
all the measurements is 0.5 percent or better. 
The two potentials given for each capillary are 
the first measurements of two successive series 
of three, as described above, carried out on the 
thirteenth day of observation. Ten other similar 
measurements were made on that day with each 
capillary. All the measurements made with this 
concentration of BaCl, are summarized in Table 
II. It was observed that the zeta-potential of 
Capillary I was initially —37.0 millivolt and on 
long standing decreased slightly to about —34 
millivolt. The zeta-potential of Capillary II was 
—33 millivolt initially, and this increased 
slightly to about —34 millivolt on standing. 
When a fresh solution was substituted after the 
measurements on the tenth day, the potential 
of Capillary I appeared to be undisturbed, while 
that of the other capillary was considerably 
increased. Because the potential of Capillary II 
soon returned to about —34 millivolt, it was 
presumed that the equilibrium had been upset 
for some unknown reason but that it was soon 
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TABLE VII. Examples of the data obtained with 
5.00 10-* M BaCl: (fifth day). 











E R ohms" 
milli- ohms cm! Ex/P milli- 
Flow cm volts 10-6 x<106 x<104 volts 
Capillary I 
LL’ 24.18 105.5 175.00 133.5 5.82 —56.3 
L—L’ 23.86 103.9 175.0 133.5 5.81 — 56.2 
Capillary II 
LL’ 24.66 104.4 650.1 133.7 566 —54.7 
LL’ 24.46 106.2 650.1 133.7 5.80 —56.1 








TABLE VIII. Summary of measurements with 
5.00 X 10-* M BaCl. 











Capillary I Capillary II 
Av. Av. 
Time devia- devia- 
days No.* «X106 ¢ mv tion « X106 ¢ mv tion 
2 18 133.2 -55.1 04% 131.9 —53.8 1.1% 
3 18 1334 —55.7 0.4 132.4 —52.7 1.3 
4 24 133.00 —55.6 0.4 131.5 —55.0> 1.3 
§ 24 1336 —56.1 0.4 133.2 —55.2 1.3 
6 18 133.6 —56.3 0.2 133.0 —55.0 1.6 
7 ® i337 —-3%5 063 133.6 —55.2 1.6 








® See Note (a) of Table II. 

b The increase in the potential of Capillary II on the fourth day may 
be attributed to the fact that, preceding these measurements, the solu- 
tion in reservoir L’ was forced into reservoir L to remix the solution in 
order to eliminate the electrode asymmetry and decrease the average 
deviation. Although this did not result, the disturbance of the double 
layer caused by flow of the solution through the capillary apparently 
had the effect of hastening the attainment of equilibrium. 


reestablished. The accepted value for the zeta- 
potential with this solution was — 34 millivolt. 
With 10-* M BaCl, the observations were 
somewhat less satisfactory. In Table III are 
given individual measurements to serve as ex- 
amples of the measured quantities, selected as 
in Table I. It will be noticed that with both 
capillaries there was a slight asymmetry of the 
observed zeta-potentials, which was the result of 
asymmetry of the silver-silver chloride electrodes. 
Such asymmetries were common occurrences, 
particularly with the more dilute solutions, and 
were ascribed to slight differences in concentra- 
tion which arose in the solutions. These differ- 
ences may have been caused by some of the silver 
chloride dissolving from the electrodes or by the 
glass of the reservoir walls dissolving, or by the 
ingress of some other impurity. The asymmetries 
could usually be eliminated, or greatly dimin- 
ished, by remixing the solution, which was done 
by forcing all the solution into one reservoir and 
bubbling a small amount.of nitrogen through it 
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for stirring. Then half the solution was forced 
back through the capillary to the empty reservoir 
so that the measurements could be continued. 
However, with 10-* M BaCl. the asymmetries 
were not large enough to make this necessary, as 
is evident from the ‘‘average deviations” from 
the mean of each series as given in Table IV, 
the summary of the measurements with this 
concentration. Earlier measurements were made, 
but it was found that the conductivity figures 
were high as compared with those observed with 
the fresh sample of the solution placed in the 
apparatus on the eighth day of observation. 
Therefore the measurements prior to the ninth 
day were discarded. The zeta-potential of Capil- 
lary II was always about 5 millivolts higher than 
that with Capillary I, and there appeared to be 
no tendency for them to converge. However, the 
two values were considered to be in close enough 
agreement so that further work was unwarranted. 
The accepted value for the zeta-potential was 
taken as —70 millivolt, the average on the 
thirteenth and fourteenth days, when it seemed 
that constant values had been reached. 

In Table V are examples of the individual 
measurements with 10-5 M BaCls. In the in- 


TABLE IX. Examples of the data obtained with 
2.50 10-* M BaCl. (second day). 


























E R ohens 1 
Pp milli- ohms cm! Ex/P milli- 
Flow cm volts 10-6 106 x 104 volts 
Capillary I 
L—L’ 27.34 16.8 36.32 643.3 3.95 —38.2 
L<—L’ 27.00 16.9 36.29 643.8 4.03 — 38.9 
Capillary II 
LL’ 28.00 13.6 134.91 644.3 3.13 — 30.2 
L<—L’ 27.80 13.9 134.95 644.1 3.22 —31.1 
TABLE X. Summary of measurements. with 
2.50 10-3 M BaCle. 
Capillary I Capillary II 
Av. Av 
Time devia- devia 
days No* «x X106 ¢mv tion « X106 {mv tion 
1 12 645.0 —38.5 0.3% 643.8 —30.6 2.0% 
2 6 643.7 —38.6 1.0 644.1 —30.7 1.6 
5 6 642.6 —38.6 0.3 642.2 —30.0 0.3 
7> 6 634.2 —37.9 0.5 635.4 —31.2 0.3 
8 12 641.2 —38.1 0.5 640.5 —31.8 2.5 








® See Note (a) in Table II. 
> A new sample of the solution was placed in both sets of apparatus 
preceding these measurements. 
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TABLE XI. Measurements with 1.00X10-? M BaCle.* 








« 
E R ohms~! t 
milli- ohms cm~! Ex/P milli- 
Flow cm volts x<10-% 106 x<104 volts 





A. Examples of the data obtained (fourteenth day) 


A—B* 37.23 34 £13.65 2369 0.216 —2.1 
A+B 40.64 3.2 13.73 2356 0.185 -—1.8 


B. Summary of the data 





t Average 

Time K E/P milli deviation 

days No.> X<106 av. volts millivolts 
1 4 2363 0.114 —2.6 0.62 
14 12 2365 0.081 — 1.84 0.21 
15 12 2369 0.081 — 1.86 0.17 








® The apparatus used for these measurements is that described in 
reference 8. 

b The ‘‘No." is the number of measurements which were averaged to 
yield the figures reported. Two such measurements are given in A. 


stances given, there is a slight asymmetry in 
the values from Capillary I, and a much larger 
asymmetry with Capillary IJ. These were ob- 
served to some extent in all the measurements, 
summarized in Table VI, and were responsible 
for the larger average deviations. Mixing the 
solutions prior to the measurements on the 
seventh day decreased the asymmetry with both 
capillaries, but did not entirely eliminate it with 
either set of apparatus. However, averaging the 
zeta-potentials causes the cancellation of a con- 
stant asymmetry so that the accuracy of the 
measurements is greater than that indicated by 
the average deviation. The accepted equilibrium 
zeta-potential value was —98 millivolt, the 
average of the measurements on the tenth day. 
The next solution measured was the 5X 10-* M@ 
BaCle, and examples of individual measurements 
are given in Table VII. In the examples there 
appeared to be no asymmetry with Capillary I, 
but there was a noticeable amount with Capillary 
II. This is also evident from the figures for the 
average deviations in Table VIII, wherein all 
the measurements with this concentration are 
summarized. Little or no variation with time 
was observed with Capillary I, but there was a 
slight increase in the zeta-potential with Capil- 
lary II. The accepted value of the zeta-potential 
was —55 millivolt. 
Measurements were then made with 2.5 
X10-* M BaCl, (see Table IX), which were 
somewhat more difficult to obtain because of the 
low streaming potential and the difficulty of 
determining the zero point for the resistance 


measurements with the direct current bridge. 
The relatively high conductivity of the solution 
resulted in considerable polarization of the elec- 
trodes when the bridge current was applied, 
causing the zero point to shift. However, the 
measurements of E and R are considered accurate 
to better than 1 percent. Essentially constant 
zeta-potentials were observed with Capillary I, 
all at about —38 mv, a higher value than was 
observed previously with a more dilute solution. 
(See Table X.) With Capillary II, the zeta- 
potential was about —31 millivolt most of the 
time, reaching —32 on the last day of observa- 
tion. Substituting a fresh solution after the 
measurements of the fifth day brought about no 
appreciable change in the zeta-potential of either 
capillary. The accepted value was —35 millivolt, 
the average of the measurements on the eighth 
day. Reference is made below to the fact that 
this is higher than the potential to be expected 
from the measurements with more dilute solu- 
tions. These results were less gratifying than 
those with the other concentrations, in that the 
two capillaries acquired potentials relatively far 
apart, and the accepted value was not as 
concordant with the other measurements as 
might be desired. 

A few measurements were made with 107 M@ 
BaCl, in the simpler apparatus constructed for 
measurements with pure water previously de- 
scribed elsewhere.? The greater ease of manipu- 


TABLE XII. Measurements with 1.0010-*§ M BaCl».* 








A 


E R ohms~! 
P milli- ohms cm Ex/P milli- 
Flow cm volts x10-6 xX 106 104 volts 





A. Examples of the data obtained (third day) 


A+B 18.52 10,007 14,790 2.187 11.81 —114.1 
A—B 15.07 8,013 14,660 2.206 11.73 —113.3 


B. Summary of the data 
Time E/P t Average 





days No. « X106 av. millivolts deviation 
1 24 2.298 520.7 — 115.7 0.6% 
2 24 2.415 504.3 —117.7 0.8 
3> 36 2.250 527.3 —114.5 y& | 
4 24 2.310 513.0 — 114.5 0.6 
5 22 2.315 523.0 —117.1 0.5 
6 12 2.339 504.7 —114.1 0.9 








® The apparatus used for these measurements is that described in 
reference 8. Bright platinum electrodes were used instead of silver- 
silver chloride electrodes to avoid introduction of silver ions and addi- 
tional chloride ions. 

b A new sample of the solution was placed in the apparatus preceding 
these measurements. 


8L. A. Wood, J. Am. Chem. Soc. 68, 437 (1946). 
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Fic. 1. Variation of the zeta-potential with time. 
Data from all-silica apparatus: O. 
Data from Capillary I: (; data from Capillary II: ©. 
A: 10-6 M; B: 10-5 M; C: 10-4 M; D: 5 X10-4 M; E: 10-3 M; 
F: 107 M. 
Selected data for 2.51073 M: @ for Capillary I; @ for 
Capillary II. 


lation of this apparatus, and the fact that the 
low streaming potential limited the accuracy of 
its measurement, rendered the use of this appa- 
ratus more desirable than the more elaborate 
apparatus employed in the above measurements. 
The results of the measurements with this 
concentration are given in Table XI. The con- 
ductance, x, is calculated from R with the cell 
constant, (length)/(cross section), of 32,345. 
The summary of the data indicates that the 
zeta-potential with this solution is about —2 
millivolt. While the percentage deviation is very 
large, because of the low streaming potential, the 
actual magnitude of the average deviation is 
comparable to that of the 10-* M solutions. 
Measurements were made with 10-§ WM BaClas, 
using the same all silica apparatus. With such a 
dilute solution, it is particularly important to 
avoid the introduction of impurities either from 
the glass or the silver-silver chloride electrodes 
of the more elaborate apparatus suitable for use 
with less dilute solutions. Therefore the bright 
platinum electrodes were used, and since the 
liquid came in contact only with these and with 
the silica, this danger was avoided. However, 
the values obtained for the specific conductance 
were larger than would be expected, as canbe 
seen from the data in Table XII. The water used 
in preparing the solutions was conductivity water 
which had been redistilled with a silica condenser, 
and should have had a « of less than 0.5 10~°. 


ROBINSON 


The conductance of 10-* M BaCl, should have 
added about 0.28 X 10-® to this value of « for the 
water, giving a total « of less than 0.8X10-°. 
However, it will be noted that the measured 
resistance of the capillary containing this solution 
was about 14,000 megohms, which gives a «x 
approximately three times larger than 0.8 x 10-°. 
This could be explained in several ways: (1) the 
solution had become contaminated, (2) the extra 
conductance was the result of surface conduct- 
ance, and (3) some path existed for the current 
other than through the solution. If (1) was the 
true reason, then the results are in error, but if 
either (2) or (3) was responsible for the high x, 
then the value for the zeta-potential is correct, 
because the quantity required for the calculation 
of the zeta-potential is in reality the resistance 
between the ends of the capillary, which is the 
quantity measured.! The accepted value of the 
zeta-potential for this concentration is —115 
millivolt, but it must be considered less certain 
than the values for the other solutions. 

In order to show the relative constancy and 
concordance of the measurements, the zeta- 
potentials are plotted in Fig. 1 as functions of 
time, except for the 2.510-* M solution which 
would have been so intermingled with the 10-* M@ 
solution data as to be confusing. The data for 
the seventh day from Table X for the 2.5 
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Fic. 2. Variation of the final zeta-potentials with con- 
centration. Barium chloride: ©; potassium chloride: 9. 
The line W is the zeta-potential of silica in contact with 
pure water, which the curves described by the points should 
approach with dilution. 
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10-3 M solution are shown immediately below 
the letter E, however. The radius of each circle 
is 2.5 millivolts, and it is evident that in general 
the agreement between the values for the two 
capillaries was better than this. 


DISCUSSION 


In Fig. 2 are plotted the accepted values for 
the zeta-potential with each of the barium 
chloride solutions, together with the zeta-po- 
tentials observed earlier? with three solutions of 
potassium chloride. The abscissa represents the 
logarithm of the molar concentration of salt. The 
solid line W is the zeta-potential of silica in 
contact with pure water, known to be about 
—177 millivolt;? both curves should approach 
this line as the concentration approaches zero 
and log c approaches — . Attention is directed 
to the fact that neither with potassium chloride 
nor with barium chloride does there appear to be 
any maximum zeta-potential, and it cannot be 
considered probable that additional data would 
reveal a maximum. This substantiates the ob- 
servations of Rutgers, Verlende, and Moorkens,? 
who reached the same conclusion about the zeta- 
potential of glass in contact with these solu- 
tions. 

Instead of going through a maximum, the 
zeta-potential for these two salts has been ob- 
served to decrease with increasing concentration 
in a more or less regular fashion. It has been 
noted above, however, that the 5107‘, 107%, 


A. J. Rutgers, Trans. Faraday Soc. 36, 69 (1940); 
Rutgers, Verlende, and Moorkens, Proc. Ned. Akad. Wet. 
Amst. 41, 763 (1938). 


2.5X10-* M points are not entirely consistent, 
because it is evident either that the zeta- 
potentials for 5X10-* M and particularly for 
2.5X10-* M BaCl, are too high, or that the 
potential for 10-* M (and possibly 10-? M) are 
too low. It will be recalled that the order in which 
the solutions were measured was 107%, 10~*, 10-5, 
510-4, and 2.510-* M, so that-a considerable 
period of time intervened between the incon- 
sistent measurements. These observations sug- 
gest that the effect may be the result of hyster- 
esis, but it seems improbable that any such 
hysteresis effect could have survived the treat- 
ment with nitric acid and steam which was given 
the capillaries when the concentration was 
changed. It seems more reasonable to conclude 
that this is the result of the inability to control 
completely for a long time the condition of the 
surface of the capillary. The fact that the valves 
for the 2.5X10-* M solution were 6 mv apart 
also indicates that the same equilibrium condi- 
tions did not prevail on each capillary surface 
with this solution. If the values for the 510-4 
and 2.5X10-* M solutions had been 5 millivolts 
lower, and the value for the 10-* M solution had 
been 5 millivolts higher than they were, then 
the points would have described a smooth curve. 
Therefore, it seems evident that the discrepancy 
can be attributed entirely to the inability to 
reproduce the surface, after a long time has 
elapsed, so that the zeta-potential is reproduced 
to better than 5 millivolts of its true value. 
However, in spite of these minor inconsistencies, 
the measurements represent a substantial im- 
provement in the state of present knowledge of 
zeta-potentials. 
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Jones and Ray observed a minimum in the surface tension-concentration curve of electrolyte 
solutions at about 0.001N. Langmuir interpreted this as an instrumental effect rather than 
a real minimum, and advanced a quantitative theory by which the surface tension data could 
be corrected. To compute the Langmuir correction factor the zeta-potentials of vitreous silica 
in contact with the solutions measured by Jones and Ray in their capillarimeter are required. 
_Since the zeta-potentials with barium chloride solutions have now been measured, the Langmuir 
theory has been tested to see whether it eliminated the minimum observed with this solute, 
and this paper presents the results of this test. It was found that Langmuir’s interpretation 


was strikingly substantiated. 





ONES and Ray! measured the capillary rise 
of water and of barium chloride solutions in 
their silica capillarimeter, and found, as with 
other electrolytes, that the first trace of added 
salt brought about an apparent decrease in the 
surface tension. The barium chloride surface 
tension-concentration curve appeared to have a 
minimum at about 0.0005 M, at which the 
surface tension was 99.982 percent of that of 
pure water. Langmuir? suggested that this mini- 
mum was not real but was the result of variation 
in the thickness of the wetting film on the 
capillary wall, and he advanced a quantitative 
theory by which the thickness could be calculated 
provided that the zeta-potential of the silica in 
contact with the various solutions was known. 
Jones and Frizzell* have discussed the application 
of this theory to 1-1 electrolytes, and Jones and 
Wood‘ made the first test of it with data on 
potassium chloride, in which it was found that 
application of the Langmuir correction elimi- 
nated the minimum. Another test® was made 
with data on the same potassium chloride solu- 
tions but with 10-* M thorium chloride added, 
but application of the correction only served to 
deepen the minimum. With the zeta-potential 
of silica in contact with various barium chloride 
solutions reported in the preceding paper, it is 
1 Grinnell Jones and W. A. Ray, J. Am. Chem. Soc. 63, 
288 (1941). 
21. Langmuir, Science 88, 430 (1938); J. Chem. Phys. 6, 
894 (1938). 
— Jones and L. D. Frizzell, J. Chem. Phys. 8, 986 
Grinnell Jones and L. A. Wood, J. Chem. Phys. 13, 


106 (1945). 
5L. A. Wood, J. Chem. Phys. 13, 429 (1945). 


possible to make another test of the Langmuir 
theory, and that is the purpose of this paper. 


LANGMUIR’S METHOD FOR CALCULATING THE 
THICKNESS OF THE WETTING FILM APPLIED 
TO BARIUM CHLORIDE SOLUTIONS 


The method by which Langmuir proposed that 
the wetting film thickness be calculated has been 
previously developed and applied? * 5 to solutions 
of 1-1 electrolytes, but this treatment is not 
applicable directly to the barium chloride solu- 
tions. Therefore, the adaptation to this class of 
electrolytes is given below. 

The Poisson-Boltzmann equation for the elec- 
trical potential in an electrolyte solution in the 
region of a uniformly charged plane is 


dy Agr > (—*) (1) 
Se N42; ’ 
wo oo § 








wherein y is the electrical potential, D is the 
dielectric constant of the solution, n; is the 
number of ions per cc of species 7, 2; is the charge 
on an ion of species 7, and ¢ is the magnitude of 
the electronic charge. Integration of this equation 
yields 


p?=[2A bin; exp (—2:n) J+C, (2) 


wherein 
n=eb/kT, 1/A=)>; N;2;", 


3 
p=dn/d0, 0=Kkx. 8) 


x is the well-known Debye-Hiickel parameter, 
and C is the constant of integration. The con- 
stant C is evaluated through use of the boundary 
condition that at the air surface, where @=0, 
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the potential is assumed to have its minimum 
magnitude y¥,» so that 7=m and p=0. Thus 


d6=dn/[(2A ¥ i ni {exp (—2:n) 


—exp (—z,m)}]'. (4) 
Let 
Bat+ be species 1; 21= +2, 
H+ be species 2; z2= +1, 
Cl- and OH~ be species 3; 23= —1. 


Then substituting K =exp (m) and y=exp (n), 





d 
ite [2A (my?+n i om 
7 ad ne K-?— oK-'—nK) }! 


Since n;=c;NX10-*, where c; is the molar 
concentration of species 7 and N is Avogadro’s 
number, Eq. (5) becomes 


(3c, +2) *dy 


dé= . 
Lary? +azy?+a*y+a, |} 





(6) 
wherein 

0, =C3=2¢,+¢2; 

d2= —[c1K-*+c.K+ (c24+2c:1)K]; = (7) 


Q3=Co,; A4=Ci. 


The constant m introduced into Eq. (4) must be 
evaluated, so that K and az may be known. 
According to Langmuir, if the film is thick 
enough so that van der Waals and other short 
range forces may be neglected, there are then 
the following three forces in equilibrium at the air 
surface of the film: (1) the pressure —(d—£)hg, 
owing to the pull of gravity, which tends to make 
the film thinner (d is the density of the solution ; 
8, the density of air; h, the height above the 
surface of the main body of the liquid; and g, 
the acceleration of gravity); (2) the pressure 
(d—8)hg/2, resulting from the surface tension of 
the cylindrical surface which tends to make the 
film thicker; and (3) the difference between the 
osmotic pressure in the bulk solution and in the 
air interface of the wetting film. From van’t 
Hoff’s simple law and the Maxwell-Boltzmann 
equation, this difference in osmotic pressure is 


Di nikT(K-**—1) 


tending to make the film thicker. At equilibrium, 
the sum of these pressures is zero, and therefore: 


Lisi mikT(K-*—1)—(d—B)hg/2=0. (8) 


Substituting n;=c;N X107: 
6yK-?+¢2K—-!+ (€2.+2¢))K 
500(d—B)hg 

~ NRT 








+3c:+20=—a. (9) 


This equation may be solved for K to give the 
value of wm. Because the potentials to be con- 
sidered are negative and m is therefore negative, 
and since K =e”, the root of Eq. (9) between 0 
and +1 is the only possible value of K. 

The thickness of the wetting film, Ar, is given 


by 
, yi yi dy 
car= f d0=(3cx+es)! f —, (10) 
vo yo Qi 


wherein Q=ay,y’+a2y’+asy+as4, you K, and y1 
is the value of y at the solid surface where y=¢. 
For numerical evaluation of Ar, Eq. (10) must 
be transformed into the Legendrian form of an 
elliptic integral. Let Q be factored such that 
€1;>€2>€3:° , 














Q=a;(y—e1)(y—e2)(y —e@s). (11) 
Substituting 
y=(p+qt)/(1+), 
wherein 
p =e1+((e1—e2)(€i1—es))', 
g =e: — ((€1—e2) (€1—e3))}; 
then 
t dt 
At=—M , (12) 
to [(t?—1)(t?—5?) }} 
wherein 
_ _ (b—e2)(p—es). 
(q—e2)(q—es) | 
unt | 3cit+Ce if 
x Lax(g—e1)(q—e2)(q—es) 
Substituting sin ¢=b/t, sin a=1/b, 
o1 d 
ar= (M/s) [ ———* __. 
¢o [1—sin? asin? ¢ }} 
(13) 


= —(M/b)[ F(a, $0) — F(a, $1) ], 


wherein F(a, @) denotes an elliptic integral of 
the first kind. Although it is not immediately 
apparent, it is true that ¢97=7/2. Thus 


Ar = —(M/b)( F(a, r/2)— F(a, 1) ], (14) 


6 ¢, é2, €3 are numerically equal to the roots of Eq. (9), 
and e2=K for the systems to be considered. 
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TABLE I. Calculation of wetting film thicknesses. 
c 0.0025 0.001 0.0005 0.0001 
1 1.1895 1.3200 1.4880 2.4840 
é2(=K) 0.84865 0.77546 0.70345 0.49077 
es — 0.49522 — 0.48824 —0.47721 — 0.40810 
m —0.16411 —0.25430 —0.35175 —0.71177 
b 2.6349 3.4321 4.4322 10.7772 
a 22.304° 16.939° 13.039° 5.324° 
¥m(mv), —4.21 — 6.53 —9.03 — 18.28 
¢(mv) (experimental) —35 — 34 —55 —70 
"1 0.25596 0.26612 0.11749 0.065517 
ty — 9.6622 — 33.226 — 20.277 — 275.05 
sin $1 — 0.27384 — 0.10330 —0.21859 — 0.039182 
pi ‘ — 15.893° — 5.929° — 12.626° —2.242° 
F(a, /2)— F(a, $1) 1.9103 1.7093 1.8118 1.6129 
1/«(X 108) 35.13 55.54 78.54 175.4 
M(A) 120.5 224.2 Siz.4 1487 
Ar(A) 87 112 152 223 














and the integral is now in a form suitable for 
computation. In some instances the angle ¢; is 
negative, and it should be noted that 


F(a, — 1) = — F(a, $1). 


The negative sign of Ar is a result of the fact 
that the potential is negative. This means that 
although |¢| >|Ym|, ¢<wWm, so that the integra- 
tion must be regarded as performed in the 
negative direction. 


COMPUTATION OF Ar 


The surface tensions of the 2.5X10-*, 110-3, 
5X10-*, and 1X10-* M BaCle solutions were 
measured by Jones and Ray,! and with the 
reliable values of the zeta-potentials reported in 
the preceding paper, their surface tension data 
can be corrected in accordance with Langmuir’s 
suggestion. The Langmuir correction requires 
the thickness of the wetting film at the meniscus 
in the capillary of the surface tension capilla- 
rimeter. The capillary radius was 0.0136 cm, and 
the meniscus was 10.862 cm above the surface 
of the main body of the liquid. In Table I is 
given an outline with numerical values of the 
computation of Ar. With each solution it was 
assumed, although the effect was negligible, that 
the concentration of the ions from the water was 
1X10-* M. The values for the fundamental 
physical constants are those of Birge ;’ the other 
required data are found in the paper of Jones 
and Ray.! The experimental values of the zeta- 
potential were used, rather than those given by 


7R. T. Birge, Rev. Mod. Phys. 13, 233 (1941); Am. J. 
Phys. 13, 63 (1945). 


a smoothed curve. The ¥» values are given to 
show that the condition is fulfilled that these 
are the minimum values, being in each instance 
less than the zeta-potential (in absolute value). 
In Fig. 1 is shown the variation of Ar with 
different values of the zeta-potential. It is evident 
that had smoothed values for the zeta-potentials 
been used, the Ar’s then obtained would have 
differed at most by 5 angstroms from the values 
in Table I. 


THE REFERENCE LIQUID 


The differential method of Jones and Ray 
yields the surface tension values which are rela- 
tive to the surface tension of the reference liquid, 
water, and this observed surface tension is called 
the “apparent relative surface tension.’”’ The 
“corrected relative surface tension” is obtained 
by multiplying the apparent relative surface 
tension by the Langmuir correction factor, 
(r—Ar.)/(r—Aro), wherein Ar, is the wetting film 
thickness of a solution at the meniscus in the 
capillary, and Aro is the wetting film thickness 
with the reference water. According to Langmuir, 
because the value of Arp depends on the concen- 
tration of ions in the reference water, it seems 
probable that the wetting film thickness with 
water should have varied somewhat with the 
purity of the water. However, Jones and Frizzell 
found that there was no variation in capillary 
rise with various samples of water with specific 
conductivities from 0.37X10-§ to 1.67 10-*. If 
the equivalent conductances of the ions in the 
water (presumably H+, HCO;-, OH-) are as- 
sumed to total about 400, this range corresponds 
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to a variation in concentration of 0.9X10-* to 
4.2X10-* M. Actually Jones and Frizzell* gave 
a value for the capillary rise of the 0.9X10-* M@ 
water which was 0.003 percent lower than that 
of the less pure samples, but they apparently did 
not consider this a significant difference, the 
probable error having been about the same 
magnitude. Since the work of Jones and Frizzell, 
it has been learned that the zeta-potential of 
vitreous silica in contact with water of «=0.2 
X10-* to 0.4X10-* is about —177 millivolt.® 
With this knowledge it is possible to reexamine 
the question of whether the Langmuir theory 
would predict a variation for these samples of 
water. Although it is improbable that the zeta- 
potential is —177 millivolt over the entire range, 
it is not likely that the variation is large, since it 
has been shown that the additicn of 10-5 M KCI 
lowers it only to —148 millivolt. It seems justi- 
fiable to assume a value of —177 millivolt over 
this range as an approximation. 

On this basis the wetting film corrections for 
the surface tension data of Jones and Frizzell 
have been calculated as a function of the ion 
concentration of the water. In Table II are 
given the wetting film corrections for four 
different samples of water, including three sam- 
ples measured by Jones and Frizzell, and a 
fourth of lower ion concentration than those. 
The Langmuir correction is given with reference 


TABLE II. Wetting film thicknesses and surface 
tensions for various water samples. 











Langmuir Apparent Corrected 
c Ar(A) correction relative s.t. relative s.t. 
0.54X10-® 366 1.000,11 . 
0.95X10-® 442 1.000,05 0.999,97¢ 1.000,02 
2X 10-6 508 (1.000,000>) (1.000,000) (1.000,000) 
4.2 10-6 543  0.999,97 1.000,00 0.999,97 








* Surface tension not available for this sample. 

>’ This sample is considered the reference liquid for the purposes of 
these calculations. Had any other sample been chosen as reference, the 
same relative values would have been held by all the numbers. 

¢ It has been considered preferable to take this value literally from 
Table VI of Jones and Frizzell (reference 5) and to reserve considerations 
of significance for the ‘corrected relative surface tensions.”’ 


to the 2X10-* M water. It is evident that the 
0.003 percent difference given by Jones and 
Frizzell is in the direction predicted by the 
Langmuir theory, whether or not any significance 


® See Table VI, p. 996 in reference 3. 
*L. A. Wood, J. Am. Chem. Soc. 68, 437 (1946). 
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Fic. 1. Variation of the wetting film thickness in the 
Jones-Ray capillarimeter as a function of the zeta- 
potentials. 


Water: A: c=0.95 X10-*; B: c=2X10-*; C: c=4.2 K1074, 
Barium chloride solutions: D: c =0.0025; E: c =0.001; 
F: c =0.0005; G: c =0.0001. 


may be attached to it. From the column of 
“corrected relative surface tensions’”’ it may be 
concluded that no significant variation should 
have been observed by Jones and Frizzell with 
the various water samples measured. It is also 
evident that had they measured water of the 
lowest ion concentration given, they would have 
detected a significant variation as a result of the 
wetting film, according to the Langmuir theory, 
because the correction factor is large enough to 
indicate that there should have been an ob- 
servable difference, 0.011 percent, easily detect- 
able within the precision of their method, in the 
apparent relative surface tensions. These argu- 
ments are of course all based upon the predicate 
that the true surface tension of the various 
samples of water is the same. It may be pointed 
out that if the zeta-potential decreases somewhat 
with increasing concentration, as it probably 
does, the evidence becomes still more favorable 
to the Langmuir theory, because the Ar’s would 
then be somewhat smaller for the 2X10~* and 
4.2X10-* M samples. It is apparent from these 
considerations that the observations of Jones 
and Frizzell with the various water samples do 
not offer any contradiction to the Langmuir 
theory. It is also apparent that the choice of the 


reference water with which to compare the BaCl, 


surface tensions does not matter greatly, pro- 
vided it is chosen from the range measured by 
Jones and Frizzell. 
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CONCENTRATION 


Fic. 2. Apparent and corrected relative surface tensions. 
Open circles are the measurements of Jones and Ray. 
Solid circles are the same measurements after application 
of the Langmuir correction. Solid line is given by the 
limiting law of Onsager and Samaras. 


APPLICATION OF THE LANGMUIR CORRECTION 
TO THE BaCl, SURFACE TENSIONS 


Whether or not it is proper to place much 
weight on the above considerations as a test of 
the Langmuir theory, a far more impressive 
substantiation can be observed with the barium 
chloride data. In Table III are given the cor- 
rected relative surface tensions of the barium 
chloride solutions with reference to the three 
samples of water. Since the Langmuir correction 
factor depends on the wetting film thickness for 
the reference water, the magnitudes of the factors 
are slightly different, depending on which Arp is 
used. It is assumed in accord with the observa- 
tions of Jones and Frizzell that the values for 
the observed apparent relative surface tensions 
are independent of the reference water, provided 
it is within the range discussed. It is clear from 
Table III that application of the Langmuir 
correction serves to eliminate the minimum from 
the surface tension curve no matter which 
reference Ary is used. The effect of the Langmuir 
correction is shown graphically in Fig. 2, wherein 
the dashed curve with open circles represents 
the observed apparent relative surface tensions, 
and the dashed line with solid circles represents 
the corrected relative surface tensions with refer- 
ence to the 0.95X10-* M water. The radius of 
each circle is 0.001 percent, and it is seen that 


the corrected surface tensions fall almost exactly 
on a straight line. Had the corrected surface 
tensions with reference to the other water sam- 
ples been plotted, the corresponding points for 
each concentration would have been 0.000,05 
and 0.000,08 units above the given ones, for the 
2X10-* and 4.2X10-* M samples, respectively. 
For the general interest, the relative surface 
tensions of these solutions have been computed 
according to the limiting law of Onsager and 
Samaras” and are given by the solid line. In 
these computations, the distance of approach of 
the ions was neglected, but the Y, and Yy series 
of Onsager and Samaras were evaluated to 
enough terms so that no approximation was 
introduced from that source. It may be observed 
that the limiting law curve appears to approach 
tangentially the line described by the corrected 
surface tensions. It is thus apparent that the 
Langmuir correction not only eliminates the 


TABLE III. Corrected surface tensions for 
BaCle solutions. 














Ar Correction Apparent Corrected 
Conc. (angstroms) factor relative s.t. relative s.t. 
Reference water: 4.2 10-§ M@ 
Water 548 (1.000,00) (1.000,00) (1.000,00) 
0.0001 223 1.000,24 0.999,85 1.000,09 
0.0005 152 1.000,29 0.999, 82 1.000,11 
0.0010 112 1.000,32 0.999,84 1.000,16 
0.0025 87 1.000,34 0.999,94 1.000,28 
0.0100 0 1.000,40 1.000,47 1.000,87 
Reference water: 2.0 10-§ M 
Water 508 (1.000,00) (1.000,00) (1.000,00 
0.0001 ° 223 1.000,21 0.999,85 1.000,06 
0.0005 152 1.000,26 0.999, 82 1.000,08 
0.0010 112 1.000,29 0.999, 84 1.000,13 
0.0025 87 1.000,31 0.99994 1.000,25 
0.0100 0 1.000,37 1.000,47 1.000,84 
Reference water: 0.95 10-*§ M 
Water 442 —(1.000,00) —(1.000,00) —(1.000,00) 
0.0001 223 1.000,16 0.999,85 1.000,01 
0.0005 152 1.000,21 0.999, 82 1.000,03 
0.0010 112 1.000,24 0.999, 84 1.000,08 
0.0025 87 1.000,26 0.999 ,94 1.000,20 
0.0100 0 1.000,32 1.000,47 1.000,79 








minimum in the surface tension curve, but yields 
surface tensions which approach infinite dilution 
with a slope in agreement with the limiting law 
of Onsager and Samaras. 


10 L. Onsager and N. N. T. Samaras, J. Chem. Phys. 2, 
528 (1934). 
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Cassie’s statistical treatment of multimolecular adsorp- 
tion on a free surface, leading to the BET equation, is 
criticized and corrected. The method is generalized to 
include adsorption when the maximum number of layers 
of adsorbate is restricted. The isotherm equation obtained 
is the same as that found by Brunauer, Emmett, and 
Teller for this case, rather than the equation proposed by 
Pickett. The BET equation, though not completely satis- 
factory, is apparently thé correct equation for the partic- 
ular model used. Real improvements in the theory should 
follow from refinements in the model rather than from modi- 
fications in the treatment of the BET model. A more refined 


model is proposed and the general method of using it is 
discussed. The BET model is shown to be a crude special 
case of the model suggested here, but it has the advantage 
of presenting no mathematical difficulties. Preliminary 
results, based on improved models, indicate that a some- 
what refined multimolecular adsorption theory is capable 
of predicting capillary condensation for suitable values of 
parameters, thus suggesting that the ideas of multi- 
molecular adsorption and capillary condensation are not 
really in conflict, as is generally assumed. Further discus- 
sion is deferred to a subsequent paper. 





EVERAL papers have appeared recently 

which direct attention to the foundations of 
the Brunauer-Emmett-Teller' theory of multi- 
molecular adsorption. Cassie* has derived the 
BET equation for adsorption on a free surface 
(Eq. (15)) using a statistical argument. Pickett? 
has proposed a modification of the BET equation 
for adsorption when the number of layers of 
adsorbate is restricted (Eq. (27)). Actually, it 
appears to the present writer that Cassie has 
obtained the correct result by an incorrect argu- 
ment and that Pickett’s modification lacks theo- 
retical justification, though it seems to extend 
the range of agreement with experiment in some 
cases. The work of Cassie and of Pickett will be 
discussed further below in the light of results 
obtained here. We consider only localized ad- 
sorption in this paper. 


MULTIMOLECULAR ADSORPTION ON A 
FREE SURFACE 


We adopt the notation of Cassie as far as 
convenient and the BET model. That is, we 
suppose that there are B localized sites per unit 
surface area of adsorbent, that X of these sites 
are occupied by adsorbed molecules (first layer 
of adsorbate), and that A—X molecules have 
been adsorbed on top of the X molecules in the 


1S. Brunauer, P. H. Emmett, and E. Teller, J. Am. 
Chem. Soc. 60, 309 (1938). 

2 A. B. D. Cassie, Trans. Faraday Soc. 41, 450 (1945). 

3G, Pickett, J. Am. Chem. Soc. 67, 1958 (1945). 

‘T. L. Hill, J. Am. Chem. Soc. (March, 1946). 


first layer (thus forming second and higher 
layers). Also, molecules in the first adsorbed 
layer are characterized by a potential energy — €: 
and those in the second and higher layers by — ez 
(assumed the same as in the liquid state), both 
referred to infinite separation as energy zero. 
The partition function for a molecule in the 
first layer is then js exp (e:1/k7) and that for a 
molecule in higher layers is jz exp (€z/kT). One 
can approximate js and jz, by harmonic oscillator 
models :® 


js=(kT/hvs)*j, (1) 
jr=(kT/hv1)* je, (2) 


where j is the partition function for all internal 
degrees of freedom. This approximation is not 
essential to the argument. 

The separate partition functions for the X 
molecules in the first layer and the AX in 
higher layers are then 





B! | i 

Os= Bw s exp (e:/kT) | ’ (3) 

pee cel tan ifr ald 
"a. 


The factor B!/(B—X)!X! is the number of 
distinguishable ways X identical molecules may 
be distributed among 8B sites. The factor 


5 See Eqs. (802,7) and (803,4) of Fowler and Guggen- 
heim, Statistical Thermodynamics (Cambridge University 
Press, London, 1939). 
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(A —1)!/(A—X)!(X—1)! is the number of dis- 
tinguishable ways A —X identical molecules may 
be distributed on top of the X molecules of the 
first layer (that is, the number of ways A—X 
objects may be placed in X boxes, there being 
no restriction on the number per box). 

The complete partition function of the system 
is then 


Q=2 QsQx 
X=1 
A'B! 
(A—X)"(B-X)\(X!)? 
XLis exp (e:/kT) Liz exp (ex/kT) ]4-*; = (5) 
k=A 
k=B 


-5 





(A <B) 
(A>B) 


in which we have neglected unity as compared 
to A and X. The difference between the present 
treatment and that of Cassie should now be 
evident. Essentially, he omitted® the factorial ex- 
pression in Eq. (4) and then incorrectly’ intro- 
duced a factor A!/X!(A —X)! on combining Qs 
and Q,. The two errors happened to cancel each 
other because unity may be neglected compared 
to A and X in Eq. (4). 

We make the usual excellent approximation of 
setting log Q equal to the logarithm of the largest 
term in the sum. The value of X which gives 
this term is found from 


0 log OsQ7 L 


ox (6) 


On taking the logarithm and differentiating as 





indicated, there results 


(A —X)(B—X) =BX?, (7) 
with . 
p= exp [(e,—e:)/AT]. (8) 
Js 


The chemical potential of the adsorbed molecules 
is 
wa _—- 9 log QsQx 
kT 0A 
A-X a1 


=lo ——-lo } 10) 
oy) Tr SIL (10) 


(9) 





where X is to be obtained from Eq. (7). For the 
gas phase, we write 


KG a 
— = — 11 
- re og p, (11) 


where a is a constant which need not be further 
specified for our purposes. Equating wa and ue: 
A-X eé 


a 
——-—log j,=—-+log p. (12) 
y aT Og JL 7 gp. (12) 





log 


For the pure liquid: 


EL 
——-lo == +o 13) 
eT jt tT & Po. ( 
Hence, we have 
x=p/po=(A—X)/A. (14) 
Combining Eqs. (7) and (14) gives the BET 
equation 
‘ A/B=cx/(1—x)(1—x+cx) (15) 
with 
c=1/8. (16) 


MULTIMOLECULAR ADSORPTION WITH A LIMITED NUMBER OF LAYERS 


In this section we shall extend the above treatment to the case of m (the maximum number of 


layers) finite. 


If A molecules are adsorbed, suppose X; are in the first layer, X2 in the second, 


the mth layer. Then 


XitXet+:: 


For given values of X1, Xo, -: 


-+, and X, in 


-+X,=A. (17) 


-, X,, the number of distinguishable arrangements for the second 


and higher layers is found by taking the product of (1) the number of ways X2 objects may be 
distributed among X; sites, (2) the number of ways X3 objects may. be distributed among X, sites, 


6 This resulted from his treatment of second and higher layers as being structureless—the BET model was not 


introduced. 


’ There is no entropy of mixing when identical molecules are exchanged between two energy levels. The argument in the 
introduction of Cassie’s paper should refer to the corresponding factor in Eq. (4). 
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etc. That is, 
X;! X,! _ Xi! 











Xo\(Xi—X2)! Xs!(X2—Xs)!  Xal(Xn1—Xa)! (Xi—Xs)(Xa—Ng)! (Kn —Xa) Xa! (18) 
Then 
B! ss 
O8= Bilis exp (a/kT) Fs (19) 
OL= XAicoe 2 ; (20) 
a FEC IS PO ee 
Q=2 vee p> QsQ1, (21) 


where the limits of summation in Eq. (21) are complicated, but are not needed here. 
We again set log Q equal to the logarithm of the largest term in the sum. The values of Xi, Xe, 
»++, X,_1 which give this term are found from 


0 log QsQ1/0X;=0, 71=1,2,---,n—1, (22) 
using Eq. (17) for X,. One finds, in the same way as Eq. (7) is obtained, 
( (B—X1)Xn=B(X1—X2)(Xn-1—-Xn), 

(X1—X2)Xn=(X2—X3)(Xn-1—-Xn), | 

é (23) 
(Xn-s—Xn-2)Xn=(Xn2—Xna)(Xna Xn), 
(Xn-2—Xn-1) Xn=(Xn-1—- Xn)’. 
Also, from Eggs. (9), (11), (13), and (17) (the steps are the same as for Eq. (14)), 








Xa=x(Xa-1—-X-). (24) 
Combining Eqs. (23) and (24), 
(B—X)x 
X,—-X,=——_———_, 
B 
Xi; -Xiu = Xi- —X; = (X —Xo x! 
1 = ( 1 oe (X, )x (25) 
(B—X))x' p 
=——— (i=2,3, ---,n-—1), 
B 
tame. 
. B 
Finally, by Eqs. (17) and (25), 
A a (X1—X2)+2(Xo—X3) +--+ > +(n—1)(Xn_-1—-Xn) +2X,, (26) 
B (B—X1)+(X1— Xo) +e + (Xn 1 — Xn) + Xn 
B-X, 





(x+2x?+---+nx") 





’ 


e-xX) 
er o<* +x") 
_ ex[1—(n+1)x"+nx"t* ] 
~ (1—x)(1—x+ex—cx*})’ 


which is the BET result. 





(27) 
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oe 


Fic. 1. 


Pickett gives an alternative equation. How- 
ever, the fact that Eq. (27) can be obtained from 
statistical mechanics (considering only possible 
states of the system, without recourse to a 
mechanism), as well as from the kinetic argument 
of the BET theory, indicates that it is certainly 
the correct equation for this model. The obvious 
implication is that real improvements in the 
theory must come from refining the model itself 
rather than from modifications of the type in- 
‘ troduced by Pickett. We’ begin the discussion of 
other models in the next section. 


INTRODUCTION TO A REFINED TREATMENT 


In this section we shall indicate how the BET 
model may be considerably refined, at least in 
principle. A more detailed treatment of the 
subject will be presented in a later paper of 
this series. We shall also treat, in future papers, 
multimolecular adsorption on a mobile first 
layer, the transition from a localized to a mobile 
first layer, and the relations between the adsorp- 
tion isotherm and the two-dimensional equation 
of state (in mobile monomolecular adsorption). 

The most general model to be considered here 
is the following. Molecules are adsorbed onto 
vacant sites of a lattice, with the restriction that 
a particular vacant lattice site cannot be occupied 
by a molecule unless at least one nearest neighbor 
site is already occupied or unless the site is in 
the first layer of the lattice (i.e., in the layer 
next to the surface of the adsorbent). The lattice 
is built up on the surface of the solid and sites 
correspond to possible equilibrium positions of 
molecules in the liquid state. Ordinarily the 
‘packing is assumed to be face-centered cubic, 
with 12 nearest neighbors (z=12). As before, 
the potential energy of interaction between the 


8 See Chapter X of Fowler and Guggenheim (reference 5) 
for an analogous discussion of monolayers. 
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adsorbent and adsorbed molecules in the first 
layer is —e, and the average potential energy of 
a molecule in the liquid state is — ez. Hence, the 
potential energy associated with each pair of 
nearest neighbors is —2€,/z. The total potential 
energy of the system is 


= —¢€,X —2e,N/z, (28) 


where X is the number of adsorbed molecules in 
the first layer and N is the number of pairs of 
nearest neighbors in the particular configuration 
being considered, including pairs involving mole- 
cules in the first layer. 

As an example, in the two-dimensional simple 
cubic configuration of Fig. la, z=4, X =3, and 
N=11. This may be contrasted with a two- 
dimensional BET configuration, illustrated in 
Fig. 1b. The BET model allows (the same is true 
in three dimensions) only vertical additions to 
the lattice and also only takes into account verti- 
cal nearest neighbors. Thus, z=2, N=A-—X, 
and, from Eq. (28), 


U=—6X—€,(A—X). 


In Fig. 1b, X=4 and N=A-—X=7. 

We return now to the general problem. For 
given values of X, A, and B, there are p(X, A, B) 
distinguishable configurations such that each 
molecule has either at least one nearest neighbor 
or is in the first layer, or both. We group these 
configurations according to the value of NV: there 
are g(N, X, A, B) configurations with exactly NV 
pairs of nearest neighbors. Hence, 


Dy g(N, X, A, B)=p(X, A, B). (30) 


We have implicitly been considering the case of 
adsorption on a free surface (n=). If m is 
finite, there is a new obvious restriction on 
acceptable configurations. In this case we may 
write 


Dn G(N, X, A, B,n)=P(X, A, B,n), 
G(N, X, A, B, ~)=g(N, X, A, B), 
P(X, A, B, ~)=p(X, A, B). 
The partition function for the system ("= ©) 
is , 
Q= > {Lis exp (e:/kT) ¥j14~* 
X=1 
X dn g(N, X, A, B) exp (2ezN/zkT)}. 
k=A (A<B), kR=B (A>B). 


(29) 


(31) 


(32) 
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For n finite replace g by G. The general pro- 
cedure is the same as in the previous sections, 
with Eq. (13) still applying. 

The above formal treatment can easily be 
modified so as to apply to capillaries with two 
parallel walls, and with the same lattice sites 
available to molecules building up on either wall. 
This is a more satisfactory but more difficult 
approach than the above (m finite), where, 





essentially, an artificial ceiling is introduced half 
way between the walls. 

The BET model is a very crude special case of 
the present one. Eqs. (5) and (21) follow from 
Eq. (32) if only vertical configurations are 
allowed and only vertical nearest neighbors 
counted (z= 2). Then all configurations (for given 
values of A and X) have the same value of N: 
N=A-X. Also, 








A'B! 
_ A, B)=g(A—X, X, A, B)= 
ia aiailian ) (A —X)!(B—X)!(X!)? (33) 
P(X, A, B, n)=G(A—X, X, A, B,n) 
B! 
a Pa Ta ee eS ee ae i 


Obvious successive improvements in the BET 
model would follow from (1) taking into account 
non-vertical configurations, (2) taking into ac- 
count non-vertical nearest neighbors, and (3) 
using a lattice more suitable than simple cubic. 
To what extent the mathematics can actually be 
carried out in these cases is, of course, another 
question. 

Although the subject will be discussed in more 
detail in a later paper, certain preliminary results 
for capillaries (n small), obtained by the use of 
improved models, may be mentioned briefly. 

If one takes into account, even though crudely, 
interactions between horizontal as well as verti- 
cal nearest neighbors, the resulting somewhat 





refined theory of multimolecular adsorption 
predicts capillary condensation quite naturally 
for suitable values of the parameters. One obtains 
regions of instability, corresponding. to a sudden 
jump in A/B from a small value to a large 
value. The phenomenon is analogous to the 
instabilities in monolayers discussed by Fowler 
and Guggenheim.* The multimolecular adsorp- 
tion and the capillary condensation points of 
view are therefore really not in conflict, as is 
generally assumed, except that Kelvin’s equation 
should not be considered valid on this micro- 
scopic scale. 

The author is indebted to Professor W. A. 
Noyes, Jr. for helpful discussions. 
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By making use of the evaporation-condensation properties of liquid mixtures, the Brunauer- 
Emmett-Teller theory of multimolecular adsorption is extended to mixtures of gases. No 
satisfactory experimental data are available at the present time with which to test the theory. 





INTRODUCTION 


ERHAPS the most important postulate in 

the Brunauer-Emmett-Teller' (BET) theory 
of multimolecular adsorption is that molecules 
of the second and higher adsorbed layers have 
the evaporation-condensation properties of the 
liquid state. If this postulate has really general 
validity (as a good approximation), one should 
be able to treat adsorption from a mixture of 
gases by assuming that molecules of the second 
and higher adsorbed layers in this more compli- 
cated case have the evaporation-condensation 
properties of a liquid mixture at the same 
temperature and composition. Hence, a com- 
parison between experiment and a theory based 
on this assumption would provide a further test 
of the fundamental ideas underlying the BET 
treatment. Such a theory is developed in this 
paper** but, unfortunately, satisfactory experi- 
mental data for comparison are apparently not 
yet available. 


ADSORPTION FROM A MIXTURE OF TWO GASES 


We shall be interested in-the adsorption of two 
gases A and B at relative pressures x4 =pa/pa° 
and xg=pz/pp°, and at a temperature 7 below 
both critical temperatures. Let ga, ds, daa, QBp; 
das, Yea, GANa, and gBNnag be, respectively, the 
heats of adsorption per mole of A on the bare 
surface of the solid, B on the bare surface, A on 
a layer of pure A, B on a layer of pure B, A on 
a layer of pure B, B on a layer of pure A, A on 
a layer of composition N4 (mole fraction of A), 
and B on a layer of composition N,. 

We consider various evaporation-condensation 
equilibria. First, if we have pure liquid A (or B) 

1S. Brunauer, P. H. Emmett, and E. Teller, J. Am. 


Chem. Soc., 60, 309 (1938). 
2@ See also T. L. Hill, J. Chem. Phys., 14, 46 (1946). 


of surface area s in equilibrium with its vapor, 
then 
adaspa°=bas exp (—gaa/RT), 


(1) 
apspp°=bpgs exp (—qse/RT), 


or 
ba 
—=pa° exp (qaa/RT), 
QA 


bp 
—= pp exp (gzzn/RT). 
ag 


Next, if a binary liquid mixture of composition 
Na is in equilibrium with its vapor, we have 
adaSpa=basNa exp (—qAna/RT), 


3 
apspp=bps(1—Na) exp (—qBna/RT). ( 


On substituting Eqs. (2) into (3), one obtains 
ba/pa®=xa=Na exp [(gaa— ana)/RT], 
bs/pe°=xp=(1—Na) 

Xexp [(ges—qana)/RT]. 


(4) 


In studying below the equilibria involved in ad- 
sorption, it will be necessary to have (¢a4—qANa) 
and (gzs—qBNa) as functions of N4. The nature 
of the dependence of these quantities on N4 may 
be found from the experimental vapor pressure 
curves, which can be expressed in the form 


xa=Na/fa(Na), fa(Na)=Na/xa, 
xp=(1—Na)/fa(Na), fa(Na)=(1—Na)/xs. 
That is, the functions fa(Na) and f(Na) may 
be tabulated from the experimental values of 
xa, Xp, and Ny. Comparing Eqs. (4) and (5): 
~ [(qana—qaa)/RT |= fa(Na), 
exp [(98N4—Qse)/RT ]=fa(Na). 


We turn now to adsorption equilibria. In 
addition to the postulates of the BET theory, 





(6) 
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MULTIMOLECULAR ADSORPTION 


we make a further assumption:?> the average 
energy of attachment of a molecule on top of 
the (¢—1)th layer (to become a member of the 
ith layer) depends on the composition of the 
(t—1)th layer only. Let s; represent the surface 
area covered by exactly i layers of adsorbed 
molecules and let Na; be the composition of the 
ith layer. Then, when equilibrium is established 
between all of the layers, we will have 


aaSoPa =basiNA, exp (—qa/RT), 
apSoPp=bpsi(1— N41) exp (—qs/RT), 
QaSi-spa=bas:Na; 
j=2,3, «+ Xexp (—ganaci-1/RT), 
apSi-1Ppp=bgsi(1—Na;) 
Xexp (—qBNaci-1)/RT). 


If we use Eqs. (2) and (6) in Eqs. (7), we can 
write 


(7) 


51NA1=Sox4Qa, 

iw — NA) =SoxpQa, 
S:NAj;=Si-ixafa(NAi-1), 
s(1—NA;) =Si-1%pfp(NaAi-1), 


(8) 
i=2,3,---| 
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where 
Qa=exp [(ga—gaa)/RT], 
Qs=exp [(¢s— ges) /RT }. 


By combining Eqs. (8) in an appropriate manner, 
there results 


$;=Si-1Ei-1 = SoHE Eo: «+ Ei-1, 


(9) 


#=1,2,--- (10) 
Eo=xaQa+xsQz, (11) 
E;=xafa(N4i)+xpfp(Nai), 
oo eer 
{N41=x4Qa/Eo, (12) 


| Na;=axafa(Nais)/Eia, i=2,3,--°. 


Let v4 and v,,4 be the volume of A adsorbed 
arid the volume of A necessary to give a complete 
unimolecular layer, respectively. Then, as in the 
BET theory, we can relate v4/vm4 to the areas 
So, 51, So, -* +. If yA; is the fraction of the area of 
the ith layer (s;+5:::+---) which is made up of 
A molecules, we have 





VA _ Gitset se tsn)yArt (Sot: ++ +5n)yAsts + +SnVAn 


(43) 


Um“ Sotsit:**+5n 


In writing Eq. (13), the number of layers is restricted to n. The case of most interest is n= , 


Substituting Eq. (10) into Eq. (13) gives 





vA _ELA+Ait+biket - ++ +E, Eo: ++ En-1)yAit (Bit: ++ +Ei Ee: ++ En-a)yAat: ++ +21: +En-1vAn] 





Um 4 1+4o(1+2£1+FiE2+:--+EiE2:- 


We assume that the area fraction vy, is related 
to the mole fraction N4, to a sufficient approxi- 
mation, by 
Na 
Na+(1—Na)(0p°/04%)* 


where v4° and v,° are molar volumes. In ordinary 
cases, the error made by replacing the ya; in 
Eq. (14) by the N4; is not serious. For example, 
if v4° and v,° differ by 20 percent, ya and Na 
differ on the average (from N4=0 to N4=1) by 





(15) 


YA 


*b It will be seen from the statistical mechanical dis- 
cussion of the BET theory in the accompanying paper 
that this assumption is the most natural one to make for 
this model. In fact, from the point of view of this statistical 
treatment (where only vertical interactions are taken into 
account), no new postulate is really necessary. That is, 
nas een under discussion follows from the BET 
model, 


Et) (14) 





~about 5 percent. In a more exact treatment of 


any given case, a curve of y4 against N4 should 
be plotted for the appropriate vz°/v4°. Then, in 
using Eq. (14), yA; may be read from the curve, 
for each Na; obtained from Eq. (12). a 
* The quantity vg/vm® may be obtained from } 


_Eq. (14) by replacing the ya; by 1—4i. 


If E21, Eq. (14) does not converge for 
n= 0: 04/Vm4 and vg/tm? become infinite. This 
corresponds physically to condensation on the 
surface of the adsorbent. The value of E,, may 
be obtained easily from 


a xafa(NAw) 
" xafa(N Ae) +xBfa(NAw) 
Ee vars xafa(N Aw) +xpfp(NAqa). 





(16) 


(17) 
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It will be noted that NV, and E,, are independent 
of the nature (Q4 and Qz) of the surface of the 
adsorbent. The influence of the surface dies out 
as the number of layers increases. As one should 
expect, if we put E=1 (condensation) in Eqs. 
(16) and (17), we obtain the vapor pressure 
equations, Eqs. (5). 

There is no difficulty in applying these equa- 
tions when liquid A and B are immiscible over 
part or all of the range in Ny. The experimental 
vapor pressure curves are used as in other cases. 

If xg=0 (or x4=0), Eq. (14) reduces to the 
usual BET equation with ya4;=Na;=1, Ei=x,4 
=x, and Ey=x4Qa=xc. A similar result is ob- 
tained in the very special case: E,;=E,=E; 
=---=F and Naéa,;=NAé.=NaA3=::-=Na. In 
this event, 


VA EvNa 
—= ’ (17a) 
m4 (1—E)(1—E+E£) 


ADSORPTION FROM A MIXTURE OF MORE THAN 
TWO GASES 





Suppose we have a mixture of gases A, B, C, 
-++, We retain the same symbols as above except 
that we use gan, fa(), etc., in which WN repre- 
sents the composition (Nu, Nz, Ne, -:-). The 
extension of Eqs. (3) to (17) to this more general 
case is obvious. In particular, Eqs. (10), (13), 
and (14) are unchanged and Eggs. (11), (12), and 
(15) become, respectively, 


Eo=xaQatxsQat+xcQct:::, 
E;=xafa(Ni) +xefp(Ni) 
txofo(N)+:*+, i=1,2,-> 
NAi=x4Qa/Eo, . 
NAj=Xafa(Ni-nr)/Ev-1, 1=2,3,-°° 
me Na 
~ Na + Na(vp"/va)i+Ne(ve?/va)i+:-: . 
In order to use Eq. (20) in Eq. (14), it is clear 
that the Nz;, Nc;, etc., must also be computed 


as in Eq. (19). However, ya;=N 4; will still be 
a fairly good approximation in most cases. 


(18) 


(19) 


(20) 





YA 


APPROXIMATIONS 
Vapor Pressure Curves 


If the vapor pressure curves of A and B show 
the usual Henry’s law and Raoult’s law char- 


acteristics at small and large mole fractions, 
respectively, it may prove convenient to repre- 
sent the experimental vapor pressure curves by 
simple analytical expressions and use these in 
place of fa(Na) and f(Na). Suitable equations 
for this purpose are® 


x4 =Na/Qas*, B=(1—Naz)? (21) 
“t xp=(1—Na)/Qpa*, a=Na? 
Ww 
Qas=exp [(gas—qaa)/RT], (22) 


Qsa=exp L(gza— ae)/RT]. 


The connection between wag of Fowler and 
Guggenheim® and our symbols is 


exp (—2wap/kT)=QasQza=Q”. (23) 


If Q=1, the solution is ideal and Raoult’s law is 
obeyed over the whole range of composition. 
If Q>1, A and B attract each other more than 
do like molecules, hence x4 and xz are less than 
Na and 1—WN4, respectively. If Q<1, the situa- 
tion is reversed. If Q£1/e?=0.135, two phases 
exist over part of the range in composition and 
Eqs. (21) can no longer be used as satisfactory 
approximations to the true state of affairs. 
In this approximation 


fa(Na)=Qan®, fa(Na)=Qzea*. (24) 

Hence, from Eqs. (11) and (12), 
E;=xaQap®™+xeQpa%™, 71=1,2,°-- (25) 
Na;=xaQap**-/Ei-1, 1=2,3,-°>. (26) 


For more than two gases, 


E;=xaQas?™Qac?© aa ake 


+xpQOpa?Opo---+-++, (27) 


etc., with 
g(A)=Na,*, o(B)=NzB;, etc. 


3 See, for example, R. H. Fowler and E. A. Guggen- 
heim, Statistical Thermodynamics (Cambridge University 
Press, London, 1939), p. 357. By not necessarily requiring, 
in Eqs. (21), that Qag=Qza, it is possible in some cases 
to fit the experimental data better (if Qas=Qzsa, the 
curves for x4 and xg are required to be identical except for 
inversion). This practical advantage is at the expense of 
thermodynamics, however, since Eqs. (21) satisfy 


dlogxa_ ,, dlog xz 
aNa BONs 


only if Qas=Qza. This suggests the policy of using Eqs. 
(21) only when a satisfactory fit may be obtained with 
Qas=Qza, and employing the experimental vapor pressure 
curves otherwise. If this policy is adopted, read Qaa=Qza 
throughout the paper. 
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Approximations of Eq. (14) 


In addition to approximating the experi- 
mental vapor pressure curves, it is highly de- 
sirable to have approximate forms of Eq. (14). 
It should be emphasized that the approximate 
equations given below are in the nature of 
suggestions only, based, however, on the calcu- 
lation of a number of widely different cases 
(using Eq. (24)). It is possible that a more 
complete study would lead to significant im- 
provements. 

The sequences NA;, NAs, NAz, --- and Ey, Es, 
E;3, «++ approach limits Na, and E,, respec- 
tively, given by Eqs. (16) and (17). In most 
cases studied, successive values of N4; oscillate 
above and below Na,, with steadily decreasing 
amplitude, while the E; do not show such 
oscillations.‘ In all cases studied, the E; converge 
much more rapidly than the N4,;. 

The method of approximation used is to break 
off the sums in Eq. (14) at a certain value of 7, 
and use mean values to obtain closed expressions 
representing the rest of the sum. For example, 


14+4£,+4£i:4£24+ Li k2k3+ --- 
.214+414+£1E(1+£+E?*+:--) 


1+£)+FiE2/(1—£), (28) 


where E is some mean value, say, (£3E.)?. 

On the basis of an incomplete study, the 
following scheme is suggested (for m=) in 
order to obtain approximate values which differ 
from those given by Eq. (14) by less than about 
3 percent: 

I. If E, differs from Ez by less than 4 percent: 


Eva 1 
cules a 
¢ tS 


Y=1/(1—E) 
E=E, 


(29) 


(30) 
(31) 


with 
and ya=YAze- 


Il. If £, differs from E2 by from 4 percent to 





about 50 percent: 
Eva Ei 


X=YA A Ag-+—— ]——_,, 
7 +(7 ityéaet —EJizE 


(32) 
Y=1+£,/(1-E), 
with E=(F2E,,)' and 


VA =YVAw(OSC.) 


(33) 


, (34) 
= (yAsY4qw)*(non-osc.). 


III. If £, differs from E, by more than about 
50 percent : 


X=yAait(y4ity42) Ei 
Eva \ 
1—E/1-E' 


+ (vatrastrast 


Y=1+£,+iE2/(1—£) 
with E=(E3E,,)' and 
YA =YAx(OSC.) 
= (yAqyA~)*(non-osc.). 
In the above equations, = «, and 


VA EoX 
= ; (38) 

Um4 1+k,yY 
The mean values suggested are, of course, rather 
arbitrary, but have some justification. Values of 
ya are obtained from the corresponding values 
of Na, using Eq. (15). In Eqs. (34) and (37), 
where two values of y4 are given, the first is to 
be used if the N4; exhibit oscillations of the type 
discussed above. Otherwise, the second value is 
to be used. 

For finite values of m, Eq. (14) itself is only a 
rough approximation since capillary condensa- 
tion is neglected. Also, approximate equations of 
the type given above become rather complicated. 
Hence, it is probably not worth while to use 
approximate equations for finite values of n 
except, perhaps, in case I above. The complete 
Eq. (14) is almost as easy to use in cases II and 
III. For case I, one finds 


Ya Eo{yas(E"*!—E*—E+1)+yal(n—1)E"*!—nE"+E]} 





: Un 
with = 
E=E, 


(39) 


(1—E)[1—E—E,(E"—1)] 


and ya>=Y4e- 


. Oscillations of this type in E£, are, in fact, not possible. For it is not difficult to show, using the thermodynamic 
equation in footnote 3, that E(Na) has a maximum (Q<1) or a minimum (Q>1) at that value of Na which satisfies 
Eq. (16). That is, EZ. is the largest or smallest possible value of E and N4.. is the associated value of Na. 















Fie. 1. 


Since Eq. (14) holds regardless of the number 
of gases in the mixture, the above approximations 
to this equation should have essentially the same 
validity in the general problem of any number 
of gases as they have in the special case of a 
mixture of two gases. 


EXAMPLE 


As far as the author is aware, no experimental 
data are available at the present time which can 
be used to test the theory properly. The following 
experimental results are desirable: (1) vapor 
pressure curves for A and B as functions of N4; 
(2) adsorption isotherms (at the same tempera- 
ture as in the vapor pressure experiments) for 
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pure A and pure B, preferably on a non-porous 
adsorbent (to reduce complications) ; (3) adsorp- 
tion isotherms (at the same temperature and on 
the same adsorbent) of A and B from mixtures 
of gases. Curves of v4/Um4 and vg/Um? should be 
obtained as functions of, say, x4 for a series of 
values of xg. The vapor pressure curves give 
fa(Na) and fe(Na) (or Q) while the adsorption 
isotherms for pure A and pure B can be used to 
obtain Qa, Oz, Um“, and v,2. Then the theoretical 
isotherms for mixtures of A and B may be 
calculated and compared with experiment. For 
more than two gases, the vapor pressures can 
be studied in pairs, as implied in Eq. (27). 

In order to illustrate certain features predicted 
by theory, we have plotted in Fig. 1 isotherms 
for a hypothetical case (using Eq. (24), v4°=v,° 
and m= «). We have chosen a fairly typical case: 
QO=4, Q4=Q2=Q’ =20. In this figure, curve By 
gives Up/Um® as a function of xg for x4 =0. Curve 
B, shows the effect of a constant pressure of A: 
Up/Um® is plotted against xg for x4=.1. At the 
lower pressures of B, the presence of A inhibits 
the adsorption of B, but this effect is reversed at 
higher pressures. Curve A of Fig. 1 represents 
the adsorption of A as a function of xg for 
constant x4 =.1. The adsorption of A at constant 
pressures of A is reduced by the presence of small 
amounts of B but greatly increased for suffici- 
ently large pressures of B. These effects are due 
to Q’ being considerably larger than Q. A and B 
compete for positions in the first layer and the 
molecules of the first layer do not attract addi- 
tional molecules to the extent that the bare 
surface does. Hence, at low pressures, the compe- 
tition reduces adsorption of both A and B. 


CAPILLARY CONDENSATION 


It has been recognized that capillary condensation cannot be treated satisfactorily by merely 
assigning an appropriate finite value to m in the original! BET theory. In a later paper,® Brunauer, 
Deming, Deming, and Teller improved the theory by taking into account the additional energy of 
adsorption in the last adsorbed layer. More recently, Pickett® has introduced a modification of 
the BET theory which appears to be an empirical improvement (for ” finite) but which lacks theo- 
retical justification (in the opinion of the present writer’). 

Before considering mixtures, several corrections to the capillary condensation equations® of the 
BET theory will be suggested. For simplicity we shall not consider the most general of these equa- 


5S. Brunauer, L. S. Deming, W. E. Deming, and E. Teller, J. Am. Chem. Soc. 62, 1723 (1940). 


6G. Pickett,'J. Am. Chem. Soc. 67, 1958 (1945). 
7T. L. Hill, J. Am. Chem. Soc. 68, 535 (1946). 
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tions (Eq. (E)) but the same argument holds. Eq. (D) is derived from 
v _ 2[s1+2s2+- ->+(n—1)Sa-1]+75, 


Um 2(sotsit: ++ +5n—1) +5 


for the case of 2n—1 layers fitting into a capillary. It appears to the writer that the contribution 
of s, to both numerator and denominator has been handled incorrectly and Eq. (40) should read 


v _ 2[s1+252+ -+>+(a—1)5,-1]+(2n—1)5, 





(40) 





























a (41) 
Um 2(sotsit: ++ +5n) 
This gives 
2n—1 2n—1 
caf 1+-( ; g—n x 2ng—g— n-th) ; et] 
Vv 
—= (42) 
Um (1—x)(1+(c—1)x+ (cg —c)x*—cgx™*"] 
instead of Eq. (D). For the case of 2m layers fitting into a capillary, similar reasoning leads to 
v = 2(S1+250+---+N5Sy 
2 Wort 2s ) sa 
Um  2(SotSit---+5n) 
_ex[l +(ng—n)x"— — (2ng—n+1)x"+ngx"t"] (44) 
(1—x)[1+(c—1)x+(cg—c)a"—cgxtt] 
Using the principle of Pickett’s modification, one finds, instead of Eqs. (42) and (44), 
2n—1 2n—1 
ca] 1+-( g—n)ar+ (n—1- z)="| 
v 2 2 
—_ (45) 
- Um (1—x) {1—x+cLx*(g—1)+x]} 
v cxl1+(ng—n)x""'+(n—1—ng)x" 
v _cx[1+(ng—n) ( g)x" | (46) 





Ym  (1—x){1—x+e[x(g—1)+2x]} 


Eq. (46) is thus a correction to Pickett’s Eq. (8). 
In dealing with a mixture of gases (say, two, for example), Eqs. (8) remain unchanged except 
fort=n: 
SnNAn=Sn—1XAafa(NAn-1) Za, 


Sn(1— NAn) =Sn—1XBfp(NAn-1)Z2, 


where g=exp (e€/RT) as in reference 5 (¢€ replacing Q). The quantities e4 and eg will depend here 
on the composition of the (n—1)th layer. The E; and Na; are unchanged except for i=n—1 and 
t=n, respectively. In these cases we write 


En-1 =xagafa(N An) +xpgbfe(NAn-r), (48) 
WAn=Xagafa(NAn—1)/En-1. (49) 


(47) 


and 


Eq. (10) holds except for i=n: 
Sn=Sn—16n-1- (50) 


Then, for 2n—1 layers fitting into a capillary, 
Va 2 (SibSet: ++ +5n)yArt (Set: ++ +5n)yaot: ++ + (GnatSn)VAn—1J+5nVAn 
Um 2(So+sit:+++5n) 


where S, is obtained from Eq. (50) and yA, is given by Eq. (15) using 94,. For 2m layers fitting 
into a capillary, the two factors of 2 in Eq. (51) should be omitted. 





(51) 
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Using Pickett’s modification, 
xagafa(NAn—1) XpZefp(NAn-1) 








6,-1= + (52) 
1—x,4 1—xpz 
and 
xagafa(NAn-1) iii 
An= (53) 
(1 —x4)En-1 


The next problem is the evaluation of the g’s. We can write® 
€4=2kaSaoNa(n—1), €8=2RpSBoN Acn-1), (54) 


where oN 4(n—1 is the surface tension of a solution of A and B of composition Na,_;, and S4 and Sz, 

are the surface areas covered by one mole of A and B, respectively, spread out into a unimolecular 

layer. As an approximation, the constants k4 and kg may be taken as equal to each other. Then, 
€A Sa Vai 


€B Se Va 


(55) 


where V, and V,z are molar volumes. This equation can be used to reduce the number of parameters 
by one, since it gives a relation between ga and gz. A further approximation is to assume ka =kg =1 
and a particular type of packing. If the surface tensions of A and B are known, and a linear relation- 
ship* between surface tension and composition assumed (or the actual relationship measured), 
then ga and gz may be calculated using Eqs. (54). 

The author is indebted to Dr. P. H. Emmett for his helpful comments. 


APPENDIX 


In this section two miscellaneous points which may be of interest will be mentioned briefly. 
1. The logic used by Brunauer, Emmett, and Teller! in deriving the equation 


Uv cx 





—a (56) 
Ym (1—x)(1—x+cx) 
can be improved a little. After obtaining the equation 
v . . . 
—=c>. ini [1 +c> x‘, (57) 
Um i=1 i=1 


it is really not legitimate to rewrite the sums in closed form (as is done in reference 1) without first 
making clear that values of'x greater than unity have no physical significance (at this point in the 
derivation we know only that x is proportional to »). After writing Eq. (57), the next step should 
be either: (1) use the equations (in the BET notation) 


x=(a/b)p exp (E1/RT), (58) 
ni asp®=bs exp (—E1/RT), (59) 
to obtain 

x=p/p%, (60) 


where Eq. (59) is the same as Eq. (1); or (2) point out that in Eq. (57) v is finite for «<1 and infinite 
for x21 and furthermore that we desire v to be finite for p<p° and to be infinite for p2 p°, hence 
x=1 should be identified with p= °. One may then proceed from Eq. (57) to Eq. (56). 

2. In the usual BET theory with n= ~, it is instructive to have some idea of the extent to which 


8N. A. Yajnik, R. K. Sharma, and M. C. Bharadway, J. Indian Chem. Soc. 3, 63 (1926). 
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the surface of the adsorbent is covered as a function of v/v, and c. If we use 
So cx v 
1 -—— =————_- = —(1--x), (61) 
A 1—x+cx Um 


and Eq. (56) (to give x as a function of v/v), we find 


v v\? v } 
Go {sone 
So Um Um Um 








0 sama (62) 
A 2(c—1) 
as the desired expression relating the fraction of the surface covered to c and v/¥m. For v/v» = 1, 
i. (= (c)} 
acura (63) 
A c-1 


and for v/vm= 2, 
So 3c—(c?+8c)! 


A 2(c—1) 





(64) 


Some values of 1—(so/A) are given in Table | for these two cases. It is fairly customary to view the 
process of adsorption, as the pressure is gradually 











_ Faste I. increased from p=0, as taking place essentially 
1 —(0/A) in the following manner: first, the first layer is 
c 2/m =1 0/0m =2 laid down and then, at about v=v,, the second 
001 031 043 layer begins to fill in, etc. It should be clear 
01 .091 .128 from Table | that, according to the BET theory, 
an .240 .333 hi . f the : isfac f 
1 500 667 this picture of the process 1s satis actory or 
A. _ 921 large values of c, but not otherwise. For very 
1000. 969 a small values of c, most of the adsorption when 


v/Um=1 takes place primarily by attaching to 
molecules already adsorbed rather than by at- 
taching to the bare surface. For example, for c=.01, when sufficient molecules have been adsorbed 
to cover the entire surface with a unimolecular layer (v=v,), actually only about one-tenth of the 
surface is covered. Hence, that part of the surface which is covered has, on the average, adsorbed 
molecules about ten layers thick. 
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It is shown that all Bose-Einstein (B.E.) fluid models 
investigated thus far undergo first-order transformations 
in phase space. The ideal fluid condenses smoothly, e.g., 
without any possibility of supersaturation, insofar as its 
thermodynamic characteristic functions and their first 
derivatives, with respect to the independent variables, 
are continuous along the transition or saturation curve. 
The non-ideal B.E. fluids considered, in which attractive 
interatomic forces averaged over the volume of the fluid 
are assumed to operate, undergo sudden condensation 
strictly similar to that exhibited by ordinary fluids in 
coordinate space. This process is accompanied by the 
occurrence of discontinuities of the first derivatives of the 
characteristic functions. By introducing a suitable repul- 
sive interatomic force, smeared over the volume of the 


fluid, the smooth condensation of the ideal fluid is changed 
to a phase transformation of the third order. Here the 
second derivatives of the characteristic functions exhibit 
discontinuities along the transformation curve. The same 
interatomic force lifts the order of the sudden first-order 
transformation by changing it to a second-order one. This 
is accompanied by the appearance of the lambda-point 
type discontinuity of the constant pressure heat capacity 
along the transition or lambda line. The interplay of forces 
of parallel and opposing tendencies in modifying the nature 
of a phase transformation without changing its order and 
in lifting or lowering their order is thus brought out in 
these studies over B.E. fluid models. The bearing of these 
results on the problems connected with the thermal prop- 
erties of liquid helium is touched upon briefly. 





INTRODUCTION 


HE possible intervention of the Bose-Ein- 
stein (B.E.) condensation phenomenon in 

the peculiar transition liquid He I—liquid He II, 
as first suggested by London! has led us to inves- 
tigate different B.E. fluid models. The ideal fluid 
undergoes what may be called a smooth conden- 
sation process in phase space. In this smooth 
transformation, beside the thermodynamic char- 
acteristic functions, their first derivatives with 
respect to the chosen independent variables re- 
main also continuous at the transition. For in- 
stance, the heat capacity at constant volume 
remains continuous at the transition line but has 
a break, i.e., its temperature derivative is discon- 
tinuous. In the modified B.E. fluid models,?* the 
phase changes exhibited were strictly of the ordi- 
nary first-order type. The condensation process 
shown by these fluids resembled strictly that 
undergone by ordinary fluids in coordiate space. 
Since the transformation undergone by liquid He 
is of the second order, these models are even for- 
mally unsatisfactory for its description. It is in- 
tended here to investigate further the B.E. fluid 


* Now with the Federal Telecommunication Labora- 
tories, Inc., New York, New York. 

1F, London, Nature 141, 643 (1938); Phys. Rev. 54, 
947 (1938). 

?F. London, J. Phys. Chem. 43, 49 (1939). 

3 L. Goldstein, J. Chem. Phys. 9, 472 (1941). 


models exhibiting first-order phase changes and 
to clarify and point out some of their properties 
which were not considered before. A modification 
of the assumed properties of these models will 
then be described. This modification has the 
effect of lifting the order of the phase change 
from the first to the second or third order. Apart 
from the theoretical interest of the B.E. fluid 
models which allow a complete and rigorous 
analysis of their phase transformations, the study 
of lifting or lowering the order of phase changes 
may also shed some light on some of the more 
elementary aspects of other cooperative phe- 
nomena. 


I. THE SMOOTH CONDENSATION PROCESS OF 
THE IDEAL B.E. FLUID 


The proof that the ideal B.E. fluid undergoes 
a smooth condensation process is already appar- 
ent on the thermodynamic diagram of this fluid 
in the pressure-volume plane.'* For volumes less 
than the saturation volume, at a given tempera- 
ture, the pressure becomes independent of the 
volume, e.g., the isothermals coincide with the 
isobars in the mixed two-phase or condensation 
region. The smoothness of the condensation of 
the ideal B.E. fluid is also apparent on the shape 
of its isothermals. These reach the saturation or 


4B. Kahn, Dissertation, Utrecht (1938). 
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transition line with a horizontal slope. The iso- 
thermal compressibility of the ideal fluid tends 
in a continuous way toward infinity as the fluid 
approaches saturation and remains infinite in 
the condensation region. It is equally worth 
noting in this connection that the zero slope of 
the isothermals of the non-condensed saturated 
ideal fluid is associated with a maximum of these 
curves. One has indeed for the isothermal curves 
of the fluid of spinless atoms: 


b+ =NkTG(a)/VF(a), V2Vo 


NRT /(T\3 @) 
- (—) G(0)/F), V< Vo. 





J 


Here Vo is the condensation volume at tempera- 
ture T, To is the condensation temperature at 
volume V defined, respectively, by 


Vo=Nh®/F(0)(2amkT)?; 


h? N \3 
To= ( ) . (2) 
2amk \ VF(0) 





In (1) and (2), @ is the negative Gibbs potential : 


per atom of mass m, divided by kT, N is the 
total number of atoms forming the fluid, and 
the functions G(x), F(x) are generalized ¢-func- 
tions; 


G(x) =>> e-**/n5!2; F(x) = —dG/dx; 


E(x) = —dF/dx; D(x) = —dE/dx, etc. G3 
G(0) =¢(5/2) =1.341; F(0)=¢(3/2) =2.612; ) 
E(x)223/x}; D(x)22'/2x?3. 

rK{1 2K 


One finds then easily 


lim (0p,/0V)r=0; 
V-Vo 


(0p_/dV)r<T)=0; 





asp avn. ct, (_NEP\Rl@))D(a) 
fen, @osiorentin (5 ray 


= —1.09NRT/V.', (4) 


which shows that the isothermal curves ~,(V)7 
have a maximum at the saturation line. In other 
words, an extrapolation of the isothermal p,(V)7r 
of the non-condensed fluid into the condensation 
region leads to a physically inadmissible portion 
of the curve along which the pressure decreases 


with decreasing volume. A physical interpreta- 
tion of the vanishing slope of the isothermals 
ps(V)r at the saturation line and the impossi- 
bility of their continuation into the two-phase 
region corresponds to the absence of super- 
saturation in phase space where the condensation 
of the ideal fluid takes place. In regard to this, 
the condensation of the ideal B.E. fluid might 
be compared to what seems to happen with 
ordinary fluids between the temperatures 7,, of 
disappearance of the meniscus and the critical 
temperature 7.5 The experimental evidence for 
the existence of a temperature 7, less than 7, is 
rather scarce.® The experimental results favor in 
the range (7,,, 7.) the same trend of the iso- 
thermal curves of ordinary fluids as that shown 
by the ideal B.E. fluid. The interpretation of the 
condensation of ordinary fluids in the tempera- 
ture interval (7, JT.) corresponds to the absence 
of surface tension. Now, in an ideal B.E. fluid 
where there are no forces acting in coordinate 
space one would naturally expect that surface 
forces be absent and that the condensation in 
phase space should proceed without the inter- 
vention of surface tension. It may be recalled in 
this connection that clusters of 3, 4 or more atoms 
play but a negligible role in the condensation of 
the ideal B.E. fluid. It seems, therefore, that 
there is here an interesting argument favoring 
the interpretation given by Mayer® to the con- 
densation of ordinary fluids in the temperature 
interval (7,,, 7.) provided that further experi- 
mental evidence be accumulated assuring the 
existence of a temperature 7, less than T.. 

This condensation accompanied by the above- 
mentioned behavior of the isothermals may be 
considered as a continuous or smooth first-order 
phase change, in which only the second deriva- 
tives of the characteristic functions of the fluid 
are discontinuous. One has indeed 


A(dp/dV)17,v=Vo=(0p4/dV)7,V=Vo 
—(dp_/dV)7,v=vo=0. (5) 


Similarly, one finds that, C, denoting the con- 
stant volume heat capacity, 


AC, = C.4( Vo) — Cy-( Vo) =0. (6) 


5]. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, Inc., New York, 1940), pp. 295-317. 
6 L. Goldstein, J. Chem. Phys. 10, 295 (1942). 
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= —1.09NkT/V,'. (7) 


We recall here that, E being the total energy of 
the fluid, 





pee ett a alle Pin ~-Med (lt 
Co4/ Wear 4 (a) / (a)—7 (a)/E(a), 
' (8) 
Nk=— = -) G(0)/F(0) 
sad Nk dT Ts 
and, consequently, 
=) (=) 
A =A(— 
dT dT? 
27 Nk (F(a))?D(a) 
= —— — lim (| —————_ 
Ss an (E(a)) 
= —3.66Nk/T, (9) 


using the asymptotic values of D(x) and E(x) as 
given by (3). 

Finally the trend of variation of the constant 
pressure heat capacity is worth noting. One finds 


C Pt 5 nee 3| 
Nk 4 F(a) L (F(a)? 

showing that as the fluid approaches saturation, 
e.g., as a0, Cp, as 1/a! according to the 


asymptotic value of E(a) for a<i. Clearly, in 
the condensation region C,_ remains infinite. 


(10) 





2. THE ORDINARY CONDENSATION PROCESS IN 
PHASE SPACE OF SOME B.E. FLUID MODELS 


In contrast with the preceding behavior of an 
ideal B.E. fluid, the modified fluid models con- 
sidered thus far undergo condensation in phase 
space as do ordinary fluids in coordinate space. 
They exhibit supersaturation in phase space and 
their heat capacity at constant volume has a 
lambda point along the saturation or transition 
curve as is the case with ordinary fluids along 
their vapor saturation curve, at volumes larger 
than their critical volume. Of these modified 
B.E. fluid models the one studied before* seems 
to have the simplest formalism generalizing that 
of the ideal fluid. It is assumed that the influence 
of the (V—1) atoms of the system acting on the 
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particular atom in question can be described by a 
potential energy well independent of the volume 
and temperature. Any given-atom can have one 
or several discrete levels in this well. The 
generality of the model is not reduced, however, 
by assuming that there is but one level —E, 
extending the individual free particle level spec- 
trum of the ideal fluid. The correspondingly 
increased cohesion in coordinate space, as ex- 
pected, hastens the condensation process in phase 
space. One finds thus that, at a given volume, 
the modified fluid starts condensing at a higher 
temperature, while at a given temperature it 
starts condensing at a larger volume than the 
ideal fluid. 

The formulas giving the pressure of this fluid 
above and below the transition temperature, 
denoted here by 7;, are the following: 


o( Wien. 
FP oenat 
NkT kT 
pP+= ’ 7 


V Eo 
F(a’ += 


_NeT oe 
oa F(Eo/kT;)' 





T<?, 


since now the statistical parameter a has to 
satisfy the inequality 


Eo 
a’ =a—— 20. 


2Eo/kT or 
aie kT 


For T2>7;2>T> the family of ee of the 
ideal and modified fluids are identical. 
The saturation curve whose equation is 


pi(V)r5=p(V)7;z 


NkT 
=, G(Es/kT 1)/F(E/kT 9), (12) 


is now displaced to larger volumes. The combined 
thermodynamic diagrams of the ideal and modi- 
fied B.E. fluid models, in the p—v representation, 
are given in Fig. 1. The energy of the system is 
given by 


3 V Eo 
at pa 2amkT TG (+=), T23T; 
(13) 
E_= — Neko; __(2amkT)RTG(E/RT), T<T; 
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where No is the number of condensed particles, or 
No=1/(e*’ —1) 
T\! F(Eo/kT 
-|1-(-) (Eo/ 1 (14) 
T;} F(Eo/kT,) 
and a’~0 for T<T;.* 


The constant volume heat capacitiés are, as a 
simple calculation shows, 


=) F( in) 
Cu. 15 a(a'+ 9 “4 sie 
MS p(wste) (w+ 
C.._ 15 G(E/kT) +3(2 (ry F(Eo/kT) 
Nk 4 F(Eo/kT;)) A T,) F(Eo/kT,) 


4G) sa 
T;} F(Eo/kT;) 




















Therefore, 

AC, (Co—Cy_)T; 

Nk ‘Nk 
9 F(Eo/kTs) -3(—) 
4 E(Eo/kT;) kT; 


2 
- ee E(E/RTs) 
kT;] F(Eo/kT;) 
which is negative and indicates that the constant 
volume heat capacity has a lambda point along 
the transition line. 
The constant pressure heat capacity is also 
discontinuous but its discontinuity is negatively 
infinite. One finds, indeed, that 


Cot _ Cut 25 (G(x))*E(«) _ 15 G(x) 9 F(x) | 
Nk Nk 4 (F(x))8 2 F(x) 4 E(x)’ 





Eo 


=a’+—. (17 
x e+ (17) 


Cp, is thus finite at the saturation line, for 
T=Ty;, but jumps to infinity in the condensation 
region. 


* In reference 3, the energy E, included also the energy 
of the very few atoms which were in the discrete level — Eo. 
However, i in order to obtain heat capacities with the dis- 
continuities predicted by thermodynamics it is necessary 
to omit this energy term in (13). The AC, given in reference 
3 is thus incomplete. 
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This behavior is characteristic to ordinary 
fluids undergoing condensation or, in general, a 
regular first-order phase change. One finds, 
indeed, that in this case and along the transition 


line 
' ; 
AC, «nies a(*) (18) 
dT \oaT/y 


in which C,_ is the constant volume heat capacity 
of the mixed phase, dV /dT is a derivative taken 
along the transition line, and (0p/dT)y is to be 
obtained from the equation of state of the fluid. 
This general thermodynamic relationship, of the 
Ehrenfest type, is indeed satisfied when ordinary 
condensation takes place. Using numerical data 
available for water’? we have computed and con- 
structed c, curves around the vapor saturation 
curve, i.e., for volumes larger than the critical 
volume. Figure 2 gives typical Ac, discontinuities 
for water. It is seen that Ac, decreases with 
increasing temperature. 

It is worth noting here that the modified B.E. 
fluids are the only fluids whose formalism allows 
a rigorous computation of the AC, discontinuity 
and provides thus the only known case where a 
direct statistical checking of the thermodynamic 
formula (18) is at hand. 

It will be shown in the next section that by 
introducing interatomic forces of the repulsive 
type, it is possible to lift the order of the phase 
change of such B.E. fluid models. 


7Cf. N. E. Dorsey, Properties of Ordinary Water Sub- 
stance (Reinhold Publishing Corporation, New York, 
1940). 
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3. SECOND- AND THIRD-ORDER TRANSFORMA- 
TIONS IN PHASE SPACE 


Consider a B.E. fluid model in which the total 
potential energy due to repulsion in coordinate 
space can be represented by a term of the type* 


U(V) =a/V, 


V being the total volume at the disposal of the 
fluid. The constant a will be considered inde- 
pendent of the temperature in a first approxi- 
mation. In this model, the number of individual 
energy levels of spinless atoms associated with 
the momentum range (p, +d) is, at tempera- 
ture 7, 





4nV *d 
dn(p) =— = (19) 
exp +(e) [it |-1 
With 
a eat 


one finds again that at volume V the fluid begins 
to undergo a phase change at temperature JT) 
given by Eq. (2). Similarly, at any given temper- 
ature 7, the fluid undergoes a phase change at a 
volume Vo given also by (2). It can also be 
shown without difficulty that 


piste, 

— av’ we 
aU 

p- “fe Ms 


*It can be seen without difficulty that the results of 
this section are quite general insofar as the potential energy 
function U(V) may contain both attractive and repulsive 
temperature independent terms. The choice of U(V) is 
restricted by the trivial condition that the model so con- 
structed behave in a physically admissible way, i.e., its 
compressibility has to be normal. 
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where pia is the pressure of the ideal fluid. Since 
U(V) is assumed to be independent of the 
temperature, the constant volume heat capacity 
of this non-ideal B.E. fluid is the same as that of 
the ideal fluid. The pressure, however, is still 
increasing with decreasing volume at tempera- 
tures and volumes less than their transition 
values, in the transition region. The introduction 
of a volume dependent repulsive force field of 
the atoms changes completely the nature of the 
phase transition of the originally ideal fluid. The 
modified fluid does not undergo a first-order 
phase change. This is already visible on the 
trend of variation of the isothermal curves in the 
transition region. In order to establish the nature 
of the phase change it is convenient to study 
the constant pressure heat capacity of the model. 

One obtains after a simple calculation, C,, ;, 
denoting the constant volume heat capacity of 


the ideal fluid, 
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It is seen now that, in the first place, the constant 
pressure heat capacity is always finite and con- 
tinuous along the transition line, with a break, 
since (dC,/dT)7) is discontinuous. One has, 
indeed, 


(<2) w por 25 (G/F) 
Nk / 1% Nk /t 8 a/kTV ' 
e. es, 25 (G(0)/F(0))2 
(=).-G Le a/kTwV 


Cp=(Cpi—Cp_)To 
= (Cy, ia — C 


In this model, therefore, both the constant 
volume and constant pressure heat capacities are 
continuous with a break, e.g., their temperature 
derivatives are discontinuous along the transition 
line. The smooth condensation process undergone 
by the ideal B.E. fluid is changed by the re- 














and 


: (22) 
v, ia—~)To=0. 
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pulsive interatomic potential energy smeared 
over the whole volume of the fluid. The transi- 
tion, in phase space, of the modified fluid is of 
the third order. 

The question arises now under what condi- 
tions will a second-order phase change occur, 
e.g., a transition which is accompanied by a 
lambda-point type discontinuity of the constant 
pressure heat capacity of the fluid along the 
transition line. The answer to this question can 
be given at once using fluid models in which the 
forces between atoms are represented by the 
potential energy terms involved in the above 
discussed models. Suppose, indeed, that in a 
B.E. fluid model there is an attractive force 
field represented by the individual negative 
energy level —E» and a repulsive force field 
represented by a term a/V considered above. 
One then finds that such a B.E. fluid model has 
a constant volume heat capacity given by Eq. 
(15), as if the repulsive force field were absent. 
However, the work function of the fluid is now 
increased by the total potential energy term 
Na/V and the pressure is 


P+ = Pro t+ Na/V?, 
p- - Pno-+Na/ vs 


where Pro and p,»— are the pressures above and 
below the transition temperature of that fluid 
in which there are no repulsive forces. The pres- 
sures Pnoy and Pro are given by Eq. (11). It is 
clear from Eq. (23) that this fluid model has no 
coincident isobars and isothermals in the transi- 
tion region. The repulsive forces giving rise to 
the pressure Na/V? have prevented the fluid 
from undergoing a first-order condensation as 
was the case above. 

A long but straightforward calculation leads 
to the following expressions of the constant 
pressure heat capacity of this model: 
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The constant volume heat capacity C,, and C,_ 
are given by Eqs. (15). 

A study of the functions ® and WV shows that 


(Cp_) Ts > (Cos) Tz 
ACp=(Cy4—Cp_)70<0, 


AC, is always finite but negative. 

One has thus the following result. Consider 
any B.E. fluid composed of atoms whose inter- 
action in coordinate space can be represented by 
an average constant attractive potential energy 
term —Ep, and an average repulsive potential 
energy term a/V, smeared over the whole vol- 
ume V of the fluid, and independent of the 
temperature. Such fluid models exhibit a second- 
order phase change in which the constant pres- 
sure heat capacity has a lambda-point discon- 
tinuity along the transition line of the fluid. The 
repulsive force field, averaged over the whole 
volume of the fluid, has lifted the order of the 
phase change. If in the absence of this repulsive 
field the B.E. fluid is ideal and exhibits a smooth 
first-order condensation in phase space, the 
switching on of the repulsive forces causes the 
fluid to undergo instead a third-order transition 
in phase space. This means continuity of the 
energy, entropy, or any other characteristic func- 
tion, with their first derivatives along the transi- 
tion line. The introduction of the repulsive force 
field in a B.E. fluid model undergoing in its 
absence a sudden first-order transformation in 
phase space, similar to that exhibited by ordinary 
fluids in coordinate space, leads to a second-order 
transformation in phase space. This is character- 
ized by the continuity of the characteristic 
functions, which have discontinuous first deriva- 
tives along the transition line. 


(26) 
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The only actual liquid exhibiting a peculiar 
second-order phase change is liquid helium. And 
a possible explanation of the thermal properties 
of this liquid assimilating it to a kind of B.E. 
fluid has been repeatedly suggested by London.” 
The formalism necessary to describe the thermal 
properties of a B.E. fluid undergoing a second- 
order phase change in phase space would have 
to be somewhat similar to the formalism of the 
just discussed B.E. fluid model. However, as 
pointed out already,’ the preceding models be- 
have rather in a normal way below their transi- 
tion line. Their expansion and pressure coeffi- 
cients, (0V/dT), and (0p/dT)y, are positive in 
contrast with those of liquid He II. It is inter- 
esting to recall in this connection that in a recent 
investigation London’ was lead to a qualitative 
argument accounting for this anomalous be- 
havior of liquid He II on the very basis of a 
B.E. fluid model. 


8F, London, J. Chem. Phys. 11, 203 (1943). 
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The difficulties inherent in the theoretical 
treatment of liquids are still further increased 
for a quantum liquid like He II. Here, the 
wave-length associated with the average thermal 
motion of the atoms may become of the same 
order of magnitude as the macroscopic dimen- 
sions of the vessel containing the liquid. This 
led Landau® to build up a continuum type or 
hydrodynamic model for this liquid. It is as yet 
not clear whether the anomalous thermal prop- 
erties of -liquid He II in bulk are correctly 
accounted for by this hydrodynamic model. 

We should like finally to add here that a closer 
analysis of phase changes of higher order ob- 
served in cooperative phenomena might disclose 
also the interplay of forces of opposing tendencies 
as those existing in the B.E. fluid models studied 
here. In particular, the raising or lowering of the 
order of a phase change might be attributed to 
the intervention of forces of opposing nature. 


9L. Landau, J. Phys. USSR 5, 71 (1941). 
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A thorough normal coordinate treatment has been carried out for the non-planar vibrations 
of benzene and its deuterium derivatives. The analysis has led to a critical survey of the assign- 
ments, and has furnished strong independent support for those advanced by Pitzer and Scott. 
It has also furnished values for all eight force constants of the complete harmonic potential 
function. For four of the constants two sets of values were obtained, both of which are physi- 
cally reasonable and which reproduce the frequencies equally well. It was not possible to 
decide between them. The force constants were then used to calculate frequencies for the 
various deuterobenzenes. Nearly sixty such frequencies were evaluated, of which about forty 
could be checked against experiment. With four exceptions the errors are less than 1.5 percent. 
Lastly, certain valence-force constants pertaining to hydrogen vibrations have been determined. 


INTRODUCTION 


ENZENE has been the subject of several 
normal coordinate treatments.'~® The first 
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1E. B. Wilson, Jr., Phys. Rev. 45, 706 (1934). 
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Note added in proof: R. P. Bell [Trans. Far. Soc. 41, 293 
(1945) ] has recently shown how to improve Lord and 
Andrew’s results for the non-planar vibrations of benzene 
while still using a simple valence force system with no 
cross terms and with only two force constants. 


was carried out by Wilson,! who suggested a 
simplified potential function with six force con- 
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stants. This function was applied to the vibra- 
tions of benzene and benzene-ds by Lord and 
Andrews,” but both the function and the data 
with which it was used are now known to be 
inadequate. Meanwhile Manneback‘ has ob- 
tained the most general potential function com- 
patible with the simple harmonic approximation 
and with the Ds, symmetry of benzene. His 
equations contain twenty-six force constants for 
the planar vibrations and eight for the non- 
planar ones. The problem then became one of 
evaluating these various force constants by 
using this potential function with the correct 
set of vibrational assignments. In a second paper 
Manneback® attempted this with the data then 
available, but because the assignments were 
incomplete, and in some cases incorrect, the 
results have little present value. Duchesne and 
Penney* attempted a similar treatment after the 
assignments of Ingold and his collaborators'® 
had been advanced. It too must be regarded as 
unsatisfactory inasmuch as several force con- 
stants had to be dropped from the potential 
function because of insufficient data. For example 
only three of the eight constants for the non- 
planar vibrations could be evaluated. At about 
this same time Bernard, Manneback, and 
Verleysen’ modified Manneback’s original equa- 
tions to apply to the planar vibrations of 
p-benzene-d2, sym-benzene-d3, and p-benzene-d, 
as well as to benzene and benzene-d,. They were 
then used with Langseth and Lord’s" assign- 
ments for the planar vibrations of all five sub- 
stances.* More recently Pitzer and Scott® appear 
to have carried out some calculations with 
Manneback’s equations, but no details are 
given. 

In summary it can be said that Manneback’s 
equations provide the complete harmonic poten- 
tial function for benzene, but that the application 
of these equations to actual data has not been 
thoroughly satisfactory. In the case of the planar 
vibrations the assignments of at least two fre- 
quencies (numbers 14 and 15 in Wilson’s ter- 
minology!) are at best very dubious. There is 
also room for doubt concerning one or two other 
assignments. Any treatment is necessarily in- 


°C. K. Ingold, et al., J. Chem. Soc. 912-987 (1936). 
"A. Langseth and R. C. Lord, Jr., Kgl. Danske Vid. 
Sels. Math. Fys. Medd 16, 6 (1938). 
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complete until these values have been estab- 
lished. No attempt to use the complete potential 
function on the non-planar vibrations has been 
described. 

It, therefore, seemed worth while to attack 
the problem again. A successful treatment should 
provide two useful results. It should furnish a 
set of force constants which are characteristic of 
the benzene ring. One may very possibly be to 
carry these force constants over to various 
derivatives of benzene in much the same manner 
that one already transfers force constants from 
one aliphatic compound to another.” Secondly, 
the analysis will serve as a severe test for the 
validity of the assignments in benzene and its 
several deuterium derivatives. If one uses ob- 
served frequencies in, say, benzene and ben- 
zene-d, to calculate a set of force constants, these 
force constants must give back correct frequen- 
cies for all the various intermediate deuterium 
compounds. If they do not do so, the original 
assignments are in error. 

The obvious starting place is the non-planar 
vibrations, for several reasons. First, the treat- 
ment of the planar vibrations by Manneback 
and his co-workers” ® is quite satisfactory except 
for those cases where the assignments are dubious. 
The need here is for more experimental data. 
On the other hand there has been no thorough- 
going treatment of the non-planar vibrations. 
Second, the non-planar vibrations make the 
most important contributions to the various 
thermodynamic functions because the corre- 
sponding frequencies are relatively low. Hence 
it will be especially useful to have a set of force 
constants for these vibrations. Third, the non- 
planar vibrations split into relatively small 


TABLE I. Symmetry properties of the non-planar vibrations 
of benzene (De, symmetry). 
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symmetry blocks, making the equations easier 
to handle. 

We have carried out a complete normal coor- 
dinate treatment of the non-planar vibrations, 
evaluating all the force constants of the harmonic 
potential function. These constants have been 
used in turn to calculate a large number of 
frequencies for the various deuterobenzenes. The 
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analysis has furnished valuable support for the 
latest assignments, those of Pitzer and Scott. 


MATHEMATICAL METHODS 


There are nine non-planar vibrations in 
benzene, of which three pairs are doubly de- 
generate. Their symmetry properties are shown 
in Table I for Dg, symmetry. We have followed 
Wilson! and Langseth and Lord" in numbering 
the vibrations and in defining the coordinate 
axes. The z axis is perpendicular to the plane of 
the molecule; the y axis is in the plane of the 
molecule and passes through atoms one and four. 

The thirteen possible deuterobenzenes possess 
in all six different symmetries. Each of the sym- 
metry classes of the various lower symmetries is 
formed by a combination of the classes of Dg). 
Fig. 1, which is taken from the paper by Langseth 
and Lord, shows these combinations for the non- 
planar vibrations." The selection rules are also 
indicated, R meaning Raman-active and J infra- 
red-active. In C2, symmetry the twofold axis 
passes through carbon atoms, as in m-benzene-d). 
In Co,* symmetry the twofold axis bisects 
carbon-carbon bonds, as in o-benzene-ds. 

To obtain the secular equation the technique 
of Wilson was used. A valence-force potential 
function which included all the possible second- 
degree interactions was first set up. Two types 
of non-planar valence-force coordinates (VFC) 
were employed. (1) The perpendicular displace- 
ment of the ith hydrogen atom out of the plane 
defined by the two adjacent carbon-carbon bonds 
was termed roy:, where 7o is the equilibrium 








VFC 
Class VFSC 
roy1 Toye Toys roy TOYs Toye Rod Rob: Rod3 Rods Rods Rods 
Au Si 673 6-3 6-3 6-3 673 6-3 
Byg Si 6+ —6-3 6-3 -6-3 6-7 -—6-3 
Ss 6-3 —6-3 6-3 -—6-} 6-3 —673 
Euwst Sita -12) 3 —12-) —12-) 3-% -12-> 
Sita 0 —} 4 0 —} 1 
Ew* Step 4 0 = 4 i 0 —} 
Sin |—374 12-4 12-4 —3-4 7% 12-9 
Ego = S0 0 } } 0 —} —} 
Eg- Sio |—-3-$ —12-% 12-4 oS Uli —12-4 





Fic. 2. U matrix for transforming VFC to VFSC. VFSC =U XVFC. 


18 Langseth and Lord seem to have reversed the ‘‘a” and ‘‘b” designations in their Figure 26, p. 15, for they are not 
consistent with the usage in the rest of the paper. We have taken the liberty of altering them in our reproduction. 
4 FE. B. Wilson, Jr., J. Chem. Phys. 7, 1047 (1939); 9, 76 (1941). 
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Class VFSC Si Ss Ss Si6a Size Sip Sin Sioa Siop 
1,P Q T Y U Zz 
Am Su bats 6 612 6 612 26 
S 48 — 344(1 +49) 
as M M 
re 4 —34(1+4p) (1+40)*, P U Zz r y 
- M M 6 642 246 612 246 
5 16 4(1+3p) 
— M M 
| 
Su 4(1+3p) (1+3e)?, S tT ba U 
» M M 4 344 4 344 
Sve 16 4143p) 
M M 
Est (Symmetric) 
Siz 4(1+3p) (1+3e)? , W UT x 
M M 12 344 12 
(i+p)?, S T 
Eos” Sioa -~ ai 334 
Ey» S: (1+p)? , W 
- - M 12 











Fic. 3. G matrix for non-planar vibrations of the deuterobenzenes. 


M =mass of carbon, 


m =mass of hydrogen (or deuterium), 


p =ro/Ro=C -H equilibrium distance/C —C equilibrium distance. 


P=1/m1+1/m2+1/m3+1/mi+1/ms+1/me, 
QO =1/m —1/me+1/m3 —1/ms+1/ms—1/me, 
S=1/m2+1/m3+1/ms+1/me, 

T = —1/m2+1/ms3 —1/ms+1/me, 

U =1/m2+1/m3 —1/m; —1/me, 


carbon-hydrogen distance and y; is an angle thus 
defined. (2) The out-of-plane ring bending vibra- 
tions were regarded as torsions around carbon- 
carbon bonds. Carbon atoms 1, 2, and 3 define 
a plane; carbon atoms 2, 3, and 4 define another. 
These planes intersect along the C2—C; bond. 
The angle between the planes is designated 8, 
and the corresponding VFC is Rode, where Ro 
is the equilibrium carbon-carbon distance. The 
corresponding coordinate for- the C,;—C, bond is 
designated Rod;, that for the C.—C; bond Rods, 
and so on around the ring. Thus it is evident 
that motions of the hydrogen atoms do not 
affect the Roé;. It is important to note also that 
the Rod; have a different symmetry with respect 
to our y and z axes than do the roy;. None of the 
former lie on the y axis, while two of the latter 
do. In other words the Rod; have the symmetry 


V = —1/m2+1/ms+1/ms—1/me, 

W =4/m +1/m2+1/m3+4/me+1/ms+1/me, 
X =4/m —1/m2+1/ms —4/me+1/ms—1/me, 

Y = —2/m1+1/m2+1/ms3 —2 /ms+1/ms+1/me, 
Z = —2/m —1/me+1/ms+2 /me+1/ms—1/me. 


of the carbon-carbon bonds, while the roy; have 
the symmetry of the atoms. 

Symmetry factoring was obtained by using 
valence-force symmetry coordinates (VFSC), 
which are orthogonal combinations of the VFC. 
The matrix U for the transformation VFSC =U 
XVFC is given in Fig. 2. For brevity the various 
VFSC are designated as Si:, Su, Ss, «++ cor- 
responding to frequencies v11, ¥4, ¥5, ***- 

The G matrix for the VFSC is given in Fig. 3. 
It contains the atomic masses and internuclear 
distances. The outlined blocks indicate the way 
the matrix factors for Ds, symmetry. All the 
terms outside these blocks are interaction terms 
which depend only on the masses of the hydrogen 
(or deuterium) atoms. In evaluating them for 
the various deuterobenzenes it is important that 
the deuterium substituents be properly located 
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TABLE II. Comparison of Manneback’s notation with that 
of this paper. 

Manneback* This paper 
m m 
M M 
c Yo 
a Ro 
k p 
on A 
Bt a 
Bs yg or o® 
B- € 
L 488 
I 344y 
Ba o or ¢? 
N 160 
J 4a 








® Reference 4. 
> Reference 15. 


with respect to the coordinate axes. The atoms 
in the benzene hexagon have been numbered one 
to six around the ring, with atoms one and four 
lying on the y axis. Then my, is the mass of the 
hydrogen (or deuterium) atom in position one 
(in atomic weight units), and so on. All these 
interaction terms go to zero for the two mole- 
cules of Ds, symmetry so that in this case the 
matrix factors as indicated. For lower sym- 
metries this is no longer true. Some of these 
quantities will differ from zero, and the cor- 
responding symmetry blocks of Dg, will be 
connected by these interaction terms. Thus for 
example class A»’’ of D3, symmetry is formed by 
combining classes As, and Be, of De, symmetry 
(cf. Fig. 1). The corresponding interaction term, 
Q/6, will be found unequal to zero for the 
deutero-benzene having D3, symmetry, i.e., 
sym-benzene-d3. Similarly the interaction term 
between E,q+ and E,,~ (or between E,,+ and 
Ew) will also be unequal to zero for sym- 
benzene-d;. All other interactions will be zero. 

There are many possible combinations of this 
type between the classes of De, to give the classes 
of lower symmetries (cf. Fig. 1). It therefore 
seemed preferable to leave the results in the 
form of the F and G matrices instead of reducing 
them to the relatively large number of formal 
equations. One must refer to Wilson’s first paper 
on the subject! for the method of reducing the 
matrices to the secular equations. Wilson’s 
method can be applied quite simply to cases 
where a single interaction term occurs, using an 
observation given in the Appendix. 
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The F matrix for the VFSC is given in Fig. 4. 
It contains force constants only. 

A discussion of the force constants for VFC 
is now in order. The diagonal force constants for 
roy; is termed A, while the interaction constants 
between the various 7oy; are designated ad», ad», 
or d,. The subscripts 0, m, and p refer to the 
ortho-, meta-, and parapositions. Thus a» is the 
constant for interaction between the out-of-plane 
motions of two hydrogen atoms meta to each 
other. Similarly B is the diagonal force constant 
for Rod;, and b,, bm, and 6, are the interaction 
constants between the various Roé;’s. The c’s 
are constants for interactions between the roy,’s 
and the Roé,;’s. The force constants and equi- 
librium distances for C—H and C—D are 
assumed to be identical. 

When the VFC have been transformed to the 
VFSC, the force constants become: 


a=A+2a,+2a,,+a, 
e=A—2a,+2dn—ay | - 
ll ly tig as pavenagen. 
g=A +Go—Gn—Gy 
B=B—2b,+2b,—b, ; 

fad-~i-h44, f™ 
n=2[—Cot¢m—Cp] |Hydrogen-ring interac- 
w=34[c,—Cy |. tion!® (see Table II). 


TABLE III. Frequencies used to evaluate the force 
constants (cm). 











No. Benzene sym-Benzene-d3 Benzene-de 
+ 691 
5 914 

10 850 » Faz 664 

11 671 533 503 

: ° 404+3 373 








% This does not completely agree with Manneback’s 
results (reference 4). He would interchange our o and ¢, 
putting o in class E,~ of De, whereas we have it in E,,*, and 
vice versa. Inasmuch as our results were checked by two 
independent workers, it seems probable that Manneback’s 
Ba and Bs are reversed in his Eqs. (44) and (44’). This 
will cause no error whatever unless one wishes to evaluate 
the valence-force constants from.the Greek constants. 
There is one other error in the second of his Eqs. (44’). 
The first minus sign should be replaced by a multiplication 
sign. 

16 7t may be wondered why Manneback’s notation has 
not been used in this paper. There are two reasons. First, 
four of his force constants differ from ours by a numerical 
factor. Our F or our G matrix would therefore have to be 
altered to compensate for this. Second, there is the dif- 
ficulty of Manneback’s force constants Bs and Ba (foot- 
note 15). It has therefore seemed advisable to use a new 
terminology. Table II compares Manneback’s notation 
with that of this paper. 
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It is interesting to note that the complete poten- 
tial function requires only eight (Greek) force 
constants, but that these are expressed in terms 
of eleven valence-force constants. The Greek 
constants are the ones which are evaluated from 
the experimental data. Consequently we shall 
not always be able to calculate the valence-force 
constants from the Greek constants. One can do 
so for the hydrogen vibrations because there are 
four of each type. For the ring constants and the 
hydrogen-ring interactions it cannot be done. 
One must regard 8, say, as being characteristic 
of a certain type of motion of the benzene ring as 
a whole. It will be impossible to say just how 
much B, b,, bm, or 6, contributes to the value of 
B. We shall therefore have no way of determining 
the relative magnitude of the ortho-, meta-, and 
paravalence-force interactions in this case. 

The problem, then, is to evaluate the eight 
(Greek) force constants from the assigned vibra- 
tional frequencies. Table III lists the assign- 
ments which were used for this purpose. They 
are taken from the paper of Langseth and Lord,"! 
and appear to be reasonably correct. Eight addi- 
tional frequencies for the para-deutero com- 
pounds were also used," but they have not been 
included in the table because they are not essen- 
tial to our results. 

The calculated value of any given force con- 
stant will vary slightly from one deuterobenzene 
to another. There are at least two reasons for 
this. First, the vibrations are not purely har- 
monic as assumed. The anharmonicity decreases 
as more deuterium atoms are introduced into 
the molecule because the average amplitude of 
vibration becomes smaller. Consequently, diag- 
onal force constants will increase with increasing 
deuterium substitution. Whenever possible we 
have used values midway between those for the 
two extremes of light and heavy benzene. A 
second reason for inconstancy is the fact that 
some of the experimental frequencies were 
measured on the liquid and some on the vapor. 
It is well known that a vibrational frequency is 
usually somewhat higher in the gaseous than in 
the liquid state. Neither effect is very marked, 
however. 

The following numerical values were employed 
in the calculations: 


m=miass of hydrogen (or deuterium) = 1.008 (or 2.016) 
atomic weight units. 
M=mass of carbon =i12.01 atomic weight units. 
ro=C—H and C—D equilibrium distance = 1.08A. 
Ro=C—C equilibrium distance = 1.39A. 
p=ro/Ro=0.777. 
v=frequency in wave numbers=const. XA‘, where A 
is a root of the secular equation. When masses are 
in atomic weight units, and force constants are 
10 times their value in dynes/cm, the constant 
is 1302.9. 


EVALUATION OF FORCE CONSTANTS 
a and 6 


The evaluation of a and g¢ scarcely needs 
comment. Frequency 11, which goes with a, is 
the only member of class Ao, in De, symmetry. 
Taking the frequency to be 671 cm™ in benzene 
and 503 cm in benzene-ds, one calculates 
values of 0.2467 and 0.2573, respectively. These 
values differ by 4 percent, which is probably 
somewhat too much to attribute to change in 
anharmonicity alone and may indicate some 
slight experimental error. The average value of 
a is 0.252. 

Four values of @ may be obtained with equal 
ease from vy in benzene and benzene-d, and 
from 104 in p-benzene-dz and p-benzene-d,. The 
following results are obtained: 


benzene v30 = 849.7 ¢g =0.3389 
p-benzene-d: 10a = 850.0 ¢=0.3391 
p-benzene-d, V10a = 662.3 ¢ =0.3405 
benzene-d, v0 = 663.5 ¢=0.3417 


y(average) =0.340 


It is to be remembered that the force constants 
must be multiplied by 10° to convert to dynes/ 
centimeter. 


3, n, and e 


Now that a and ¢ are known, the evaluation 
of B, n, and e can be readily carried out in either 
of two ways. One method is to use the equations 
for class A»”’ or Ds3,. This class contains fre- 
quencies 4, 5, and 11 and force constants a, 8, 
n, and e. (Figs. 1 and 4.) The value of a@ is now 
known, so that the three frequencies will suffice 
to determine the three unknown force constants. 
The assignment of 691, 914, and 533 cm to the 
three vibrations is reasonably certain. On car- 
rying out the calculation one finds an interesting 
result; two sets of values are obtained for the 
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Fic. 4. F matrix for non-planar vibrations of 
the deuterobenzenes. 


force constants (Table IV). The value for e, 
which is the constant for the hydrogen vibration, 
is unique, but there is no obvious means of 
choosing between the two sets of values for 6 
and 7. The difficulty arises because 8 (or 7) is 
obtained from a quadratic equation, and both 
roots of the quadratic are physically reasonable. 

One can also evaluate 8, 7, and e from class 
B;, of V, symmetry, which contains frequencies 
4, 5, and 100 and force constants ¢, 8, , and e. 
Using Langseth and Lord’s assignments of 634, 
966, and 738 cm™ for p-benzene-d2 and of 605, 
927, and 765.3 for p-benzene-d,, and using the 
known value of gy, one gets again two sets of 
solutions as indicated in Table IV. 

It is apparent that sets 1 and 3 are com- 
parable, and sets 2 and 4. The differences be- 
tween comparable values are large, and may be 
due to any of the following reasons: (a) an error 
in yg or, more probably, in a@ (incidentally the 
average values of these constants were used); 
(b) a small error in the assignments for the three 
compounds; (c) the fact that the frequencies for 
sym-benzene-d3; are vapor-phase frequencies, 
while those for the para derivatives are liquid- 
phase frequencies. (This, however, probably 
should make the diagonal force constants e and 
8 higher for the former than for the latter, since 


F. A. MILLER AND B. L. CRAWFORD, 


JR. 


frequencies are usually higher in the vapor state. 
The converse is found.) It should be mentioned 
that sets 1 and 2 give identical results when 
used to calculate frequencies, and so do sets 3 
and 4. One must then decide which particular 
set to use, or whether to take the average of, say, 
sets 1 and 3. We arbitrarily chose to use set 1 in 
later calculations, as it seemed to give somewhat 
better agreement with observed results than did 
either set 3 or the average of sets 1 and 3. 

One is now in a position to calculate v4 and »; 
in benzene. These frequencies have never been 
assigned with certainty. Pitzer and Scott 
selected the weak Raman line at 685 cm~ for v4, 
and then used the product rule to calculate »; 
as 1016 cm~. We have calculated v4 and v5; from 
our values for e, 8, and n, using both sets 1 and 3. 
The results are: 


Set 1 Set 3 
V4 693 669 
V5 993 1041 


It will be observed that the calculated values for 
v, bracket the weak line at 685 cm—. This offers 
strong support for the choice of Pitzer and Scott, 
and makes an otherwise weak assignment very 
probable. We then follow these authors in taking 
vs to be about 1015 cm™. 


6, o, and w 


The obvious place from which to calculate 
these force constants is class E” of D3, symmetry, 
which contains frequencies 10, 16, and 17 and 
force constants ¢, 0, ¢, and w. The value of ¢ has 
already been found, so that one has three fre- 
quencies with which to evaluate three unknown 
force constants. Langseth and Lord assigned the 
values 711.7, 373, and 815 cm to vio, vig, and 
viz, respectively, in sym-benzene-d;. Unfor- 
tunately these authors are uncertain about the 
assignment for viz, and we shall see later that 
it must be revised upward by about 100 cm™. 
But using their values for the moment, one can 


TABLE IV. Values for e, 8, and 7. 











From D; From Va 
1 2 3 a 
€ 0.532 0.532 0.585 0.585 
B 0.203 0.339 0.205 0.367 
n 0.258 0.373 0.279 0.415 
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TABLE V. Location of 947. 











Trial Freq. No. Bz(calc.)  $-Bz-ds(assumed) 3Bz-de(calc.) 

10 -— 711.7 _ 

1 16 398 373 350 
17 1002 950 805 
10 — 711.7 — 

2 16 403 373 346 
17 958 920 789 
10 — 711.7 — 

3 16 409 373 342 
17 912 890 . 772 








calculate some tentative force constants. Once 
again two sets of results are obtained. Both sets 
reproduce Langseth and Lord’s assignments in 
m-benzene-dz, sym-benzene-d3, vic-benzene-d3, 
and m-benzene-d,, but they do not give back the 
observed frequencies of benzene and benzene-d 
nor of benzene-ds; and benzene-ds. Thus for 
example instead of getting 404 cm for v6 in 
benzene, one calculates 433 cm=. It must. be 
concluded that Langseth and Lord’s assignments 
are not internally consistent through this series 
of isotopic molecules. 

Zeigler and Andrews” were led by a study of 
the heat capacity of solid benzene-d, to suggest 
that Langseth and Lord’s value for v7 was too 
low. Pitzer and Scott® have discussed the assign- 
ment at some length. They say that 815 cm™ 
for v17 in sym-benzene-d; leads to an unreason- 
able potential function, and also that it leads to 
a calculated heat capacity for benzene which is 
too high. By a consideration of weak bands they 
assign 985 cm in benzene and 790 cm in 
benzene-ds. From these values they find v7; to 
be 950 cm-! in sym-benzene-d3. 

We arrive at essentially the same results, but 
by a different procedure. There are three relevant 
cases where one can be fairly certain that rio 
and vy. are correctly assigned, namely 712 and 
373 cm=! in sym-benzene-d; and 404 cm in 
benzene. The correct set of force constants must 
reproduce these three values. Only two of these 
frequencies belong to sym-benzene-d3, however. 
Our procedure was to estimate a value for the 
third frequency in sym-benzene-d; (v7), and then 
to use these three frequencies to calculate a trial 





1” W. Zeigler and D. H. Andrews, J. Am. Chem. Soc. 
64, 2482 (1942). 


set of values for 6, ¢, and w. These in turn were 
then used to calculate v1, and »;; in benzene, with 
the hope that they would give 404 cm for 746. 
It soon became apparent that in order to obtain 
404 cm—, one has to choose 717 for sym-benzene-d; 
within the range 890-950 cm~. The results for 
three such trials are indicated in Table V. Ob- 
viously v7 is confined within the following 
ranges: (a) 890-950 cm= in sym-benzene-d;; 
(6b) 910-1000 cm= in benzene; (c) 770-805 cm 
in benzene-ds. One cannot locate »:7 more pre- 
cisely in this way because the calculated value 
of v1¢ in benzene is rather insensitive to changes 
of v7 in sym-benzene-d3. Furthermore the nu- 
merical value of vig in benzene has an uncer- 
tainty of about 3 cm—.!8 

The next step is to look in these regions of the 
spectra for weak bands that may be attributed 
to v17. It may be noted that v0a, vise, and vi7¢ 
will have the same numerical values in the fol- 
lowing groups of compounds because the cor- 
responding portions of the G matrices are 
identical: (a) benzene, benzene-d;, and p-ben- 
zene-d2; (b) m-benzene-d2, sym-benzene-d3, vic- 
benzene-d3, m-benzene-d,; (c) p-benzene-d,, ben- 
zene-ds, benzene-ds. Thus the search may be 
extended to the spectra of these additional 
compounds. 

Starting with benzene, one finds a very weak 
infra-red band in the vapor at 962 cm-, and in 
the liquid at 985 cm-!. Inasmuch as 77 is for- 
bidden to appear in benzene, and as selection 
rules are rigorous in the vapor phase, the first of 
these cannot be accepted. The second would be 
a possibility if it really belonged to a vibration 
different from the 962 cm one, and if liquid 
forces were breaking down the selection rules. 
This 985 cm! band appears as a shoulder on 
the strong 1033-cm— absorption, so that it is 
difficult to estimate the frequency and intensity 
accurately. One concludes that the spectrum of 
benzene offers little help. In benzene-ds a very 
weak line is observed at 790 cm in the Raman 
spectrum of the liquid. If this is v7, it appears 
in violation of the selection rules. Pitzer and 
Scott were the first to assign 985 and 790 to v7. 
Going now to the intermediate isotopic deriva- 
tives, we find a possibility in the weak Raman 


18 See p. 928 of reference 10. 
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TABLE VI. Values of the force constants 
(kX 10-* dynes/cm).* 











1 2 
a 0.252 _ 
7) 0.340 _- 
€ 0.532 -—4 
B 0.203 0.3398 
n 0.258 0.373° 
o 0.402 — 
0 0.157 0.757 
w — 0.154 —0.515 








* Set 1 was used in most of the calculations. 
a See also Table IV. 


line at 920 cm™ in vic-benzene-d3, observed by 
Langseth and Lord. They had assigned it to 
vp, but questioned their choice. There are no 
other observed frequencies in these ranges than 
these three (with the exception of the very 
intense breathing frequency »;, which cannot be 
mistaken). We therefore take 1; to be 920 cm 
in sym-benzene-d;, and 790 cm in benzene-dg. 
In benzene it will be somewhere near 985 cm~", 
the uncertainty being due to its appearance as a 
shoulder on another strong band. 

The justification for this assignment will be 
found in its results. Using 711.7, 373, and 920 
cm! for vio, vig, and 47 in sym-benzene-d3, and 
taking g=0.340, one calculates the following 
two sets of values for the force constants. 


1 2 
o 0.402 0.402 
0 0.157 0.757 
w — 0.154 —0.515 


Again there is no obvious reason for preferring 
one set over the other. We arbitrarily chose set 1 
for later calculations. 

Using these force constants, one calculates v17 
to be 789 cm in’ benzene-dg, in good agreement 
with the observed 790 cm-!. For benzene, one 
calculates 958 cm in place of the (rather 
indefinite) 985 cm. Pitzer and Scott have 
already pointed out that this assignment explains 
the rather strong infra-red bands in benzene at 
1810 and 1965 cm™, which appear in both gas 
and liquid, as the allowed combination tones 
Viotvi7 and vg+717. A value of about 955 cm 
is indicated for 17 from these considerations. 

This assignment of v1; is further strengthened 
by some results for p-benzene-d2. Frequencies 
16b and 17) become infra-red active in this 
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derivative. Their values were calculated from the 
force constants o, 6, and w to be 360 and 868 
cm~!, respectively. A sample of p-benzene-d, 
was prepared and the infra-red spectrum of the 
vapor obtained.!® A strong band was observed at 
872 cm—, which may very well be 717. This offers 
valuable support for the assignments and for the 
values of the force constants. 

Thus we have arrived at virtually the same 
conclusions concerning v1; as did Pitzer and 
Scott, but by quite a different method. Indeed, 
our assignments for all the out-of-plane vibra- 
tions are essentially the same as theirs. 

This completes the determination of the force 
constants. For convenience they have been 
gathered together in Table VI. Where two or 
more sets of values were obtained, the first set 
was the one which was used in our calculations. 


FURTHER RESULTS 


Once all the force constants are at hand, one 
can proceed to calculate many of the non-planar 
frequencies for the various deuterium deriva- 
tives. This has been done for those cases where 
the symmetry blocks are not too large to make 
the labor prohibitive. Table VII gives the 
results. Thirty-nine of the fifty-nine calculated 
frequencies may be checked against experiment. 
The maximum error is 4.3 percent or 16 cm™. 
With but four exceptions, calculated and ob- 
served values agree within 1.5 percent. The 
results may, therefore, be considered satis- 
factory. 

One is also now able to evaluate the valence- 
force constants which are associated with the 
hydrogen vibrations—i.e., A, @o, Gm, and @p. 
Introducing the numerical values of a, ¢€, ¢o, 
and g, one calculates A =0.378, a,=—0.057, 
Qn =0.007, and a,= —0.026. It is interesting to 
note that the force constant connecting the 
out-of-plane displacements of - two hydrogen 
atoms para to each other is about half of that 
for the case where they are ortho. From the signs 
of these force constants it is evident that the 
potential energy required for a “‘dishing”’ of the 
hydrogens (i.e., all the hydrogens above the 
carbon ring) is less than that for a “puckering’”’ 
of the hydrogens (hydrogen atoms alternately 


19 Detailed results of this investigation will be published 
separately. 
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TaBLe VII. Calculated vs. observed frequencies. 




















ee Benzene-di o-benzene-d2 m-benzene-d2 p-benzene-d2 sym-benzene-ds 
Freq. iff. 
No. Calc. Obs. (%) Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. Cale. Obs. _— Diff. 
4 693 685s +1.2 633 634 —0.2 691 (691) _ 
5 993 Vv. 959 966 —0.7 914 (914) — 
11 678 (671) +1.0 579 582 —0.5 600 594° +1.0 533 (533) as 
10a 851 852 —-0.1 783 782 +40.1 712 711 +0.1 851 850 +0.1 
851 (850) +0.1 712 (712) _ 
10b 729 738 —1.2 
16a 403 403 0.0 373 374 -—0.3 403 400> +0.8 
403 (404) —0.7 373 (373) -- 
165 391 384 +41.8 360 itr. 
17a 958 — 920 — 958 Vv. 
958 ~9854 ? 920 d. 
176 945 — 868 872° 
. vic-benzene-ds o-benzene-ds m-benzene-d4 p-benzene-d Benzene-ds Benzene-d« 
req. 
No. Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. Calc. Obs. Diff. 
4 615 605 +1.7 600 v. 
$ 918 927 —1.0 812 v.e 
11 534 — 549 ir. 498 (503) —1.0 
10a 712 «712 0.0 737 739 «~-—0.3 712 711 = +0.1 662 662 0.0 662 664 -0.3 
662 (664) —0.3 
10b 759 765 0.8 
loa 373 +374 —0.3 373 375 +0.5 346 v. 346 — 
346 337% +2.7 
166 353 369 —4.3 377° ir. 
17a 920 920 0.0 920 789 Vv. 739 — 
789 790% —0.1 
17b 863 — 918 itr. 








Observed values taken from references 10 and 11. 
)—Used to evaluate force constants. 
ae numbers indicate values which appear in violation of selection 
rules 
v.—Forbidden in both Raman and infra-red. 
ir.—Allowed in infra-red only; hasn't been studied. 


above and below the carbon ring). This is borne 
out by the frequencies of the corresponding vibra- 
tions. The ‘‘dishing’’ frequency (v1) is 671 cm—; 
the ‘‘puckering”’ frequency (v5) is about 1015 
cm (calculated). 

One might deduce from this that attractive 
forces prevail between the hydrogen atoms, 
making it energetically cheaper to displace two 
hydrogen atoms in the same direction out of the 
equilibrium plane than in opposite directions. 
And indeed this is compatible with the geometry 
of the molecule. The van der Waals radius of 
hydrogen is 1.2+0.1A,”! while the closest hy- 
drogen atoms in benzene are approximately 
2.47A apart. Thus they may well be far eriough 
removed so that repulsive forces no longer pre- 
dominate. However, it must be admitted that 
according to this concept a, cannot be more 
positive than both a, and a,—although actually 





°° The frequencies depend also on the reduced masses 
and on interaction with other vibrations which may be in 
the same symmetry class, but these effects alone could 
hardly produce the large difference of 340 cm which is 
observed. 

1 L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, New York, 1940), second edition, p. 189. 


* Reference 10. 

b Reference 11, p. 21, p. 67 (footnote). 

¢ Footnote 19. 

4 Would be hidden by strong i.r. band at 914, 
© Would be hidden by strong i.r. band at 811. 


it is. One must, therefore, conclude that this 
simple hypothesis is untenable. We have no 
adequate explanation to suggest. 

In conclusion it may be pointed out that the 
numerical values for these force constants should 
be applicable to derivatives of benzene. Crawford 
and Brinkley” have already outlined methods 
for such a transfer of force constants between 
aliphatic compounds 


ACKNOWLEDGMENT 


One of us (FAM) wishes to express his appre- 
ciation to the National Research Council for the 
grant of a fellowship. 


APPENDIX 


In most of the isotopic molecules of lower symmetry, 
two blocks of the De, G matrix (Fig. 3) were connected by 
only one interaction term, so that the complete block of 
the G matrix for the isotopic case could be written 











gi) eee ut 0 eee 0 \ 
gi™ SO d coe 0 
» Cf) 
0 eee d 2m4i™*! eee 2uya™*! 
0 0 Lema os gaya 
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with d representing the interaction. If we denote the 
secular equations for the separate blocks (setting d=0) by 


0=A(A)= Da;an-i, (A2) 
, i=0 

0=B(A) = D dari, (A3) 
im0 


and their product by 
0=C(d)=A(A)-B(A) = B cadmtr-t, 
k=0 
(A4) 


k 
Ce= & agdn-i, 
i=0 


then the secular equation with d~0 is simply obtained 
from C(A): wherever the product gmn™gm4i1"*! appears, one 
has only to substitute therefor (gmn”gm4i"*!—d?). 

This is easily seen. For, of all possible minors of the G 
matrix (A1) containing d, those which do not obviously 
vanish will be of the same form as (A1)—i.e., will have 
d's “bridging’’ the principal diagonal—or will be of the 
form : 





gi) ey gn} £m! 
gk gn* gm" X 
, (AS) 
0 0 d 
0 X 
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or of the form 
gu" riches gm" 
gn" eee £m™ 
(A6) 
0 cen an 
S tos & 





The determinants of the form (A6) will clearly not con- 
tain d, while determinants of the form (A5) will correspond 
to vanishing minors of the F matrix; hence no contribution 
from d will arise in either of these cases. So we need con- 
sider only minors of the form (A1). Such a determinant 
may be expanded, by Laplace’s development, in the 
second minors of second-order determinants based on the 
mth and (m+1)th columns. To be non-vanishing, these 
second minors must be obtained by striking out one of 
the first m rows and one of the last 2; this restricts the 
second-order determinants which will actually contribute 
to those of the types 


iA | | a © X d 
lo x|’ d X ad x| 
There is but one of this last type, and it will contribute 
the term (gm"Zm41"*! —d?)—the only way in which d can 
appear in the final secular equation. 
Similar considerations can be applied if more than one 


interaction couples two blocks, but the increased com- 
plexity makes it scarcely worth while. 


? 
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The Infra-Red Spectra of -Benzene-d, and p-Benzene-d, 
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INTRODUCTION 

N a paper by the authors which appears in 
this issue! reference is made to some experi- 
mental work on the infra-red spectrum of p- 
benzene-d2. The purpose of this note is to present 
the data on this compound, together with that 
on p-benzene-d,, in detail sufficient for the above- 
mentioned paper. The work was originally under- 
taken for the purpose of locating v14 and v15 and 
for completing the analysis of the vibrational 
spectrum of benzene. This aspect will be dis- 

cussed in a later report. 


PREPARATION OF SAMPLES 
The deuterium compounds were prepared by 


1F, A. Miller and B. L. Crawford, Jr., J. Chem. Phys. 
14, 282 (1946). 


the method of Langseth and Klit,? which involves 
the continuous formation of a Grignard reagent 
and its decomposition with deuterium chloride. 
The last step is the separation of the deutero- 
benzene from the great excess of diethyl ether 
solvent by distillation. In the case of the p- 
benzene-d, an unknown intermediate fraction 
boiling at 70-75° was obtained in addition to 
the desired p-benzene-d, fraction at 76.3-77.0°.* 
On exposure to air this intermediate fraction 


2 A. Langseth and A. Klit, Kgl. Danske Vid. Sels. Math. 
fys. Medd. 15, No. 13 (1938). 

3 Because these temperatures are uncorrected for atmos- 
pheric pressure (745 mm) and for errors in calibration of 
the thermometer (at least 1.2°), they have only relative 
significance. The boiling point of benzene was observed 
to be 77.0° under the same conditions, although the correct 
value is 80.1°. 
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gradually turned pink. To eliminate any absorp- 
tion bands in the spectrum of the p-benzene-d, 
sample which may be due to this impurity, the 
spectrum of the 70—75° fraction was obtained 
under identical conditions. It was found that all 
absorption bands but one in the spectrum of the 
70-75° fraction were of the same order of in- 
tensity as those in the spectrum of the p-benzene- 
d, sample. They were therefore thought to be 
due to p-benzene-d,, since the 70—75° fraction is 
certainly rich in this. The one exception (1213 
cm~!) was very intense in the 70—75° fraction, 
but was weak in the p-benzene-d, sample. 
Consequently it was considered to originate with 
the impurity and was rejected from the spectrum 
of the deuterobenzene. The diminution of in- 
tensity was so great that there is little likelihood 
of having overlooked any other such inter- 
ferences. 


SPECTROSCOPIC PROCEDURES 


The spectrum of p-benzene-dz was obtained 
on two different spectrometers: (1) one con- 
structed at the University of Minnesota, using 
fluorite (CaF), rocksalt, and KBr prisms. (2) a 











TABLE I. 
p-benzene-d2 p-benzene-ds 
cm~! Int. cm=! nt. 
594 vs 780 w 
657 w 817 vs 
693 m 832 Ss 
766 m 865 Ss 
778 vs 925 vs 
792 m 938 Ss 
815 Ss 1058 m 
829 m 1122 m 
873 vs 1180 w 
924 Ss . 1207 m 
1034 m 1281 w 
1097 vw 1301 w 
1227 w 1341 w 
1325 vw 1396 vw 
1416 s 1447 vs 
1473 vs 1563 m 
1609 w ~1722 vw 
1719 m 1831 vs 
1814 Ss 2283 vs 
1924 Ss ~2320 m 
2101 vw 
2194 w 
2276 vs 
2314 m 
2866 w 
2992 m 
3044 m 
3083 vs 
3279 w 

















TABLE II. 
p-benzene-d2 p-benzene-d« 
Calc. Obs. Calc. Obs. 
Vil 600 594 549 
V16b 360 377 
V7 868 873 918 925 








Perkin-Elmer instrument with rocksalt optics at 
the University of Illinois. Results agreed well, 
so only averaged or preferred values will be 
tabulated here. 

The spectrum of p-benzene-d, was obtained in 
the rocksalt region with the Perkin-Elmer instru- 
ment only. 

All the spectra were measured in the gas phase 
with a vapor pressure of 85 mm of mercury. The 
cell length of the Minnesota spectrometer was 
20 cm; that of the Perkin-Elmer instrument 
10 cm. 


RESULTS AND COMMENT 


The results are tabulated in Table I. 

Several assignments can immediately be made 
to the non-planar vibrations which are infra-red- 
active. These are indicated in Table II.‘ 

There is little doubt that the absorptions of 
873 and 925 cm~ originate from the same vibra- 
tion, for both are very intense and the band 
contours are similar. Even if calculated values 
were not already at hand, one would probably 
want to assign the two bands to »:7 because the 
frequency in p-benzene-d:z is lower than that in 
p-benzene-d,. It is nevertheless noteworthy that 
the value of »:7, can be calculated so precisely. 
This provides convincing support for the assign- 
ment of v;7 advanced by the authors in the 
accompanying paper, and leads one to feel that 
force constants o, 6, and w cannot be far wrong. 

More complete experimental details and an 
analysis of the remainder of the spectra will be 
presented in a later paper. 
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4The calculated values are taken from Table VII of 
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On the Force Constants of the XY. Molecule 


DELIA M. SIMPSON 
Laboratory of Physical Chemistry, Cambridge, England 
February 25, 1946 


N a recent paper Glockler and Jo-Yun Tung! have 

suggested a method whereby the four force constants 
required to describe the behavior of the triatomic molecule 
XY2 (X=M, Y=m) can be determined from its three 
vibrational frequencies. They employ force constants 
C1, C2, C3, C4 defined by the potential field: 


2 V=ai[Ari?+Ar2? ]+co-re[A(2a) 
+2c3-1%o9-A(2a)[Ari+Are ]+2c4-Ari- Are, (1) 


where Ar;, Are are the changes in length of the XY bonds, 
equilibrium value ro, and A(2a) is the change in the apical 
angle. These authors discuss the conditions necessary for 
all four force constants to be real, and suggest that the 
real value of c, should be chosen which makes c; approxi- 
mately the same as the force constant of the corresponding 
diatomic molecule XY. The limits of reality of c, are given 
by the relationship: 


0¢4/0¢2= 0c4/dc3=0, (2) 


which also defines its maximum and minimum values. 
Glockler and Jo-Yun Tung show that if the maximum 
value of cy is chosen, then for molecules such as H2O and 
HS, ¢: is very similar to the force constants of the radicals 
OH and SH. The force constants calculated for these 
triatomic molecules using the three equations for the 
vibrational frequencies and the maximum value of c, are 
used to derive the frequencies of the isotopic molecules 
D.0 and D.S, when good agreement with the observed 
values is obtained. 

The set of force constants defined by Eq. (1) is not 
unique. Kohlrausch,? for example, uses the potential field: 


2V=f[Ari2+Ar2?]+f’-AR?+d-re[A(2a) 
+d’ -ro-Ro[A6:°+AB82"], (3) 


where AR is the change of the YY bond length, equilibrium 
value Ro, and AB; and Afz are the changes in the basal 
angles. Both c; and f define the behavior of the XY bond, 
and cz and d that of the YXY angle, but they are not 
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identical for: 
ai=f+f’ sin? a+2d’ cos? a, 
c2=d+d’'/2+f"’ cos? a. 


The method of Glockler and Jo-Yun Tung may be 
criticized on the following general grounds. 

(i) In a polyatomic molecule the force constant of a 
particular bond will not necessarily be the same as that 
of the corresponding diatomic molecule, even if the number 
of electrons involved is identical, for the bond environment 
must be different in the two cases. 

(ii) The force constant of a bond in a polyatomic 
molecule is not uniquely defined, as is shown by (4) above. 
There is thus no fundamental reason for preferring to make 
c; rather than f equal to the force constant of the diatomic 
molecule. 

(iii) The relationship (2) gives other relationships if 
different sets of force constants are used, for example for 


(3): 


(4) 


ad’ /af’ =tan? a/2. (5) 


Neither (2) nor (5) seems to have any obvious physical 
justification. 

(iv) Since c, is essentially an interaction constant it is 
surely extraordinary that it should be given its largest 
possible value. 

(v) For the particular case m= WM the relationship (2) 
gives the remarkable condition that f’=f whatever the 
value of 2a. This conclusion is clearly ridiculous unless 2a 
is approximately 7/3. 

Although Glockler and Jo-Yun Tung have obtained 
good results for a number of molecules, it would, therefore, 
appear that their method is not one of general applicability 
and that the information obtained by using it should not 
be considered entirely reliable. 

1G. Glockler and Jo-Yun Tung, J. Chem. Phys. 13, 388 (1945). 


2K. W. F. Kohlrausch, Der Smekal-Raman-Effekt (Verlagsbuch- 
handlung, Julius Springer, Berlin-1938), p. 65. 





Comment on “On the Force Constants of the 
XY, Molecule’’ 


GEO. GLOCKLER AND Jo-YUN TUNG 
University of Iowa, Iowa City, Iowa 
March 1, 1946 


E appreciate the fact that the Editor has permitted 
us to see the above note by Miss Simpson on 
“Force Constants of the XY2 Molecule.” 

We would like to point out that the idea of the paper by 
Glockler and Tung was first to indicate a simple representa- 
tion of all the real values of the force constants of the XY: 
molecule and secondly to make an attempt to find a simple 
empirical method which might be useful in determining a 
unique set of force constants from all the possible ones. 
So far the empirical condition studied was shown to hold 
with seven molecules. We ourselves pointed out. that 
similar considerations would not apply to the COS mole- 
cule. Since we are continuing our calculations we will 
refrain at this time from taking valuable space and will 
answer the criticism of Miss Simpson at a later date. 


294 





1946 


iy be 


of a 
; that 
imber 
ment 


tomic 
bove. 
make 
tomic 


ips if 
le for 


(5) 
ysical 


- it is 
irgest 


p (2) 
r the 
ss 2a 


ained 
fore, 
bility 
1 not 


sbuch- 


the 


itted 
n on 


er by 
enta- 
XY; 
mple 
ing a 
ones. 
hold 
that 
nole- 

will 

will 





LETTERS TO 


The Photoelectric Effect in Metal-Ammonia 
Solutions 


RICHARD A. OGG, JR. 
Department of Chemistry, Stanford University, California 
March 14, 1946 


HE previous studies! of inner photoelectric effect 
(i.e., of photo-conductivity) of metal-ammonia solu- 
tions have been extended to include the external photo- 
electric effect—i.e., the emission of electrons into the 
vapor phase. Extremely dilute (some 10-* to 10~* molar) 
sodium solutions were employed, in a special photo-cell 
maintained just above the freezing point of ammonia. 
The light source was a tungsten filament projection lamp, 
in conjunction with a carefully tested set of light filters. 
The quantum efficiency of photoelectric emission appeared 
to be similar in magnitude to that previously observed? 
for relatively concentrated solutions, but the spectral 
dependence proved entirely different. Whereas the long 
wave-length threshold? varies from some 7700A for a 
saturated solution to some 8500-8700A for a 4X 10-? molar 
solution, the highly dilute systems here studied displayed 
their maximum sensitivity in the near infra-red. Exact 
determination of the corresponding threshold requires 
more precise work with an infra-red monochromator, but 
with the filters employed it is estimated to be in the region 
around 15,000A. The extremely dilute solutions in question 
would thus appear to possess a work function of the order 
of one-half of that of the concentrated solutions. 

The extraordinary difference in photoelectric behavior 
of extremely dilute and relatively concentrated sodium- 
ammonia solutions is best explained by the previously 
proposed! * equilibrium between trapped electron pairs 
and trapped single electrons, and by the striking difference 
in absorption spectrum’ between these two components. 
At the very great dilutions employed in the present studies, 
the proportion of trapped pairs is negligible, and the 
estimated photo-threshold at A~15,000A corresponds to 
the work function of the process 


e(trapped in NH3)+/v—e(gas), @~0.8 volt. (1) 


However, measurements of magnetic susceptibility* ® indi- 
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cate that even the most dilute solution previously studied? 
is characterized by a great predominance of trapped 
electron pairs over trapped single electrons. Further, in 
the red and near infra-red regions, the absorption coeffi- 
cient? of trapped pairs greatly exceeds that of trapped 
single electrons. Consequently, in such solutions the above 
process (1) is quantitatively negligible, and the observed 
photo-threshold at \<8700A corresponds to the work 
function of the process 


é2(trapped in NH3)+hy—> (2) 
e(trapped in NH;)+e(gas), #>1.42 volts. 


The possibility of precise experimental determination of 
the work function for process (1) is especially interesting 
as affording an upper limit to the “solvation energy” of 
single electrons in liquid ammonia solution. (The excess 
over the “true” energy results from application of the 
Franck-Condon principle to the non-adiabatic removal of 
the electron from the solvent cavity). This allows a 
strictly experimental approach to the long standing theo- 
retical problem® of separation of sums of solvation energies 
of positive and negative ions into the respective contribu- 
tions of the individual ionic species. 

The difference in work functions between processes (2) 
and (1) provides an upper limit (for reasons similar to 
those outlined above) to the dissociation energy of the 
trapped electron pair into two electrons trapped in separate 
cavities. The provisional value obtained here seems not 
unreasonably greater than that roughly estimated from 
the temperature dependence of the magnetic suscepti- 
bilities,“5 and from current calorimetric studies which 
show dilute (some 107 to 10-! molar) sodium-ammonia 
solutions to display a definitely negative heat of dilution— 
in contrast to the positive-effect expected for simple 
electrolytes, and experimentally observed with ammonia 
solutions of sodium bromide.’ 

1R. A. Ogg, Jr., J. Am. Chem. Soc. 68, 155 (1946). 

2 J. Hasing, Ann. d. Physik [5] 37, 509 (1940). 

3R. A. Ogg, Jr., J. Chem. Phys. 14, 115 (1946). 

4E. Huster, Ann. d. Physik 33, 477 (1938). 

5S, Freed and N. Sugarman, J. Chem. Phys. 11, 354 (1943). 

6 For literature references, see O. K. Rice, Electronic Structure and 


Chemical Binding (McGraw-Hill Book Company, Inc., New York, 1940). 
7C. A. Kraus and F. C. Schmidt, J. Am. Chem. Soc. 56, 2297 (1934). 





